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Abstract 

The spectral method for building first integrals of ordinary linear differential systems is 
elaborated. Using this method, we obtain bases of first integrals for linear differential systems 
with constant coefficients, for linear nonautonomous differential systems integrable in closed 
£SJ \ form (algebraic reducible systems, triangular systems, the Lappo-Danilevskii systems), and 

for reducible ordinary differential systems with respect to various transformation groups. 
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Introduction 

One of the most important problems of the general theory of differential systems is the 
problem of finding first integrals. In 1878 the French mathematician J.G. Darboux showed 
how general integrals of the first-order ordinary differential equations possessing sufficient in- 
variant curves are constructed [1]. His investigation gave the classical problem (the Darboux 
problem) about building of first integrals by known partial integrals. The review of the litera- 
ture and the current situation of the theory of integrals are given in the monographies [2 - 5] . 

Linear differential systems are of interest for mathematicians both per se and as a tool 
for studying nonlinear differential equations by means of the method of linearization. Note 
also that systems of linear differential equations play a broad and fundamental role in electri- 
cal, mechanical, chemical and aerospace engineering, communications, and signal processing. 
There exist many good works on linear differential systems (see, for example, [6 - 24]). 

In this paper we study the Darboux problem of the existence of first integrals for main 
classes of linear ordinary differential systems. Using the method of partial integrals for poly- 
nomial differential systems [3, pp. 187 - 226; 28 - 31], we obtain the spectral method for 
building first integrals of linear differential systems [25 - 27]. 

At a later time, the spectral method has been applied to linear multidimensional differen- 
tial systems [32 - 38] and to nonlinear Jacobi's differential systems [39 - 44]. 

The material of this paper is made on the base of our articles [25 - 27; 36; 37; 45; 46]. 

To avoid ambiguity, we stipulate the following notions (using the article [47]). 

Consider an ordinary differential system of n-th order 

^ = X(t,x), (0.1) 

where x G R n , tGK, the column vector ^ = colon (-7^ , . . . , — irO , the coordinates of the 

at V at at J 

vector function X(t, x) = colon (X^t, x), . . . , X n (t, x)) are the continuously differentiable on 

a domain G C R n+1 scalar functions X i : G — > R, i = 1, . . . , n. 

We recall that by domain we mean open arcwise connected set. 

A continuously differentiable scalar function F: G' — > R is called a first integral on a 
domain G' C G of the ordinary differential system (0.1) if 

XF(t,x) = for all (t,x) G G', 

n 

where the linear differential operator X(t,x) = d t + ^ X i (t,x)d x for all (t,x) G G is the 

i=l 

operator of differentiation by virtue of system (0.1). 

A continuously differentiable scalar function F: G' — > R is a first integral on a domain 
G' C G of the ordinary differential system (0.1) if and only if the function F: G' — > R is 
constant along any solution x: t — > x(t) for all t G J of system (0.1), where x: t — > x(t) is 
such that (t,x(t)) G G' for all t G J C R, i.e., F(t,x(t)) = C for all t G J, C = const. 

A continuously differentiable scalar function w: G' — > R is said to be a partial integral 
on a domain G'cG of the ordinary differential system (0.1) if 

Xw(t,x) = <f>(t,x) for all (t, x) G G', 

where <I>: G' — > R is a scalar function such that &(t,x)\ = for all (t,x) G G'. 

A continuously differentiable scalar function w : G 1 — > R is a partial integral on a domain 
G' C G of the ordinary differential system (0.1) if and only if the function w: G' — >■ R 
vanishes identically along any solution x: t — > x(t) for all t G J CK of system (0.1), where 
x: t — > x{t) is such that (t,x(t)) G G' for all t G J, i.e., w(t,x(t)) = for all t G J. 
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A set of the functionally independent on a domain G'cG first integrals F^. G' — > R, 
I = 1, . . . , k, of system (0.1) is called a basis of first integrals (or integral basis) on the domain 
G' of system (0.1) if for any first integral G' — > R of system (0.1), we have 

V(t,x) = $(F 1 (t,x),...,F k (t,x)) for all (t,x)€G', 

where is some continuously differentiable function on the range of the vector function 
F: (t,x) — > (F x (t, x), . . . , F k (t, x)) for all (t,x) G G'. The number k is said to be the 
dimension of basis of first integrals on the domain G' for system (0.1). 

The ordinary differential system (0.1) on a neighbourhood of any point of the domain G 
has a basis of first integrals of dimension n [48, pp. 175 - 177; 49, pp. 256 - 263]. 

If the ordinary differential system (0.1) is autonomous (the vector function X : x — > X{x) 
for all i £ I C R™ not depends on the independent variable t), then the ordinary differential 
system (0.1) in an n-dimensional integral basis has n — 1 functionally independent on the 
domain X autonomous first integrals [50, pp. 161 - 169]. 

The paper is organized as follows. 

In Section 1 the spectral method for building first integrals for linear homogeneous and 
nonhomogeneous differential systems with real constant coefficients is developed. Here we pay 
special attention to the construction of bases of real first integrals for these systems. 

The Euler method [51, pp. 350 - 360; 52, pp. 93 - 101] and the matrix method [51, pp. 
329 - 349; 53, pp. 166 - 172] are the main methods for finding solutions of linear homogeneous 
differential systems with constant coefficients. For building first integrals we know the method 
of integrable combinations [48, pp. 171 - 173] and the N.P. Erugin — N.A. Zboichik method 
[54, pp. 464 - 469; 55]. These approaches show only ways of constructing first integrals for 
linear differential systems, but they do not build first integrals in explicit form. 

Using the method of partial integrals [3, pp. 187 - 226; 28 - 31] for polynomial differential 
systems, we obtain the spectral method of building first integrals in explicit form for linear 
differential systems with constant coefficients [26; 27; 45]. 

Note also that the results of Subsection 1.1.2 are co-ordinated with the articles [20 - 24]. 
The approach of finding autonomous first integrals with using the method of Jordan cano- 
nical forms was developed by M. Falconi and J. Llibre [20]. The existence of common first 
integrals for two coupled third-order linear differential systems (which satisfy the Frobenius 
compatibility condition) is discussed in [21]. The problem of building rational first integrals for 
linear systems was considered by the French mathematician J. A. Weil in the paper [22]. The 
class of linear autonomous differential systems with no rational first integrals is determined 
by A. Nowicki [23]. First integrals are constructed using integrating factors in [24]. 

In Section 2 we investigate the problem of the existence of first integrals for linear nonau- 
tonomous differential systems integrable in closed form. 

In the general theory of differential systems we know some results about reducibility of 
differential systems to some special forms. For example, any linear nonautonomous differential 
system is reducible to an algebraic reducible system or to a triangular system [56; 8, p. 33]. 

Three classes of linear nonautonomous differential systems integrable in closed form (al- 
gebraic reducible systems, triangular systems, the Lappo-Danilevsky systems) are considered. 
The regular method of building first integrals for these systems is elaborated [36; 46]. Also, 
the N.P. Erugin problem of the existence of autonomous first integrals [54, p. 469] for an 
nonautonomous homogeneous Lappo-Danilevskii differential system is solved [25]. 

In Section 3 we consider reducible linear nonautonomous differential systems with respect 
to various transformation groups (periodic, polynomial, orthogonal, Liapunov, exponential, 
and other groups) to the linear differential systems with constant coefficients. 

The spectral method for building first integrals of reducible systems is elaborated. 

The notion of reducible systems has been entered by A.M. Lyapunov [57]. The develop- 
ment of the theory of reducible systems is associated with the name of N.P. Erugin [58]. 

In addition, in this article some examples are given to illustrate the obtained results. 
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1. Integrals of ordinary linear differential system 
with constant coefficients 

1.1. Linear homogeneous differential system 

Consider an ordinary linear homogeneous differential system with constant coefficients 

* = (i.i) 

where x = colon(x 1; . . . , x n ) G M. n and A = ||a^|| is an square constant matrix of order n 
with entries s^el, £ = 1, . . . , n, j = 1, . . . , n. The differential system (1.1) is induced the 
autonomous linear differential operator of first order 

n 

5=1 

where the vectors = (a^, . . . , a^ n ), £ = 1, . . . , n. 

The differential system (1.1) on a domain G C W l+1 has a basis of first integrals of 
dimension n. Moreover, the differential system (1.1) has also n— 1 autonomous functionally 
independent first integrals (autonomous integral basis) on a domain X C W 1 . 

The aim of Section 1 is to build these bases of first integrals for system (1.1). 

1.1.1. Linear partial integral 

A complex-valued linear homogeneous function 

n 

p: x ^J2 U (, X Z for all a; eM" (i^ G C, £ = 1, . . . , n) (1.2) 
5=1 

is a partial integral of system (1.1) if and only if 

2lp(x) = Xp(x) for all x G W 1 , A G C. (1.3) 

The identity (1.3) is equivalent to the linear homogeneous system 

(B-XE)u = 0, (1.4) 

where the column vector v = colon(z/ 1 , . . . , v n ) G C n , E is the n x n identity matrix, and 
the matrix B is the transpose of the matrix A. 

The linear homogeneous system (1.4) has a nontrivial solution if and only if 

det(5 - \E) = 0. (1.5) 

We shall say that the equation (1.5) is the integral characteristic equation of system (1.1), 
and its roots are integral characteristic roots of system (1.1). Besides, a solution v of the 
linear system (1.4) is an eigenvector of the matrix B corresponding to the eigenvalue A. 

Thus we have proved the following statement. 

Lemma 1.1. The linear homogeneous function (1.2) is a partial integral of the differential 
system (1.1) if and only if the vector v G C n is an eigenvector of the matrix B. 

Therefore a linear partial integral of the linear homogeneous differential system (1.1) is 
generated by an eigenvector of the matrix B. 

Lemma 1.1 is base for ours spectral method of building first integrals of system (1.1). 
In addition, first integrals of the differential system (1.1) are building by eigenvectors and 
eigenvalues of the matrix B with using orders of elementary divisors. 

1.1.2. Autonomous first integrals 

Case of real eigenvalues. If the matrix B has two linearly independent real eigenvectors, 
then we can find an autonomous first integral of the differential system (1.1) by using following 
assertions (Theorem 1.1, Corollaries 1.1 and 1.2). 
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Theorem 1.1. Suppose u l and v 2 are real eigenvectors of the matrix B corresponding 
to the distinct eigenvalues X 1 and X 2 (X 1 / A 2 ), respectively. Then the ordinary linear 
homogeneous differential system (1.1) has the autonomous first integral 

F: x -> \v l x\ hl \v 2 x\ h2 forallx^X, (1.6) 

where X is a domain of the domain of function DF C W 1 , the numbers h l and h 2 are a 
real solution to the equation X 1 h 1 + X 2 h 2 = with |/t 1 | + \h 2 \ ^ 0. 
Proof. By Lemma 1.1, the linear functions 

p k : x — > v k x for all x G R n , k = 1, k = 2, 

are partial integrals of system (1.1). Hence, 

SlzAr = A fc v k x for all x G W 1 , k = 1, k = 2. 

Using these identities, we have 

2LF(:r) = l^ 1 ^!^ 1 1 |z^ 2 x| ft2 1 (/i 1 sgn(z^ 1 x) \v 2 x\ %v l x + h 2 sgn(z^ 2 x) | z^ 1 ^ | 2tz/ 2 x) = 

= (A 1 /i 1 + A 2 /i 2 ) F(x) for all iflc DF. 

If /i x and /i 2 are real numbers such that X 1 h 1 + X 2 h 2 = with j/j-J + \h 2 \ / 0, then 

the scalar function (1.6) is an autonomous first integral of system (1.1). ■ 

Corollary 1.1. If u is a real eigenvector of the matrix B corresponding to the eigenvalue 
A = 0, then the linear homogeneous function 

F: x^vx for all x £ R n (1.7) 

is an autonomous first integral of the ordinary linear homogeneous differential system (1.1). 

Indeed, suppose v is a real eigenvector of the matrix B corresponding to the eigenvalue 
A = 0. Then, from the identity (1.3) it follows that 2lz/x = for all x G W 1 . This yields 
that the function (1.7) is an autonomous first integral on the space W 1 of system (1.1). I 

Corollary 1.2. Let A / be an eigenvalue of the matrix B corresponding to two 
real linearly independent eigenvectors v 1 and v 2 . Then the ordinary linear homogeneous 
differential system (1.1) has the autonomous first integral 

F:x^f V -^- forallxeX, (1.8) 
u z x 

where X is a domain from the set {x: v 2 x / 0} C M". 

Let us remark that the first integrals (1.7) and (1.8) are algebraic. Note also that alge- 
braicity of the first integral (1.6) depends on the numbers h-y and h 2 . For example, if the 
numbers h 1 and h 2 are rational, then the first integral (1.6) is algebraic. But this condition 
isn't necessary for algebraicity of the first integral (1.6). At the same time we have 

Property 1.1 (sufficient condition for algebraicity of basis of autonomous first integrals). 
If all eigenvalues of the matrix B are simple and rational, then the differential system (1.1) 
has a basis of autonomous algebraic first integrals. 

Example 1.1. Consider the fourth-order ordinary linear differential system 



dx-^ dx<2 a n 

—j-^ — x 1 — 2x 2 — x 4 , —j^- = — x± + 4x 2 — x 3 + 2x 4 , 

dx^ dx ^ 
dt 2 6 4 dt 



(1.9) 
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We claim that the matrix 



B 



1-10 2 
-2 4 2-4 

0-11 2 

- 1 2 1 - 2 



has the eigenvalues X ± = 0, A 2 = A 3 = 1, and A 4 = 2. Indeed, the characteristic equation 
1-A -10 2 



2 4 - A 2 -4 

-1 1-A 2 

1 2 1 -2-A 







A(A-l) 2 (A-2) = 0. 



The rank of the matrix B — X 2 E is equal 2. Therefore the double eigenvalue A 2 = 1 of 
the matrix B has k = 4 — 2 = 2 simple elementary divisors A — 1 and A — 1 . 

The matrix B has four simple elementary divisors A, A — 1, A — 1, and A — 2. 
The linear homogeneous system 

v \ ~ v 2 + 2z/ 4 = °' 
- 2v x + 4z/ 2 + 2i/ 3 - 4z/ 4 = 0, 

- v 2 + ^3 + 2z/ 4 = °> 

— ^ 1 + 2v 2 + ^3 — 2z/ 4 = 



(B — A 1 £^) colon (z/ l5 . . . , i/ 4 ) = 



4- 



'4> 



V "3 



I/„. 



Hence i/ 1 = ( — 1, 1, — 1,1) is a real eigenvector corresponding to the eigenvalue A x = 
of the matrix B. 

The scalar function (by Corollary 1.1) 

F 1 : x — > — x 1 + x 2 — x 3 + x 4 for all x£l 4 

is an autonomous linear first integral of the ordinary differential system (1.9). 
The linear homogeneous system 

(B — X 2 E) colon (i^, . . . , i/ 4 ) = <^=^ 



"1 = "3 + v A-> 
v 2 = 2z/ 4 . 



- v 2 + 2z/ 4 = 0, 
- 2v x + 3z> 2 + 2^3 - 4z/ 4 = 0, 

- v 2 + 2z/ 4 = 0, 
— v 1 + 2f 2 + 1/3 — 3z/ 4 = 

Hence v 2 = (2, 2, 1, 1) and u 3 = (1, 0, 1, 0) are two linearly independent real eigenvectors 
corresponding to the double eigenvalue A 2 = 1 of the matrix B. 
The scalar function (by Corollary 1.2) 



^23 : x 



x-^ ~\~ x^ 



for all x £ X 1; 



where a domain X 1 C {x : x 1 + x 3 7^ 0}, is an autonomous first integral of system (1.9). 
The linear homogeneous system 

- v l - v 2 + 2u 4 = 0, 



(B — X 4 E) colon(i^ 1 , . . . , v 4 ) = 



- 2v x + 2v 2 + 2u 3 - 4v 4 = 0, 

- v 2 ~ v 3 + 2 ^4 = °) 
— v 1 + 2v 2 + f 3 — 4^ 4 = 



"1 = ^3> 
v 2 = 2v 4 . 



Hence v A = (0,2,0,1) is a real eigenvector corresponding to the eigenvalue A 4 = 2 of 
the matrix B. 

6 



V.N. Gorbuzov, A.F. Pranevich 



First integrals of ordinary linear differential systems 



Using the linearly independent real eigenvectors v 2 and u A of the matrix B, we can 
build (by Theorem 1.1) an autonomous first integral of the differential system (1.9). Since 
the equation h 2 + 2h 4 = 0, we have, for example, h 2 = 2, h 4 = — 1. The scalar function 

(2x 1 + 2x 2 + x 3 + x 4 ) 2 , „ „ 
F 24 : x -)• - — — — for all x G X 2 , 

where a domain X 2 C {x: 2x 2 + x 4 / 0}, is an autonomous first integral of system (1.9). 

The functionally independent first integrals F 1 , F 23 , and F 24 are an autonomous integral 
basis of system (1.9) on any domain X from the set {x: x l + x 3 / 0, 2x 2 + x 4 / 0} C R 4 . 

Case of complex eigenvalues. Let the function (1.2) be a complex- valued partial integral 
of system (1.1). Then, from the identity (1.3) it follows that 

2lRep(x) = \ Rep(x) - A Imp(x) for all x G R n , 

(1.10) 

2llmp(x) = ARep(x) + A Imp(x) for all x G R n , 

where the real numbers A = Re A, A = ImA. Thus we have the following criteria for the 
existence of a complex- valued partial integral of system (1.1). 

Lemma 1.2. The function (1.2) is a complex-valued partial integral of system (1.1) if and 
only if the system of identities (1.10) holds. 

Using Lemma 1.2, we may establish the following propositions. 

Property 1.2. If the system (1.1) has a complex-valued partial integral (1.2), then the 
complex conjugate function p is a complex-valued partial integral of system (1.1) and 

Slp(x) = A p(x) for all x G K™, 

where A is the conjugate of the complex number A from the identity (1.3). 
Proof. Using the system of identities (1.10), we get 

* , — * 

2lp(x) = 2lRep(x) — z2llmp(x) = ARep(x) — A Imp(x) — i (XKep(x) + A Imp(x)) = 

= (A - i\)(Rep(x) - i1mp(x)) = \p(x) for all xeR n .M 

Let K be either the real number field R or the complex number field C. 

Property 1.3. The product u 1 u 2 of the polynomials u 1 : R n — > K and u 2 : R™ — >■ K is 
a polynomial partial integral of system (1.1) «/ and only if the polynomials u l and u 2 are 
polynomial partial integrals of system (1.1). 

Proof. From definition of partial integral and 

2l(n 1 (x)n 2 (x)) = u 2 (x)^lu 1 (x) + u 1 (x)%Lu 2 (x) for all x G R n , 

we get the assertion of Property 1.3. I 

Property 1.4. The system (1.1) has the complex-valued polynomial partial integral (1.2) 
if and only if the real polynomial 

P: x -)• Re 2 p(x) + Im 2 p(x) /or all x G R n 

is a partial integral of system (1.1). Moreover, the following identity holds 

2lP(x) = 2AP(x) /or a// x G R n , 

where the number A = A + i A (A = Re A, A = ImA) is <?u>en by the identity (1.3). 
Proof. Using Properties 1.2 and 1.3, and the identity 

p(x)p(x) = Re 2 p(x) +Im 2 p(x) for all x G R n , 
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we have the polynomial P is a real partial integral of system (1.1) if and only if the linear 
function (1.2) is a complex- valued partial integral of system (1.1). 
By the system of identities (1.10), it follows that 

2t(Re 2 p(x) + Im 2 p(x)) = 2Rep(x) 2lRep(x) + 2lmp(x) %1mp(x) = 
= 2Rep(x) (A Rep(x) — A Imp(x)) + 2Imp(x) (A Rep(x) + A Imp(x)) = 

= 2 A (Re 2 p(x) + Im 2 p(x)) = 2 A P(x) for all 

Property 1.5. Let the function (1.2) 6e a complex-valued partial integral of system (1.1). 
T/ien t/ie Lie derivative of the function 

Im j?(x) 

99: x — > arctan — — ^— y /or a// iflcji: Rep(x) 7^ 0} 

6y virtue of system (1.1) is egua/ to 

2l<£>(x) = A /or a// x£l, 

where the number A = A + i A (A = Re A, A = Im A) is given by the identity (1.3). 
Proof. Using the identities (1.10), we obtain 

lmp(x) 1 Kep(x) 2llmp(x) — lmp(x) 2lRep(x) 



2l</?(x) = 21 arctan 



Rep(x) Im 2 p(a;) Re 2 p(x) 



Re 2 p(x) 

Rep(x)(ARep(x) — A Imp(x)) — Imm(x) (A Rep(x) + A Imp(x)) 



A for all x £ X. 



Re 2 p(x) + Im 2 p(x) 

If the matrix B has an imaginary eigenvalue, then we can build an autonomous first 
integral of system (1.1) by using following assertions (Theorems 1.2, 1.3, and 1.4). 

* ~_ * ^ . 

Theorem 1.2. Suppose v = v + vi [Rev = v, Imv = v) is an eigenvector of the matrix 

5 corresponding to the imaginary eigenvalue A = A + Ai (Re A = A, ImA = A / 0). T/ien 
t/ie ordinary differential system (1.1) /ias i/ie autonomous first integral 

* 

* I \ vx \ 

F: x -> ((vx) 2 + (?x) 2 ) exp - 2 £ arctan— for all x G X, (1.11) 

V A z/x / 

where X is a domain from the set jx: 

Proof. Taking into account Properties 1.4 and 1.5, we get 

* ~ * ~ 
/ \ ^/j^ \ # * I \ vx 

2LF(x) = exp(-2^ arctan— 2l((i/x) 2 + (?x) 2 ) + ((^x) 2 + (?x) 2 ) 2lexpl-2 0= arctan — 
V A vx) V A vx 

* * 

(* \ ~ \ * /A z^x \ 

2A - 2 ^ Aj (( vx) 2 + ( vx) 2 ) exp^ - 2 Z arctan — J = for all x G X. 

This implies that the function (1.11) is an autonomous first integral of system (1.1). ■ 
Transcendency of the first integral (1.11) of systems (1.1) depends on the imaginary eigen- 
value A = A + Ai (if A = 0, then the first integral (1.11) of system (1.1) is algebraic). 
The proof of Theorems 1.3 and 1.4 is similar to that one in Theorem 1.2. 
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Theorem 1.3. Let v 1 = v 1 +u 1 i (Re^ 1 = v l , \m.u l = u l ) be an eigenvector of the matrix 
B corresponding to the imaginary eigenvalue A 1 = + (ReA 1 =Ai, ImA x = A 1 7^ 0), v 
be a real eigenvector of the matrix B corresponding to the eigenvalue A 2 / 0. Then the 
ordinary differential system (1.1) has the autonomous first integral 

„ 9 / A 2 v 1 x \ 

F : x — > v x exp — — arctan - — lor all x G X, 

V X 1 fax J 

where X is a domain from the set jx: v^x 7= 0}. 

Theorem 1.4. Let u 1 = v l +V 1 i and v 2 = v 2 + z7 2 i (Re v T = v T , Im v T = v T , r = 1, 2) 

6e fauo linearly independent eigenvectors of the matrix B corresponding to the imaginary 

* * — - * . — • 

eigenvalues A 1 = Ai + A 1 i and A 2 = A2 + A 2 i (Re A T = A T; Im A r = A T , r = 1, 2) X ± 7^ A 2 . 

T/ien i/ie ordinary differential system (1.1) /ias i/ie autonomous first integral 

~ 2 ~ 1 

' V X ~ V X 

F : x — > Ai arctan - — — A2 arctan - — for 
9 ~ 1 

z/ z x z^x 

where X is a domain from the set {x: z^x 7^ 0, z^ 2 x 7^ 0}. 
Example 1.2. The autonomous differential system 

^- = 2x 1 +x 2 , = x 1 + 3x 2 - x 3 , ^ = - x 1 + 2x 2 + 3x 3 (1.12) 

has the eigenvalues X 1 = 3 + i and A 2 = 2 corresponding to the eigenvectors v l = (1, i, —1) 
and v 2 = (3,-1,-1), respectively. The first integrals (by Theorem 1.2) 

F 1 : x — > ((x 1 — x 3 ) 2 + x 2 ) exp ^ — 6 arctan — — — ^ for all x G X 

and (by Theorem 1.3) 

F 2 : x — >■ (3x 1 — x 2 — x 3 ) exp^ — 2 arctan — ^ for all x G X 

are an autonomous integral basis of system (1.12) on any domain X C {x: x 1 — x 3 / 0}. 
Example 1.3. The autonomous differential system 

(1.13) 



— - 1 - = — 3xi + x 9 + 4xq + 2x 4 , — - = 8x1 — 3x 9 — 2xo + 6x 4 , 
dt dt 



—— = — 9x 1 + 3x 2 + 4x 3 — 4x 4 , — — = 6x 1 — 3x 2 — 4x 3 + 2x 4 

tit (Jbb 

has the eigenvalues A x =i, X 2 = 2i corresponding to the eigenvectors u 1 = (l — i, — l + 2i,2i,2), 
v 2 = (i, — l,i, 1 + 2i). The functionally independent first integrals (by Theorem 1.2) 

F l : x ->■ (x x - x 2 + 2x 4 ) 2 + ( - x 1 + 2x 2 + 2x 3 ) 2 for all x G M 4 , 

F 2 : x -)• ( - x 2 + x 4 ) 2 + (x 1 + x 3 + 2x 4 ) 2 for all x G M 4 , 
and (by Theorem 1.4) 

Xi ~\- Xq 2X/i Xi "I - 2Xo ~t~ 2Xq 

F 3 : x ->■ arctan — - - 2 arctan ± - for all x G X 

x<2 ~\~ x ^ X-^ x 2 — \~ 2x 4 

are an autonomous integral basis of the ordinary differential system (1.13) on any domain X 
from the set {x : x 1 — x 2 + 2x 4 / 0, x 2 — x 4 7^ 0} C R 4 . 
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Case of multiple elementary divisors. 

Definition 1.1. Let u° be an eigenvector of the matrix B corresponding to the eigen- 
value X with the elementary divisor of multiplicity m. A non-zero vector u k G C n is called 
a generalized eigenvector of order k if the vector u k satisfying 

(B - XE) u k = his*' 1 , k = 1, . . . ,m - 1. (1.14) 

In this case we can build first integrals of system (1.1) by using following assertions. 

Theorem 1.5. Let X be the eigenvalue of the matrix B with elementary divisor of 
multiplicity m (to ^ 2) corresponding to the real eigenvector v° and to the real 1-st order 
generalized eigenvector v 1 . Then the system (1.1) has the autonomous first integral 

F: x — > u°x expf — A ) for all x G X, Icfi: u°x / 0}. 
\ u v x J 

Proof. Using the equalities (1.14) under the condition k = 1, we get 

21 v x x = X^x + u°x for all x G W l . 

Then, the Lie derivative of the function F by virtue of system (1.1) is 

~, n /\ AXv 1 x + i> x)v°x — Xv°xv 1 x \ ^, . „ _ 

VLF(x )= A -A- )F(x) = for all x G X. ■ 

v (v u xy j 

From Theorem 1.5, we have the following 

Corollary 1.3. Let X = X + Xi (Re A = A, ImA = A / 0) be the complex eigenvalue 
of the matrix B with elementary divisor of multiplicity m (m ^ 2) corresponding to the 
eigenvector u° = u° + v° i (Rev = t>°, Imv° = u°) and to the 1-st order generalized eigen- 
vector v 1 = u 1 +u 1 i (Re 2^= is 1 , lmv 1 = u 1 ). Then the ordinary differential system (1.1) on 
a domain X C {x: v^x / 0} has the autonomous first integrals 

F 1 :x-,U>xf+ (v°x) 2 ) expf - 2 A°(s) - Afls) \ 

and 

u°x Xa(x) + Xp(x) 
Fr, : x — > arctan - — — — for all x G X, 

>x (>x) 2 + {v*xf 

where the polynomials a : x — > v°x v 1 x+v °xu 1 x, (3: x — > i/ xv 1 x — v x v x x for all x G R n . 
The proof of Theorem 1.6 is similar to that one in Theorem 1.5. 

Theorem 1.6. Suppose v° and v 1 are a real eigenvector and a real 1-th order generalized 
eigenvector of the matrix B corresponding to the eigenvalue X 1 = with elementary divisor 
of multiplicity m ^ 2, and v 2 is a real eigenvector of the matrix B corresponding to the 
eigenvalue A2. Then the system (1.1) on a domain X has the autonomous first integral 

F: x — > v 2 x expf — A 2 —^r- ) for all x G X, X C {x: u°x 7^ 0}. 
\ v u x J 

If A 2 = 0, then from Theorem 1.6, we have Corollary 1.1. 

Example 1.4. The autonomous system of ordinary linear differential equations 

^ = \ %1 - 5x 2 + 2x 3l = 5x 1 -7x 2 + 3x 3 , = 6x 1 - 9x 2 + Ax 3 (1.15) 

has the eigenvalue A x = with elementary divisor A 2 of multiplicity 2 corresponding to the 
eigenvector u 1 = (1, —2, 1) and to the 1-st order generalized eigenvector v 2 = (0, —1,1), and 
the simple eigenvalue A 3 = 1 with elementary divisor A — 1 corresponding to the eigenvector 
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v 3 = (3, — 3, 1). The scalar functions (by Theorem 1.5) 

F 1 : x — > x 1 — 2x 2 + x 3 for all x G M 3 

and (by Theorem 1.6) 

F 2 : x — > (3x l — 3x 2 + x 3 ) exp — — for all x G X 

are an autonomous integral basis of system (1.15) on a domain X C {x: x 1 — 2x 2 + x 3 / 0}. 
From Theorem 1.6, we obtain 

Corollary 1.4. Suppose v° and v 1 are a real eigenvector and a real 1-th order generalized 
eigenvector of the matrix B corresponding to the eigenvalue A 1 = with elementary divisor 
of multiplicity m ^ 2, and v 2 = v 2 + z? 2 i (Rez/ 2 = is 2 , Imi/ 2 = v 2 ) is an eigenvector of the 
matrix B corresponding to the complex eigenvalue \ 2 = \ 2 + X 2 i (ReA 2 = \ 2 , ImA 2 = A 2 / 0). 
Then the system (1.1) on a domain X has the autonomous first integrals 

F 1 : x — > ^( v 2 x) 2 + [v 2 x) 2 ^ exp^ - 2 A2 -5— ^ for all x e X 

and 

~i 1 

V X ' — V X * 

F 9 : x — >■ arctan - — — A 9 -—— for all x £ X, X C {x: u°x ^ 0, z^ 2 x 7^ 0). 

2 l,2 x 2 U°X 1 

Theorem 1.7. Suppose v 01 and u 11 are a real eigenvector and a real 1-th order gen- 
eralized eigenvector of the matrix B corresponding to the eigenvalue X 1 with elementary 
divisor of multiplicity m 1 ^ 2, and v 02 and u 12 are a real eigenvector and a real 1-th order 
generalized eigenvector of the matrix B corresponding to the eigenvalue \ 2 with elementary 
divisor of multiplicity m 2 ^ 2. Then the system (1.1) has the autonomous first integral 
11 12 

F: x -»• V -^- - V -^- for all x G X, X C {x: u 01 x / 0, u 02 x / 0}. 

V UL X U yJZ X 

Indeed, the Lie derivative of the function F on a domain X by virtue of system (1.1) is 

_ (X^x + v m x)v m x - X 1 u m x u u x (X 2 u 12 x + iy 02 x)u 02 x - X 2 v 02 x v 12 x _ 
W ~ (^V) 2 (^p " °- 



Theorem 1.8. Suppose X is the eigenvalue with elementary divisor of multiplicity 
m ^ 2 of the matrix B corresponding to an eigenvector v° and to generalized eigenvec- 
tors v k , k = 1, . . . , m — 1. Then the system (1.1) on a domain X C {x: v°x 0} has the 
functionally independent autonomous first integrals 

F ( : x -> * f (x) forallxeX, ( = 2,...,m-l, (1.16) 

where the functions ^ : X — >■ R, £ = 2, . . . , m — 1, are £/te solution to the system 
k 

v k x = Y J ( k T ~l)^T(x)v k ~ T x forallxeX, k = 1, . . . , m - 1. (1.17) 

T = l 

Proof. Using the identities (1.3) and the equalities (1.14), we obtain 

2lz/ x = A z/°x for all x G ]R n , 

(1.18) 

21 z/x = A z/ fc x + k v k 1 x for all x G R n , fc = 1, . . . , m — 1. 

The system (1.17) has the determinant {v®x)' m ~ l such that (z/°2;) m_1 7^ for all x 
from a domain X C {x: u°x / 0}. Therefore there exists the solution r = 1, . . . , m — 1, 
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on the domain X of the functional system (1.17). Let us show that 

1 for all x G X, k = 1, 



m k (x) 



for all x G X, fc = 2, . . . , m — 1, 



(1.19) 



The proof of identities (1.19) is by induction on m. 

For m = 2 and m = 3, the assertions (1.19) follows from the identities (1.18). 
Assume that the identities (1.19) for m = e is true. Using the system of identities (1.18) 
and the system (1.17) for m = e + 1, m = e, we get 



e-1 



21 i/ £ x = A ^( £ : 1 1 )^ r (x)z/ £ - r a; + (e - l)^(^)* r (x) z^-^x + v e ~ x x + i/°a;atf e (x) 



T=l 



T=l 



for all x G X, X C {x: z/°x / 0}. 

Now taking into account the system (1.17) with k = e — 1 and k = e, the identity (1.18) 
with k = e, and z/°x ^ for all x G X, we have 

m e (x) = for all x G X. 

This implies that the identities (1.19) for m = e + 1 are true. So by the principle of 
mathematical induction, the statement (1.19) is true for every natural number m ^ 2. 

Thus the functions (1.16) are functionally independent first integrals of system (1.1). ■ 

Theorem 1.8 is true both for the case of the real eigenvalue A and for the case of the 
complex eigenvalue A (Im A / 0). In the complex case, from the complex- valued first integrals 
(1.16) of system (1.1), we obtain the real-valued first integrals of system (1.1) 

F£ : x ->■ Re* ( (x), F f 2 : x -)• Im <& f (x) for all x G X, £ = 2, ...,m-l, 

where X is a domain from the set {x: (z/°x) 2 + (z?°x) 2 ^ 0} C R ra . 
Example 1.5. Consider the linear autonomous differential system 



dx 1 
~~dt 



The matrix B 



Ax i x 2 , 



4 
1 




dx 2 
~dt 

3 1 
1 
1 1 



dx 3 
~~dt 



x l + x 3- 



(1.20) 



has one triple eigenvalue A 1 = 2. 



The rank of the matrix B — X 1 E is equal 2. Therefore the eigenvalue A 1 = 2 has 
k 1 = 3 — 2 = 1 elementary divisor (A — 2) 3 of multiplicity 3. 

2u l + 3v 2 + 1/3 = 0, c 
The system (f? — A 1 £J) colon ( f 2 , 1/3) = <^ ^ — u 1 — u 2 = 0, <^ 

— l/n — ZA> = 



'3> 



Hence v° = (1, - 1, 1) is an eigenvector of B corresponding to the eigenvalue A 1 = 2. 

The 1-st order generalized eigenvector u l of the matrix B corresponding to the eigen- 
value Ai = 2 is a solution of the system 

2u l + 3v 2 + ^3 = 1) 

(B — XiE) colon ( u 1 , u 2 , u 3 ) = colon(l, — !,!)<*=>< 



- Vl - u 2 = - 1, 4^ 
- i/ 2 - 1/3 = 1 



1. 



Hence z/ 1 = (1,0, - 1) is a generalized eigenvector of the 1-st order of the matrix B corre- 
sponding to the eigenvalue \ l = 2. 

12 



V.N. Gorbuzov, A.F. Pranevich 



First integrals of ordinary linear differential systems 



The 2-nd order generalized eigenvector v 2 of the matrix B corresponding to the eigen- 
value \ = 2 is a solution of the system 

' 2v x + 3u 2 + u 3 = 2, ( 

(B- \E) colon (z/ l5 f 2 , v 3 ) = 2colon(l, 0, — 1) <^ < 



- v x - v 2 = 0, <S> 



"2) 

z/ 9 + 2. 



Hence z/ 2 = (0,0,2) is a generalized eigenvector of the 2-nd order of the matrix B corre- 
sponding to the eigenvalue X 1 = 2. 

The functionally independent scalar functions (by Theorem 1.5) 

-2 - ^ x ^ j for all (x l5 x 2 ,x 3 ) G X, 

F 2 : (x^x^Xg) -> ( x i ~ ~ 2X3^ - x 2 + x 3 ) for all (x 1; x 2 ,x 3 ) G X (by Theorem 1.8) 

(x 1 — x 2 + x 3 ) 

are a basis of autonomous first integrals on a domain X C {(x 1; x 2 , x 3 ) : x 1 — x 2 + x 3 / 0} of 
the ordinary differential system (1.20). 

Example 1.6. The sixth-order ordinary autonomous linear differential system 

dx dx 2 

= x 1 — 2x 2 + x 3 — 2x 6 , = 3x 2 — x 3 — x 5 + 2x 6 , 

~^f = ~ x 1 + x 3 + 2x 4 + 2x 5 , = ~ x i + x 4 + x 5 + x 6i ( L21 ) 

dxr dxs? 

^ /-y> I V • J > • - 1 ' -y* I -'y v i t ■ 

has the triple complex eigenvalue A 1 = 1 + i with the elementary divisor (A — 1 — i) 3 of 
multiplicity 3 corresponding to the eigenvector v = (1, 1, 0, 0, i, 0) and to the generalized 
eigenvectors u 1 = (0, 1, 0, i, i, 1), z/ 2 = (0, 1, i, 0, i, 0). The functions 

F 1 : x ->■ P(x) exp( - 2</?(x)) for all x G X (by Theorem 1.2), 
F 2 : x -> P(x) exp ^ - 2 ^p^y^ ) for all x G X (by Corollary 1.3), 

F 3 : x -»• t^(x) - a ( x ) + for a ll x G X (by Corollary 1.3), 

F 4 : x -> + /3 2 (x)-a 2 (x) for ^ ^ g ^ Theorem Lg) 

P 2 (x) 

6(x)P(x) - 2a(x)/3(x) , „ _ 1 c , 

F. : x ->• o / \ for all x G X (by Theorem 1.8), 

where the polynomials 

P: x — >■ (x 1 + x 2 ) 2 + x 2 , a: x — >■ (x x + x 2 )(x 2 + x 6 ) + x 5 (x 4 + x 5 ), 
0: x-> (x 1 +x 2 )(x 4 + x 5 )-x 5 (x 2 + x 6 ), 7: x^x 2 (x 1 +x 2 )+x 5 (x 3 + x 5 ), 
5 : x — > (x 1 + x 2 )(x 3 + x 5 ) — x 2 x 5 for all x G M 6 , 
the scalar function ip : x — > arctan — for all x G X, are first integrals on a domain X 

X ^ ~\~ Xt^ 

from the set {x : x 1 + x 2 / 0} C M 6 of the ordinary differential system (1.21). 

The functionally independent first integrals F 1 , . . . , P 5 are an autonomous integral basis 
on a domain X of the ordinary differential system (1.21). 
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1.1.3. Nonautonomous first integrals 

The ordinary differential system (1.1) is induced the nonautonomous linear differential 
operator of first order 93 (t, x) = d t + 2l(x) for all (t, x) G R n+1 . 

Adding one nonautonomous first integral of system (1.1) to an autonomous integral basis 
of system (1.1), we can construct a basis of first integrals for system (1.1). 

Such procedure every time can be carried out on the base of the following statements. 

Theorem 1.9. Suppose v is a real eigenvector of the matrix B corresponding to the 
eigenvalue X. Then the ordinary differential system (1.1) has the first integral 

F: (t,x) ->• vx exp( - Xt) for all (t,x) G R n+l . 

Indeed, the Lie derivative of the function F by virtue of system (1.1) is 

*BF(t, x) = d t F(t, x) + OF(t, x)= - XF(t, x) + XF(t, x) = for all (t, x) G R n+1 . ■ 

Example 1.7 (continuation of Example 1.1). By Theorem 1.9, using the eigenvalue 
A 2 = 1 corresponding to the eigenvector v 2 = (2,2,1,1), we can build the first integral 
F 2 : (t, x) ->• (2xi + 2x 2 + x 3 + x 4 ) e~* for all (t, x) G R 5 of system (1.9). 

The functions F 1 , F 23 , F 24 (constructed in Example 1.1), and the function F 2 are an 
integral basis of the autonomous differential system (1.9) on a domain Kxl, where X is a 
domain from the set {x: x 1 + x 3 / 0, 2x 2 + x 4 / 0} C R 4 . 

Example 1.8 (continuation of Example 1.2). The autonomous differential system (1.12) 
has the eigenvector v 2 = (3, — 1, — 1) corresponding to the eigenvalue A 2 = 2. 
By Theorem 1.9, the differential system (1.12) has the first integral 

F 3 : (t, x 1; x 2 , x 3 ) — > (3x 1 — x 2 — x 3 ) e~ 2t for all (t, x 1 , x 2 , x 3 ) G R 4 . 

The functionally independent first integrals F 1: F 2 (constructed in Example 1.1), and 
F 3 are a basis of first integrals for the differential system (1.12) on a domain Rx 1, where 
X is a domain from the set {(x 1 ,x 2 ,x 3 ): x 1 — x 3 / 0}. 

Corollary 1.5. Let v = v + vi (Rev = v, \m.v = v) be an eigenvector of the matrix B 

corresponding to the complex eigenvalue X = X + Xi (Re A = A, ImA = A / 0). Then the 
ordinary differential system (1.1) has the first integrals 

F 1 : (t, x) -»• ((^x) 2 + (Vx) 2 ) exp( - 2A t) for all (t, x) G R n+1 

and 

jyX — ' 

F 2 : (t, x) — > arctan — — Xt for all (i, x) G R x X, 
vx 

where X is a domain from the set {x: fx / 0} C R n . 
Proof. Using Properties 1.4 and 1.5, we get 

fSF^x) = exp(-2At)5S((^x) 2 + (?x) 2 ) + ((^x) 2 + {Vx) 2 ) 5Sexp( - 2\t) = 
= 2A((^x) 2 + (Vx) 2 ) exp(-2At) + ((^x) 2 + (Vx) 2 )d t exp( - 2\t) = for all (t,x)eR n+1 , 

5SF 2 (t,x) = 55 arctan- 93 (At) = A - d t (Xt) = for all (t,x) G R x X. 

vx 

Therefore the scalar functions F 1 : IR ra+1 — > R and F 2 : R x X — > R are first integrals of 
the ordinary linear autonomous differential system (1.1). ■ 
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Example 1.9 (continuation of Example 1.3). Using the eigenvalue A 1 = i corresponding 
to the eigenvector v 1 = (1 — i, — 1 + 2i, 2i, 2), we can construct the first integral 

^ I 2^7 I 2^7 

F 4 : (t,x) -> arctan l - ^— ^ - t for all (t,x) Glxlj (by Corollary 1.5) 

of the differential system (1.13) on a domain X 1 from the set {x: x 1 — x 2 + 2x 4 7^ 0}. 

The functionally independent first integrals F x , F 2 , F 3 (constructed in Example 1.3), and 
the function F 4 are an integral basis of the differential system (1.13) on a domain IxX, 
where X is a domain from the set {x: x 1 — x 2 + 2x 4 / 0, x 4 — x 2 7^ 0}. 

Theorem 1.10. Lei A 6e i/ie eigenvalue with elementary divisor of multiplicity m ^ 2 
0/ t/ie matrix B corresponding to a real eigenvector v° and to a real generalized eigenvector 
v 1 of the 1-st order. Then the system (1.1) has the first integral 

F: (t,x) -»• for all (t,x) eKxl, X C {x : z/°x / 0}. (1.22) 

Indeed, the Lie derivative of F on a domain R x 1 by virtue of system (1.1) is 

„ . „ . . _ (Xv 1 x + z/°x)^°x — z/ 1 x Ai/°x „ _ 

<8F (t, x) = d t F(t, x) + 2tF(t, x) = - 1 + ^ -Vto = 0. ■ 

(z^ u x) z 

Example 1.10 (continuation of Example 1.4). The system (1.15) has the eigenvalue 
X 1 = corresponding to the eigenvector v 1 = (1, — 2, 1) and to the 1-st order generalized 
eigenvector z^ 2 = (0, — 1, 1). By Theorem 1.10, the scalar function 

F 3 : (t,x 1 ,x 2 ,x 3 ) ->• X3 ^ t for all (t, x 1; x 2 , x 3 ) G M 4 

X-j^ ZiX<2 \ X^ 

is a first integral of the linear differential system (1.15). 

The functionally independent first integrals F 1 , F 2 (constructed in Example 1.4), and 
F 3 are an integral basis of the differential system (1.15) on a domain RxX, where X is a 
domain from the set {(x 1 ,x 2 ,x 3 ) : x 1 — 2x 2 + x 3 / 0}. 

Example 1.11 (continuation of Example 1.5). The system (1.20) has the eigenvalue 
X 1 = 2 corresponding to the eigenvector v° = (1, — 1, 1) and to the 1-st order generalized 
eigenvector v 1 = (1,0, — 1). Using Theorem 1.10, we can build the first integral 

F 3 : (t,x 1 ,x 2 ,x 3 ) — > — — — — — t for all (t, x l ,x 2 , x 3 ) G R 4 . 

of the linear autonomous differential system (1.20) 

The functionally independent first integrals F ± , F 2 (constructed in Example 1.5), and 
F 3 are an integral basis of the differential system (1.20) on a domain KxX, where X is a 
domain from the set {(x 1 , x 2 , x 3 ) : x 1 — x 2 + x 3 7^ 0}. 

If A is a complex number, then the function (1.22) is a complex-valued first integral of 
system (1.1). In this case, from the complex- valued first integral (1.22) of system (1.1), we 
obtain two real-valued first integrals of system (1.1) 

V X f x + V X V X 



F 1 : (t, x) -> — t for all (t, x) G R X X 

(v°x) + (v°x) 



and 



~1 ~0 1 

V u x V x — V XV x 



F 9 : (t, x) -»• — — for all (t, x) G K x X, 

(>x) 2 +(?o x ) 2 

where X is a domain from the set {x: (y °x) 2 + (z?°x) 2 7^ 0}. 
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F 6 :(t,s)-> ^1^^2^-r^4-r^ _ ^ for ^ ( tjX ) €RxX . 



Example 1.12 (continuation of Example 1.6). Using the eigenvector v° = (1, 1,0,0,2,0) 
and 1-st order generalized eigenvector v l = (0, 1, 0, i, i, 1) corresponding to the eigenvalue 
X 1 = 1 + i, we can construct the first integral of system (1.21) 

{x 1 + x 2 ){x 2 + x 6 ) + x 5 (x 4 + x 5 ) 

The functionally independent first integrals F X ,...,F 5 (constructed in Example 1.6), and 
F 6 are a basis of first integrals for the differential system (1.21) on a domain Rx 1, where 
X is a domain from the set {x: x 1 + x 2 7= 0}. 

Let us remark that using Theorems 1.8, 1.9, 1.10, and Corollary 1.5, we can always build 
an integral basis of system (1.1) on the base of nonautonomous first integrals. 

Example 1.13. The fifth-order ordinary linear autonomous differential system 



dx-^ . . dx<2 . . o o dx^ 

—jj- = 2{x 1 — Sx 2 + 2x 3 + £5), —j^- = 2x 1 — x 2 + 2x 3 + 2x 4 , = x 2 — x 3 , 



(1.23) 



^*^4 . dx^ 

— oX-^ "I - 3 5' ~dt~ — 2 3 ~~ ^" 



has the eigenvalues A 1 = — 1, A 2 = 1 — i, and A 3 = 1 + i with elementary divisors (A + l) 3 , 
A — 1 + i, and A — 1 — i, respectively, corresponding to the eigenvectors u 01 = (1,0,1, 1,0), 
v 2 = (1, 1 + i, 0, 1 + i, 1 — i), v 3 = (1, 1 — i, 0, 1 — i, 1 + i) and to the generalized eigenvectors 
u n = (1,1/2,1,1,0), u 21 = (0,1,1,0,0). 

Using the eigenvector u 01 = (1,0, 1,1,0), the generalized eigenvector u n = (1,1/2, 1,1,0) 
of the 1-st order, and the generalized eigenvector u 21 = (0, 1, 1, 0, 0) of the 2-nd order corre- 
sponding to the eigenvalue A 1 = — 1, we can build the first integrals of system (1.23): 

F 1 : (t, x) -»• (x 1 + x 3 + x 4 ) e* for all (t, x) G M 6 (by Theorem 1.9), 

2r -4- t A- 2r -I- 2t 

F 2 : (t,x) -> V " * for a11 S M x Xi (by Theorem 1.10), 

2 V ; 2(x x + x 3 + x 4 ) v 7 7 ' 

and (by Theorem 1.8) 

F 3 : (t,x) -> 4( " 2 + ^ + * 3 + ^ " (2Xl + * 2 + 2X3 + 2X4)2 for all (t,x) G R x X l5 

where X 1 is a domain from the set {x: x 1 + x 3 + x 4 7^ 0} Cl 5 . 

Using the eigenvector zv 2 = (1, l + i, 0, 1 + i, 1 — corresponding to the complex eigenvalue 
A 2 = 1 — i, we can construct (by Corollary 1.5) the first integrals of system (1.23): 

F 4 : (t, x) ->• ((x x + x 2 + x 4 + x 5 ) 2 + (x 2 + x 4 - x 5 ) 2 ) e -2 * for all (t, x) G IR 6 

and 

F 5 : (t, x) -»• arctan — ^ ^1^5 — + f for aU ^ a;) G R x X 2 , 

X-^ "t - X 2 "T x 4 -|- x^ 

where X 2 is a domain from the set {x : x 1 + x 2 + x 4 + x 5 7= 0} CK 5 . 

The functionally independent first integrals F 1 ,...,F 5 are an integral basis of the linear 
autonomous differential system (1.23) on a domain Rxl, where X is a domain from the 
set {x : x 1 + x 3 + x 4 / 0, x l + x 2 + x 4 + x 5 / 0} C R 5 . 
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1.2. Linear nonhomogeneous differential system 

Consider an ordinary linear nonhomogeneous differential system with constant coefficients 

- = Ax + f(t), (1.24) 

where /: t — >■ colon(/ 1 (t), . . . , f n (t)) for all t G J is a continuous function on an interval 
Jcl. The differential system (1.1) is the corresponding homogeneous system of (1.24). 
The system (1.24) is induced the linear differential operator of first order 

n 

t(t,x) = dt + ^2{A ( x + f ( (t))d Xe for all [t,x)eJxR n , A^ = (a^, . . . ,a^ n ), £ = 1, . . . ,n. 
5=1 

1.2.1. Case of simple elementary divisors 

If the matrix B has simple structure, then we can build first integrals of the differential 
system (1.24) by using following assertions (Theorem 1.11 and Corollary 1.6). 

Theorem 1.11. Let u be a real eigenvector of the matrix B corresponding to the eigen- 
value A. Then a first integral of system (1.24) is the scalar function 

t 

F: (t,x) -»• vx exp(- At) - J i//(C) exp( - \()d( for all (t,x) G J x R n , (1.25) 

*o 

where t is a fixed point from the interval J. 
Proof. Using Lemma 1.1, we get 

n 

€F(t,x) = d t F(t,x)+$lF(t,x) + ^2fz(t)d x F(t,x) = -Xvxexp(-Xt)- i//(i)exp(-At) + 

5=i £ 

n 

+ A vx exp( -\t) + ^2 "{/{W ex P( - At) = for all (t, x) G J x R n . 
5=1 

Therefore the function (1.25) is a first integral of the differential system (1.24). ■ 
Corollary 1.6. Suppose v = v + vi (Rev = v, \mv = v) is an eigenvector of the matrix 

B corresponding to the complex eigenvalue A = A + Ai (Re A = A, ImA = A 7^ 0). Then first 
integrals of the differential system (1.24) are the scalar functions 

t 

F e : (t,x) ^ a e {t,x) - J a d {(,f(())d( for all (t,x) G J xR", 9 = 1,2, (1.26) 

*o 

where t is a /ixed point /rom t/i.e interval J, the functions 

a 1 : (t, x) — > ( vx cos At + vx sin At) exp( — A t) for all (t, x) G IR ra+1 , 

a 2 : (t,x) — >■ (zat cosAt — ^x sinAt) exp( — At) for all (t, x) G IR n+1 . 

Proof. From Theorem 1.11 it follows that the function (1.25) is a complex- valued first 
integral of the differential system (1.24). Then, the real and imaginary parts of this function 
are the real- valued first integrals (1-26) of the differential system (1.24). ■ 

Example 1.13. Let us consider the linear nonhomogeneous differential system 

^P- = 2aN + x 2 + 2e 2t , ^ = x, + 3x 2 - x, + 10, ^ = - x, + 2x 2 + 3x, + e 3 '. (1.27) 
dt 12 dt 1 2 6 dt i2i \ j 

The differential system (1.12) is the corresponding homogeneous system of (1.27). 
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Using the eigenvectors v 1 = —1), v 2 = (3, —1, —1) corresponding to the eigenvalues 
A x = 3 + i, A 2 = 2, respectively, we can build (by Corollary 1.6 and Theorem 1.11) the 
functionally independent first integrals (an integral basis) of system (1.27) 

F l : (t, x l ,x 2 ,x 3 ) — > ((x 1 — x 3 + 1) cos t + (x 2 + 3) sint)e~ 3 ' + (cost — sint)e~* + sint, 
F 2 : (t, x 1 ,x 2 ,x 3 ) — > ((x 2 + 3) cost + (x 3 — x l — l) sint)e~ 3 * — (cos t + sint)e~* + cost, 
F 3 : (t, x 1 , x 2 , x 3 ) — > (3x 1 — x 2 — x 3 — 5)e~ 2t + e* — 6t for all (t, x l5 x 2 , x 3 ) G M 4 . 

Remark. Under thr conditions of Corollary 1.6, we have the scalar function 

t 

F: (t,x)-> ((^) 2 + (?x) 2 )exp(-2At) - 2^(4, a) | a 1 (C,/(C))dC + 



t 



+ a 2 (t,x) y"a 2 (C,/(C))dc) + (/«i(C/(C))dc) + (/a 2 (C,/(C))dc) 

*0 *0 *0 

is also a first integral on the domain JxM™ of the linear differential system (1.24). 
1.2.2. Case of multiple elementary divisors 

If the matrix B has multiple elementary divisors, then we can construct first integrals of 
system (1.24) by using following assertions (Theorem 1.12 and Corollary 1.7). 

Theorem 1.12. Let A be the eigenvalue with elementary divisor of multiplicity m ^ 2 
of the matrix B corresponding to a real eigenvector u° and to real generalized eigenvectors 
u h , k = 1, . . . ,m — 1. Then the system (1.24) has the functionally independent first integrals 

k-l 

F k+1 : (t,x) -> ^xexp(-At) - £g) t k ' T F T+1 (t, x) - C k {t) 

T = 

(1.28) 

for all (t,x) G J x E n , fc = 1, . . . , m - 1, 

where the integral F 1 : (t,x)— >■ z^°xexp(— At) — C (t) /or a// (t,x)G Jxl" (fry Theoreml.ll), 
the scalar functions 
t 

C k :t^ J (is k f(()exp( - AC) + fcC fc _i(C)) d C for all t £ J, k = 0, . . . , m - 1, t G J. 

*o 

Proof. The proof of Theorem 1.12 is by induction on m. 
By the equalities (1.14), it follows that 



£(z/ e x exp(-At)) = {ev^x + ^fit)) exp( - At) 

for all (t, x) G J x R n , e = 1 . . . , m - 1. 
Let m = 2. Using the system of identities (1.29), we get 

£F 2 (t,x) = <t{y x x exp( - At) - tF x (t,x) - C x (tj) = 
= (z.°x + ^ 1 /(t))exp(-At)-F 1 (t,x)- (^/(t) exp( - At) + C (t)) 
= (z/>x exp( - At) - C (t)) - F x (t, x) = for all (t, x) G J x M n . 

Therefore the function F 2 : J x R n — > R is a first integral of system (1.24). 
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Assume that the functions (1-28) for m = fi are first integrals of system (1.24). Then, 
from the identity (1.29) for the function F +1 it follows that 

fi-i 

(t, x) = £ exp(- At) - ^ (^^^ 



T=0 

11-1 



-V ECr 1 )*" T F T+1 (t,x)- t iF lx (t,x)-{^f(t)eM-M) + »C l ,_ 1 (t)) 



r=0 

/i-2 



= ^- 1 x-Y J { fl r 1 ) t ^ " 1 i ? T+ i(^)-C (1 _ 1 (t))-^(f,x)=0 for all (i,i)€JxR n . 

T = 

This implies that the scalar function F +1 : J x R" R (for m = is a first integral 

of the linear nonhomogeneous differential system (1.24). 

So by the principle of mathematical induction, the scalar functions (1.28) are first integrals 
of the differential system (1.24) for every natural number m ^ 2. ■ 

Example 1.14. Consider the linear nonhomogeneous differential system 

-± = 4x 1 - x 2 + e 3t , = 3x 1 +x 2 -x 3 + 8t, = x x + x 3 + 4. (1.30) 

The differential system (1.30) has the eigenvalue Ai = 2 corresponding to the eigenvector 
u° = (1, — 1,1) and to the generalized eigenvectors v l = (1,0, — 1), v 2 = ( — 2,2,0). 
An integral basis of system (1.30) is the functions (by Theorems 1.11 and 1.12) 

F 1 : (t, x) -»• (xi - x 2 + x 3 - 4t)e" 2 * - e* for all (t, x) G M 4 , 
F 2 : (t, x) -»• (xi - x 3 + 2t - l)e" 2 * - t F x (t, x) - 2e* for all (t, x) G M 4 , 
F 3 :(t,x)-> 2(x 2 -xi+3t + 2)e~ 2 *-t 2 F 1 (t,x)-2tF 2 (t,x)-2e' for all (t,x)GlR 4 . 

Corollary 1.7. Lei A = A + Ai (Re A = A, ImA = A ^ 0) 6e t/je complex eigenvalue of 
the matrix B with elementary divisor of multiplicity m ^ 2 corresponding to an eigenvector 
i>° = v° + T/°i (Re^° = v°, Imu = u°) and to generalized eigenvectors v k = v k + u k i 
(Rei ,k = v k , lmu k = v k ), k = 1, . . . , m—1. Then first integrals of the linear nonhomogeneous 
differential system (1.24) are the functions 

k-l 

F Q k+1 : (t, x) a Qk (*,*)-£ ffl ^ *i >T+ i (*, *) - C Qk (t) 

T=0 

(1.31) 

for all (t,x) G J xR", fc = 1, . . . , m - 1, (9 = 1, 2, 
where the first integrals (by Corollary 1.6) 

(t,x) -> a fl0 (t,x) - C gQ (t) for all (t, x) G J x R n , = 1,2, 
t/ie scalar functions 

a lk : (t,x) (v k x cos\t + v k x sin At) exp( - \ t) for all (t, x) G R n+1 , k = 0, . . . , m - 1, 

a 2fc : (t,x) ->■ (z/ fc x cos At - v k x sin At) exp( - At) for all (t, x) G IR n+1 , k = 0, . . . , m — 1, 
t 

C efc :t^ J{a ek ((J(C)) + kC ek l (C))d( forall t£j, fc = 0,...,m-l, = 1,2, t G J. 
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Proof. Formally using Theorem 1.12, we get the complex-valued first integrals (1.28) of 
the differential system (1.24). Then, the real and imaginary parts of the functions (1.28) are 
the real- valued first integrals (1.31) of the differential system (1.24). ■ 

Example 1.15. The linear nonhomogeneous differential system 



dxi dx 2 
—— = — Xq + x A + 4a cos t, —r~ =ii + x q + 4a sin t, 
dt ^4 dt i d 



dx 3 dx 4 c 

—j— = — x 2 + x 3 — x 4 + 6t, —r- = — x 1 + x 2 + x 3 — x 4 H . 2 , > 
at " at sin i 



(1.32) 



where a, 6, and c are some real numbers, has the complex eigenvalues A 1 = i, A 2 = — i 
with elementary divisors (A — i) 2 , (A + i) 2 , respectively. 

The number A x = i is the eigenvalue corresponding to the eigenvector v° = (1,2,1,0) 
and to the generalized eigenvector u 1 = ( — 1 + i, 0, 0, i) of the 1-st order. 

Using the real numbers Ai =0, A 1 = 1, the real vectors is = (1, 0, 1, 0), v° = (0, 1, 0, 0), 
v x = ( — 1,0,0,0), v x = (1,0,0,1), and the scalar functions 

a 1Q : (t, x) — > cos t (x 1 + x 3 ) + sint x 2 , « 20 : (t, x) — > cos t x 2 — sint {x l + x 3 ), 

a u : (t,x) — > — cos t x 1 + sint (x 1 + x 4 ), a 2± : (t, x) — > cost (x-y + x 4 ) + sint x 1 

for all (t, i)eJ,x M 4 , J, = (ttZ; tt(Z + 1)) for all / G Z, 

C 1Q (t) = J (4a+bt cos t)dt = 4at+6(cost+t sint), C 2Q (t) = — J btsintdt = 6(t cost— sint), 

C u (t) = J ( - 4a cos 2 t + 4a cos t sin t + + 4at + 6(cos t + t sin t)) dt = 
= 2a(t 2 — t + sin 2 t — sin t cost) + 6(2 sin t — t cost) + c In tan ^ for all i£ J ( , 

if c COS ^ \ 

C 21 (t) = / ^4a cos 2 t -I — r^-j- + 4a cos t sin t + b(t cos t — sin t) J dt = 

= 2a (t + sin 2 t + sin t cos t) + 6(2 cos t + t sin t) — for all t G J, , 

smt 4 

we can build (by Corollaries 1.6 and 1.7) the first integrals of system (1.32) 

F n : (t, x) — > cos t (x 1 + x 3 ) + sin tx 2 — Aat — 6(cos t + t sin t) for all (t, i) £ J ( x M 4 , 

F 21 : (t, x) — > cos tx 2 — sin t (x x + x 3 ) + 6(sin t — t cos t) for all (t, x) G J z x M 4 , 
F 12 : (t, x) — >■ — cos tx 1 + sin t (x x + x 4 ) — tF u (t, x) — 2a(t 2 — t + sin 2 t — sin t cos t) — 

— 6(2 sin t - t cos t) — cln tan - for all (t, x) G J l x R 4 , 

F 22 : (t, x) — >■ cost (x x + x 4 ) + sintxj^ — tF 21 (t,x) — 2a(t + sin 2 t + sin t cost) — 

- 6(2 cos t + t sin t) + for all (t, x) G J, x M 4 . 

smt 

The functionally independent first integrals F u , F 21 , F 12 , and F 22 are a basis of first 
integrals for the differential system (1.32) on a domain J t x R 4 . 
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2. Integrals of ordinary linear nonautonomous differential systems 

integrable in closed form 

2.1. Algebraic reducible systems 

Consider an ordinary real linear nonhomogeneous differential system of the n-th order 

UT 

- = A(t)x + f(t), (2.1) 

where x = colon(x 1 , . . . ,x n ) G W 1 , the continuous on an interval J C M coefficient matrix 
A: t — > A(t) for all £ G J is diagonalizable by a constant semilarity matrix [59, p. 61], the 
vector function /: £ — > colon(/ 1 (t), . . . , f n (t)) for all t G J is continuous. 

The corresponding homogeneous system of the nonhomogeneous system (2.1) is 

d i=A {t)x . (2.2) 

By [48, p. 186], the nonhomogeneous differential system (2.1) and the corresponding 
homogeneous differential system (2.2) are called algebraic reducible. 

2.1.1. Partial integrals 

A complex-valued linear homogeneous function 

p:x^ux forallxGM n , v G C n , (2.3) 

is a partial integral of the algebraic reducible system (2.2) if and only if the identity holds 

2lp(x) = \(t)p(x) for all (t,x) G J x W\ (2.4) 

where the linear differential operator 2t(t, x) = d t + A(t)xd x for all (t, x) G J x IR n , the 
scalar function A: J — > R. The identity (2.4) is equivalent to the linear system 

(B(t) - X(t)E)u = 0, (2.5) 

where E is the nx n identity matrix, and the matrix B: t — > B(t) for all t G J is the 
transpose of the matrix A: t — > A(t) for all £ G J. 

The following basic propositions (Lemmas 2.1 and 2.2) are base for the method of building 
first integrals of the algebraic reducible systems (2.1) and (2.2). 

Lemma 2.1. Suppose v is a real eigenvector of the matrix B : t — > B(t) for all t G J 
corresponding to the eigenfunction A: t — > A(i) for all t G J. Then the linear function (2.3) 
is a partial integral of the algebraic reducible differential system (2.2). 

Indeed, if v G W l is an eigenvector of the matrix B : t — » B(t) for all £ € J corresponding 
to the eigenfunction A: t — > X(t) for all t G J, then v is a solution to the system (2.5). This 
implies that the identity (2.4) holds. Thus the linear function (2.3) is a partial integral of the 
algebraic reducible differential system (2.2). ■ 

Lemma 2.2. Let v = v + v i [y = Rez/, v = Im^ ^ 0) be a complex eigenvector of the 

matrix B: t — > B(t) for all t G J corresponding to the eigenfunction A: t — > \(t) + A(t) i 

/or all t £ J (\: t ReA(£), A: t — > ImA(£) for all t G J), let the functions P and ip 

be defined by P: x — > {vx) 2 + (yx) 2 for all x G R n anc? tp: x ^ arctan — for all x G X, 
where X is a domain from the set {x: / 0} C R n . Then, we have 

2lP(x) = 2A(£)P(x) /or a// (£,x) G J x M n and 2t^(x) = A(£) for all (£,x) G J x X. 

Proof. Formally using Lemma 2.1, we get the complex- valued function (2.3) is a partial 
integral of the differential system (2.2) and the following identity holds 
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VL{vx + i vx) = (\(t) + i X(t)) {tx + i vx) for all (t, x) G J x R n . 
This complex identity is equivalent to the real system of identities 

%vx = X(t)vx-X(t)vx, %vx = X(t)vx + X(t)vx for all (t, x) G J x R n . 
Using this system of identities, we obtain 

2lP(x) = 2t((zV) 2 + (zat) 2 ) = 2^x21^ + 2vx%vx = 2vx (x(t) vx - X(t)vx) + 

+ 2vx (X(t) vx + ~X(t) tx) = 2X(t) {(vx? + {vx?) = 2 A(f) P{x) for all (t, x) G J x R« 
and 

= ^(AW^ + A(t)^-^(A(t)^-A(t)^) = for ^ M £jxxa 

(vx) 2 + (zat) 2 

2.1.2. First integrals 

Using Theorems 2.1 and 2.2, we can obtain first integrals of the algebraic reducible non- 
homogeneous differential system (2.1). 

Theorem 2.1. Suppose v is a real eigenvector of the matrix B : t — > B(t) for all t G J 
corresponding to the eigenfunction A: t — > X(t) for all t £ J. Then the algebraic reducible 
differential system (2.1) has the first integral 

t 

F: (t,x) -»• vx<p(t) - j vf(()<p(()d( for all (t,x) G J x R n , (2.6) 
*o 

where the exponential function 

t 

<p: t ->■ exp^- y A(C) dcj / or «^ i G J, 
*o 

and i is a fixed point from the interval J. 

Proof. From Lemma 2.1, we get the following 

t 

<BF(t, x) = 2LF(i, x) + f(t) 8 x F(t, x) = vx 8 t <p(t) -d t j vf(C) <p(() d( + <p(t) 21 vx + 

*o 

+ f{t)d x vxip(t) = - X(t)vxip(t)-vf(t)ip(t) + X{t)vxip{t) + vf(t)ip(t) = 

for all (t,x) G J x W 1 , 

where the linear differential operator 93 (i, x) = d t + (-A(t) x + f(t)) d x for all (t, x) G J x R n 
is the operator of differentiation by virtue of system (2.1). 

Therefore the function (2.6) is a first integral of the algebraic reducible system (2.1). ■ 

Example 2.1. Let us consider the linear differential system of the third-order 

= tx 3 (t)x 1 + (a^t) - a 3 (t))(x 3 -x 2 ) + f 1 (t), 

^ = (a 2 (t) - a 3 (t)) (x 3 - Xl ) + a 3 (t) x 2 + f 2 (t), (2.7) 

^ = (a 3 (*) " " 2 (*)) *i + " «i(*)) ^2 + («i(*) + a 2 (*) " <*s(*)) *3 + / 3 (*)> 
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where the scalar functions a^: J — > R and : J — >■ R, £ = 1, 2, 3, are continuous on an 
interval J C R. The coefficient matrix 



A: t 



"3^) «3( i ) - Q l( £ ) -a 3 (t) 

a 3 (t)-a 2 (t) a 3 (i) -a^t) a x (*) + a 2 (t) - a 3 (t) 



for all t G J 



such that the transposed matrix B: t A T (t) for all t G J has the real eigenvectors 
v l = (0, — 1, 1), v 2 = ( — 1, 0, 1), and z^ 3 = (1, 1, — 1) corresponding to the eigenfunctions 
A^: t — > ag(t) for all t G J, £ = 1,2,3, respectively. 

Therefore the differential system (2.7) is algebraic reducible. 

By Theorem 2.1, we can build the first integrals on the domain J x R 3 of system (2.7) 

t t C 

Fl :(t,x)^ (x 3 -x 2 )e W (^-Ja 1 (OdC S j - /(/ 3 (C) " / 2 (0) «q>( - J a x {9)dB^, 

*o *o Co 

t t C 

F 2 : (t,x)-> (x 3 -x 1 )exp^-Ja 2 (OdCj ~ J {MO ~ fi(0) exp( - | a 2 (9) dfi) d(, 

*0 *0 Co 

t t C 

F 3 : (t,s)-> (x 1 +x 2 -x 3 )ex P ^y > « 3 (C)dc)-/(/i(C)+/2(O-/3(O)exp^-y'a 3 (0)^ dC, 

*0 *0 C() 

where t and £ are fixed points from the interval J. 

The functionally independent first integrals F l5 F 2 , and F 3 are an integral basis on the 
domain J x R 3 of the algebraic reducible differential system (2.7). 

Theorem 2.2. Let v = v + vi (y = Re ^, v = Imu ^ 0) 6e a complex eigenvector of the 
matrix B: t —> B{t) for all t G J corresponding to the eigenfunction A: t — >■ A(i) + A(t) i 

/or a// £ G J (A: t — > ReA(i), A: t — > ImA(t) for all t G J). Then the algebraic reducible 
differential system (2.1) has the first integrals 

t 

F T : (t, x) -> 7 T (t, s) - y 7t (C, /(C)) M a// (t, x) G J x R n , r = 1, r = 2, (2.8) 
*o 

where t is a fixed point from the interval J, the scalar functions 
t t t 

7l : (t ,x) - (&«./X(0« + «-./a(0«)^(-/X({)-c) «,«> EJ xr, 

t Q t Q <Q 

t t t 

72 : <«„) ( feM ./A(0«-fc-n/X(0*)«p(-/X(0-c) /«-« <«.«) ^ ^xr. 



Proof. Formally using Theorem 2.1, we have the complex-valued function (2.6) is a first 
integral of the algebraic reducible system (2.1). Then the real and imaginary parts of this 
complex- valued first integral are the real- valued first integrals (2.8) of system (2.1). ■ 
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Example 2.2. The algebraic reducible differential system 

= — tanh tx 1 — cosh t (x 2 — x 3 ) — 3, 



^ = - 1 - t&nhtjx 1 + (l-j- coshi)x 2 + cosh£x 3 + 2t 3 + sint - e*, (2.9) 
= ^- — 1 — tanh t — cosh t^j x 1 + ^1 — — + tanh t — cosh t^j x 2 + (cosh t — tanh t) x 3 + 



(it 

+ 2t 3 + sin t - e* - sinh 2 1 



has the eigenvectors v 1 = ( — 1, 1,0), z> 2 = — i), and v 3 = (1, — i,i) corresponding to 
the eigenfunctions \ l : t — > 1 — — , A 2 : t — > — tanh i + cosh t i, and A 3 : t — > — tanh t — cosh i i 
for all i G J, respectively, where J is an interval from the set {t: t ^ 0}. 

Using the eigenvector v 1 = ( — 1,1,0) and the corresponding eigenfunction A x : t — > 1 — j 
for all i G J, we can build (by Theorem 2.1) the first integral 

t 2 / 11 \ 

F x : (t,x) ->• te _i (x 2 -x 1 ) + — + (^2t 4 + 8t 3 + 24t 2 + 45t + 45+- cost + - t(cost + sint)J 



on the domain Jxl 3 of the algebraic reducible differential system (2.9). 

Using the eigenvector v 2 = — i) and the eigenfunction A 2 : ^ — > — tanht + coshti 

for all t G J, we can build (by Theorem 2.2) the first integrals of system (2.9) 

F 2 : (t,x) — > cosh t (cos (sinh t) x 1 + sin(sinht) (x 2 — x 3 )) + 

+ (sinh 2 t - 2) cos(sinht) + (3 - 2 sinh t) sin(sinht) for all J x M 3 , 

F 3 : (t,x) — > coshi( — sin(sinht) x 1 + cos(sinhi) (x 2 — x 3 )) + 

+ (3 - 2sinht)cos(sinht) + (2 - sinh 2 1) sin(sinht) for all J x E 3 . 

The functionally independent first integrals F v F 2 , and F 3 are an integral basis on the 
domain JxR 3 of the algebraic reducible differential system (2.9). 

In the case of the homogeneous algebraic reducible differential system (2.2), we have the 
following statements (Corollary 2.1 and Theorem 2.3). 

Corollary 2.1. Under the conditions of Theorem 2.1, we have the scalar function 

t 

F : (t, x) -»• vx exp ^ - J A(C) d(J for all (t, x) £ J x W\ t G J, 

*o 

is a first integral of the homogeneous algebraic reducible differential system (2.2). 

Theorem 2.3. Let the assumptions of Theorem 2.2 hold, then the scalar functions 

t 

F 1 : (t, x) -> ( (vx) 2 + (?x) 2 ) exp ^ - 2 ^ A(C) ^ for all (t , x) £ J x R n 

*o 

and 

t 

z^x /" — - ^ 
F 2 : (i, x) — > arctan — — / X(Q d( for all (t, x) G J x X, t G J, lc{i:ra^0}, 

ZAT J 



*0 



are first integrals of the homogeneous algebraic reducible differential system (2.2). 
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Proof. From Lemma 2.2, we get 

t t 
f^F 1 (t,x) = exp(-2 J A(C)dc) 2l((£x) 2 + (£x) 2 ) + ({Ixf + (£x) 2 ) 21 exp 2 ^ A(C) 



i t 







= 2A(t) (ipx) 2 + (i/x) 2 ) exp (-2 J \{Q d^j - 2\{t) ({yx) 2 + (ux) 2 ) exp 2 ^ A(C) <%j = 

for all (t,x) eJxr, 
t 

Z^X /" — - ~ ~ 

2lF 2 (t,x) = 21 arctan- 21 / A(C)d( = X(t) - X(t) = for all (t,x) 6 J x X. ■ 

ZA7! J 



Z^X 

*o 

Example 2.3. Consider the algebraic reducible differential system 
dx dx 

—± = a 2 (t) x x + a 3 (t) (x 3 - x 2 ), —^ = (a 2 (t) - a^t)) x x + (a^t) - a 3 (t)) x 2 + a 3 (t) x 3 , 

(2.10) 

_3 = (Q, 2 (t) _ ai (t) _ a3 (t)) Xl + (a^t) - a 2 (t) - a 3 (t)) x 2 + (a 2 (t) + a 3 (t)) x 3 , 

where the scalar functions : J — >■ R, £=1,2,3, are continuous on an interval JcM. 

The system (2.10) has the eigenvectors v 1 = (—1,1,0), v 2 = —i), and v 3 = 
corresponding to the eigenfunctions A 1 : t — > a 1 (t) for all t G J, A 2 : t — > a 2 (t) + ct 3 (t) i for 
all i G J, and A 3 : t — > a 2 (t) — a 3 (t) i for all t £ J, respectively. 

Using the eigenvector v 1 = ( — 1,1, 0) and the corresponding eigenfunction A 1 : t — > a 1 (t) 
for all t G J, we can build (by Corollary 2.1) the first integral 

t 

F 1 : (t,x) ->• (x 2 -x^exp^- Ja^Qd^j for all (t,x) £ J xR 3 

*b 

of the algebraic reducible differential system (2.10). 

Using the eigenvector v 2 = (1, z, — i) and the corresponding complex-valued eigenfunction 
A 2 : t — > a 2 (t) + a 3 (t) i for all t G J, we can construct (by Theorem 2.3) the first integrals of 
the algebraic reducible differential system (2.10) 

t 

F 2 :(t,x)-> (x? + (x 2 -x 3 ) 2 )exp^-2 Ja 2 (QdCj for all (t,x)eJxR 3 

and 

t 

F 3 : (t,x) ->■ arctan- 2 - * - / a 3 (() d( for all (t,x) G J x X, 

where X is a domain from the set {x: x 1 / 0} C K 3 . 

The functionally independent first integrals F l5 F 2 , and F 3 are an integral basis on the 
domain J x X of the algebraic reducible differential system (2.10). 
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2.2. Triangular systems 

Consider an triangular linear nonhomogeneous differential system of the n-th order 
j n 

- 1 ± = J2a ij (t)x j + f i (t), i = l,...,n, (2.11) 

j=i 

where a „- : J — > R and f ,■ : J — > R are continuous functions on an interval J C R. 
The corresponding homogeneous system of the nonhomogeneous system (2.11) is 

By Perron's theorem [56], any linear differential system can be transformed to an upper 
triangular system by an orthogonal transformation 1 . A triangular system is integrated in 
closed form and its fundamental matrix can be chosen triangular [8, pp. 32 - 33]. 

Using Theorem 2.4 and Corollary 2.2, we can build integral bases for the triangular dif- 
ferential systems (2.11) and (2.12). We now state the main results of this subsection. 

Theorem 2.4. First integrals of system (2.11) are the functions 

T- 1 

F T : (t,x) -> X n+1 „ r ^ n+1 „ r (t) - J2 A r^) F ^ X )- B r(t) 

(2.13) 

for all (t, x) G J x R n , r = 1, . . . , n, 

where the scalar functions 



1 r-i 



A -t 



i v k=i 7 



(2.14) 

£ = 1, . . . , t — 1, t = 1, . . . , n, : i — > 1, £ = 1, . . . , n — 1, Aio : i — ► for all t G J, 



/ (/ n+ i-,(o+E a n + i-,» + i-^(ot + i-^(o # _^oV +1 _ r (o^c, 



T 



</5 r : t ->• exp^ - J a TT {Qd(^j, Yv : * ~> ex P ^ a rr (C)<^C /or a// t £ J, r = 1, . . . ,n. (2.15) 



Proof. The system (2.11) is induced the linear differential operator of first order 

n n 

£(t, s) = c\ + ( Oij(t) + /»(*)) 9* for all (t, x) £ J x R n . 



If r = 1, then the scalar function F x : (t, x) — > x n ip n (t) — B 1 (t) for all (t, x) £ J x R n is 
a first integral of the differential system (2.11): 

£F x (t,x) = x n d t (p n {t) + (a nn (t)x n +f n (t))d Xn (x n ip n (t)) - d t Bx(t) = - a nn (t)ip n (t)x n + 



1 The transformation x = U(t) y is called orthogonal if the matrix U: t —¥ U(t) for all t £ J is orthogonal, 
i.e., U(t)U T (t) = U T {t)U(t) = E for all i€ J, where E is the identity matrix. 
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+ {(a nn (t)x n + f n (t))ip n (t)-d t jf n (() i p n (()d( = for all (t,x) e J xR n . 
Let r = 2. Then the Lie derivative of function F 2 by virtue of system (2.11) is equal to 

ZF 2 (t,x) = £(*„_! P„_i(*) " ^(t^ftz) - B 2^)) = 
= ~ Vl,„-l(^n-lWVl + (Vl,n-l(*)Vl + °n-l,n(*) X n + /n-l(*)) ¥> n _i(*) " 



re— 



t 

^(OWO^CR)^ 



o 



= Vl,n (*)^n (*)V„_i (*) (*n¥>n (*) " ^ (*) " ^ (*, *)) = for all (t, x) G J X M n . 

Therefore the function F 2 : J X R n — > R is a first integral of system (2.11). 
Suppose the functions F T : JxR" — >■ R, r = 1, ...,£— 1, are first integrals of the 
triangular differential system (2.11). Then, for r = e, we have 

£F e (t,x) = £x , 09 , (t) + x , £u> (t) - 

eV ' ' n+l-e ^ n+l-e y ' n+l-e ^n+l-e^ > 

e-1 e-1 

5=1 5=1 

= ( Y, a n+l-e,n+l-eH (*) X n+l- C -K + (*)Vn+l- e (*) " a n+l- £ ,n+l- e (*) X n+l-^n+l- e (*)- 

V 5=0 7 

5=1 v fc=l 7 
n+1— e,ri+l— e+5 

v 5=1 7 

(£-1 £-1 x :-5 \ 

E a n+l-£,n+l -e+^K+l-e+Z " E (<Vl-£,n + l-£+5 W t+1 -£+5^E A e-^ F k^ ^J" 
5=1 5=1 v fc=i 7 

£-1 V 

5=1 7 

£-1 , £-5-1 

C— 1 V ! 1 



5=1 v fc=i 

- B e _ € (i) - A _^(t)F £ _^t,x)\<p n+1 _ E (t) = for all (t,x) G J x R n . 

This yields that if r = e, then F £ : J x R™ — ► R is a first integral of system (2.11). 
Thus the functions (2.13) are functionally independent first integrals of system (2.11). 
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Example 2.4. The linear nonliomogeneous differential system 

^ = -\ x i + -y( 6 ~ e "*) x 2 - (6 + e~')x 3 + 8t 5 + 4t + 2(t 3 + 3t 2 + 6/ + 6)e~\ 

= \ (l + \ - 5e<)x 2 + \[}-\ + 5e*)x 3 + 2^2 (e« - t 3 ), (2.16) 



^ = i(l-^-5e^ 2 + i(l + ^ + 5e*)x, + 2V2 ( 
can be reduced by the orthogonal transformation 



»i = yi. ^2 = ^(2/2-2/3)' ^3 = ^(2/2 + 2/3) 



to the triangular differential system 



^ = " J 2/i - e~*y 2 - 6y 3 + 8t 5 + At + 2(t 3 + 3t 2 + 6t + 6)e" 
^ = 2/ 2 + 5e<y 3 - 2t 3 + 4e*, ^ = J 2/ 3 + 2t 3 . 



(2.17) 



By Theorem 2.4, using the functions 
ip l : t — > t, ip 2 : t — >■ e - *, </j 3 : i — > ^ > V'i : t — > — , ip 2 : t ^ e f , ^3 : * ~~ ^ ^ 2 f° r an * £ ^ 

A n :t->1, A 21 :^^ 3 , A 31 : t^I(9-2t)t 4 , A 32 :t->-^ 2 for all t & J, 

B^.t^t 2 , B 2 : t -> t 5 + 4t + 2(t 3 + 3t 2 + 6t + 6)e~* , B 3 : t -»• (t - 1) t 6 for all i G J, 

where J is an interval from the set {t: t 7^ 0}, we can build the functionally independent 
first integrals on the domain Jxl 3 of the triangular differential system (2.17): 

F 1 : (t, y) -> 1 y 3 - t 2 , F 2 : (t, y) -> e -'y 2 - | ty 3 + 1 i 5 - At - 2(t 3 + 3t 2 + 6t + 6)e"*, 
F 3 : (t, y) ^2/i + ^ * 2 e"*y 2 + ^ (3 - t)t 2 y 3 - - t 3 (2t 4 + t 3 + 4) - t 2 (i 3 + 3t 2 + 6t + 6) e -*. 

Now using the inverse transformation y 1 = x 1 , y 2 = (x 2 + x 3 ), y 3 = (x 3 — x 2 ), we 
obtain the first integrals of the differential system (2.16): 
~ \/2 

F 1 : (t, x) ->• ^ ( x 3 ~ x 2) ~ t2 for a11 (*> x ) G J x r3 > 
F 2 : (i, x )-> (5t + 3e"')x 2 - ^- (5t - 3e"*)^3 + - t 5 - 4t - 2(t 3 + 3t 2 + 6t + 6)e -t , 

Do 

~ \/2 \/2 

F 3 : (t, x) -»• te x + ^- t 2 (t - 3 + e~*)x 2 + ^- t 2 (3 - t + e~*)x 3 - 

- - t 3 {2t A + t 3 + 4) - t 2 (t 3 + 3t 2 + 6/ + 6)e~* for all (t, x) G J x M 3 , J C {t : t / 0}. 

The functionally independent first integrals F±, F 2 , and F3 are an integral basis on any 
domain JxM 3 of the linear differential system (2.16). 
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From Theorem 2.4 with f T (t) = for all t G J, r = 1, . . . , n, we have the following 
Corollary 2.2. First integrals of system (2.12) are the functions 

T-l 

F T : (t, x) -> x n+1 _ T <p n+1 _ T (t) - £ A ri (t)i^ (t, x) for all (t,x) G J x R», r = 1, . . . , n, 

5=1 

where the functions : J — > R, £ = 1, . . . , r, r = 1, . . . , n, anc? 99 T : J — >■ R, r = 1, . . . , n, 
are given by the formulas (2.14) and (2.15), respectively. 

Example 2.5. The second-order linear homogeneous differential system 

= (cos 2 t a{t) — costsint &(t) + sin 2 t c(t)) x 1 + (cost sin t (a (t) — c(t)) + cos 2 1 b(t) — l) x 2 , 

(2.18) 

= (cost sin t(a(t) - c{t)) - sin 2 t b(t) + 1) x 1 + (sin 2 t a(t) + cost sin t b(t) + cos 2 t c(t)) x 2 , 

where a: J — > R, 6: J — > R, and c: J — >■ R are continuous functions on an interval J C R, 
can be reduced by the orthogonal transformation x 1 = cos t y 1 — sin t y 2 , ^2 = srn * 2/i + cos * 2/2 
to the triangular differential system 

^-=a(t)y 1+ b(t)y 2 , = c(t)y 2 . (2.19) 

By Corollary 2.2, the scalar functions 

t 

F\ - (t,y 1 ,y 2 ) -> y 2 exp ^ - J c(() d(^J for all (t,y 1 ,y 2 ) G J x R 2 

*o 

and 

t t C 

F 2 : (t, yi ,y 2 )^ y ieX p(- Ja(QdCj - J b(() exp^J (c(9) - a(6))de^ d( F^y^y^ 

for all (t, yi ,y 2 ) G J x R 2 , t , C G J, 

are first integrals of the triangular differential system (2.19). 

Using the inverse transformation y 1 = costx 1 + sintx 2 , y 2 = — sintx 1 + costx 2 , we 
get the integral basis of the linear differential system (2.18): 

t 

? ' :(Ml ' IlH (c-t*,--«x l) «p(-/«(0««c) forall h,« ! ) 6J xl', 





t 



F a : (*,«,.«,)-» (cos t * 1+S i„ 4 x 2 )ex P ( -/«««)- 
t c 
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2.3. The Lappo-Danilevskii systems 
2.3.1. Linear homogeneous differential system 

Consider an nonautonomous linear homogeneous differential system 

dx m 

= A(t) x, x G R n , A: t -> aji^Aj for all t G J, (2.20) 

where continuous functions ay. J — > R, j = 1, . . . , m, are lineary independent on an interval 
Jet, and Aj, j = 1, . . . , m, are a commuting family of real constant n x n matrices: 

A j A k = A k A j' j = 1 >---,™, k = l,...,m. 

The differential system (2.20) is a Lappo-Danilevskii system [8, pp. 34 - 36], i.e. the 
coefficient matrix A of system (2.20) is commutative with its integral [15]: 

t t 

A(t) j A{t) dr = j A(t) dr A(t) for all t G J, t G J. 

to to 

In this subsection we study the problem of building first integrals for the Lappo-Danilevskii 
differential system (2.20). Using the approaches of constructing first integrals for linear partial 
differential systems [25], we also solve the N.P. Erugin problem of the existence of autonomous 
first integrals [54, p. 469] for the Lappo-Danilevskii system (2.20). 

Partial integrals. A complex-valued linear homogeneous function 

p-.x^rvx for all x G K™, v G C™, (2.21) 

is a partial integral of the system (2.20) if and only if 

2lp(x) = X(t)p(x) for all (t,x) eJxK", (2.22) 

where the linear differential operator 2l(i, x) = d t + A(t) xd x for all (t, x) G J x M. n , the 
function A: J — > C. The identity (2.22) is equivalent to the linear homogeneous system 

( £ a A t ) B j - v = ° f ° r an 1 e j > ( 2 - 23 ) 

V j=1 / 

where E is the identity matrix of order n, the matrices B- are the transpose of the matrices 
Aj, j = 1, . . . , m, respectively. Since A-, j = 1, . . . , m, are a commuting family of matrices, 
we see that there exists a relation [60, pp. 203 - 207] between common eigenvectors and 
eigenvalues of the matrices j = 1, . . . , m. 

The following basic statements (Lemmas 2.3 and 2.4) are base for the method of building 
first integrals of the Lappo-Danilevskii differential system (2.20). The proof of Lemmas 2.3 
and 2.4 is similar to that one in Lemmas 2.1 and 2.2. 

Lemma 2.3. Suppose v is a common real eigenvector of the matrices B- corresponding 
to the eigenvalues X J , j = 1, . . . ,m. Then the function (2.21) is a linear partial integral of 
system (2.20), where the scalar function A: J — > M in the identity (2.22) is given by 

m 

A : t -> ^2 X j aj(t) for all t G J. 



3 



=1 



Proof. If v is a common real eigenvector of the matrices B- corresponding to the eigenva- 

m 

lues X 3 , j = 1 , . . . , m, then v is a common real eigenvector of the matrix B : t — ^ Yl a j (i) Bj 

m 

for all t G J corresponding to the eigenfunction A: t — > X 3 aAt) for all t G J. 

3 = 1 
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Therefore v is a solution to the functional system (2.23). This yields that the identity 
(2.22) is satisfied. Consequently the linear function (2.21) is a partial integral of the Lappo- 
Danilevskii homogeneous differential system (2.20). ■ 

Lemma 2.4. Suppose v = v + v i (v = Re v, v = lmv / 0) is a common complex eigenvec- 

* . ~ . * . . ~ . 

tor of the matrices B- corresponding to the eigenvalues X 3 =\ 3 +\ 3 i (X J = Re A J , X 3 = ImA J ), 
j = 1, . . . , m. Then the Lie derivatives of the scalar functions 

P: x — > (vx) 2 + (vx) 2 for all x G W 1 and -ip: x — > arctan — for all i£l 

vx 

by virtue of system (2.20) are equal to 

m 

2lP(x) = 2^X j aj{t)P{x) for all (t,x) £ JxR™ 

3=1 

and 

m 

2t^(x) = ^ "jW f° r aU (*' x ) 6 J x X ' 

where X is a domain from the set {x: vx / 0} C M n . 

Proof. Formally using Lemma 2.3, we obtain the complex-valued function (2.21) is a par- 
tial integral of the differential system (2.20) and the following identity holds 

%(vx + i Vx) = (*X(t) + i X(t)) {yx + i vx) for all (t, x) S J x 1". 
This complex identity is equivalent to the real system of identities 

%vx = X(t) vx - X(t)vx, 21 vx = X(t) vx + X(t) vx for all (t, x) G J x R n , 

* m * „ m _ 

where the scalar functions A: t — > Yl X 3 aAt) for all t G J, A: t — >■ ]P X 3 uAt) for all i G J. 
Using this system of identities, we have 

2lP(x) = 2t((^x) 2 + (vx) 2 ) = 2vx%Lvx + 2vx$lvx = 2vx (x(t) vx - X(t)vx) + 

+ 2Dx (X{t)vx + X(t)vx) = 2*X(t) ((vx) 2 + (Vx) 2 ) = 2*X(t)P(x) for all (t,x) eJxK", 

, i^x ( A(t) z?x + X(t) vx) - vx ( \(t) vx — X(t) vx) ~, s 
21 ^(x) = K — V V ' — ^ = X(t) for all (t,x) G J x X. ■ 

(ra) 2 + (z/x) 2 

Nonautonomous first integrals. Using Theorems 2.5, 2.6, and 2.7, we can obtain first 
integrals of the the Lappo-Danilevskii homogeneous differential system (2.20). 

Theorem 2.5. Let the assumptions of Lemma 2.3 hold, then the scalar function 

F: (t, x) -»• vxexp(^- J ^ A J ' q^t) dr^j for all (t, x) G J x W a , t G J, (2.24) 

is a /irst integral on the domain J x M n o/ £/te Lappo-Danilevskii system (2.20). 
Proof. From Lemma 2.3, we get 

SLF(t,x) =exp( - ! YX 3 aj (T)dT \ 21 vx + vx 21 exp ( - f ^X 3 a^dr) = 
t 3=1 I 3=1 

m K m 

= ^ X 3 aj (t) F(t, x) - F(t, x) d I ^ X 3 a, (r) dr = for all (t, x) G J x R n . ■ 

i=i / j=i 
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Example 2.6. The linear homogeneous differential system 



dx 1 
~~dt 



= (t + 2 sin t ) x l + sin t x 2 , 



dx 2 
~dt 



= s'mt x 1 + t x 2 



(2.25) 



has the coefficient matrix A: t 



t + 2 sin t sin t 
sini t 



t € 1, where the constant matrices A 1 



such that A(t) = tA 1 + sint A 2 for all 



. Since the matrices 



1 




2 1 


and A 2 = 




1 




1 



A 1 and ^4 2 are commutative, we see that the system (2.25) is a Lappo-Danilevskii system. 

The matrices B x = A\ = A 1 and S 2 = A 2 = A 2 have the common real eigenvectors 
v l = (l — a/2, l) and u 2 = (l + \/2, l) corresponding to the eigenvalues A} = 1, \\ = 1— \/2 
and A 2 = 1, A 2 = 1 + a/2 , respectively. 

Using the functions a 1 : t — > t for all t £ I and a 2 : t — > sin t for all i 6 I, we can 
build (by Theorem 2.5) the first integrals of the Lappo-Danilevskii system (2.25) 

/ t 2 \ 
F 1 : (t,x 1 ,x 2 ) ->• ((1 - y/2)x 1 + x 2 ) exp(^ - — + (1 - y/2)costj for all (t, x x , x 2 ) G R 3 

and 

ft 2 \ 
F 2 : (t,x 1 ,x 2 ) ->• ((1 + \/2)xj +x 2 ) exp(^ - — + (1 + y/2)costj for all (t, x v x 2 ) G M 3 . 

The functionally independent first integrals and are an integral basis of the Lappo- 
Danilevskii differential system (2.25) on space R 3 . 

Theorem 2.6. Let the assumptions of Lemma 2.4 hold, then first integrals of the Lappo- 
Danilevskii differential system (2.20) are the scalar functions 

/ \ m \ 
F 1 : (t,x) -»• ((^) 2 + (ra) 2 ) expf -2 / ^ a j( T ) dT ) f or dl i^ x ) ^ xR " 



and 



F 2 : (t, x) — > arctan 



ux 



vx 



\ m 



r)dr for all (t,x) G J x X, 



where t is a fixed point from the interval J, and a domain X from the set jx: w/O}. 
Proof. Taking into account Lemma 2.4, we obtain 

%F 1 (t,x) = exp(-2 f a^dT^^txf + {vx) 2 ) + 

^ t i =1 

-2 / ^Va^rjdr) = 2 A j ^(i) ^(t, x) + 



+ ((^) 2 + (^) 2 ) d t exp( - 2 / ^ i j a 3 (T)dr] = for all (t,x) G J 
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%F 2 {t,x) = 21 arctan 



ux 



ux 



I m 

21 / ^A'a^r) 
i J=1 



dr 



m „ m 

= ^2\ j aj (t) - d t / p J a 3 (r)(lT = for all (t,x) 6 Jxl 

Thus the functions F 1 and F 2 are first integrals of system (2.20). ■ 
Example 2.7. Consider the second order Lappo-Danilevskii differential system 



dx. 



acosujtx 1 + bsmutx 2 , 



dXr, 



d sin ojtx 1 + a cos uit x 2 , 



where a, b, and d are positive constants, and uj is an non-zero real number. 



(2.26) 



The matrices B 



a 




-d 




and B 2 = 




a 




b 



have the common complex eigenvectors 



v 1 = (\fd, Vbi) and u 2 = (Vd, — Vbi) associated with the eigenvalues \\ = a \f = — Vbdi 

and A3, = a, X 2 = Vbdi, respectively. 

Using the numbers A* = a, A? = 0, X\ = 0, Af = — Vbd, and the scalar functions 
ct-y '. t — r cos ujt, a 2 : t — > smut for all tel, we can find (by Theorem 2.6) the basis of first 
integrals for the Lappo-Danilevskii differential system (2.26) 

2a 

UJ 



F 1 : (t, x 1 , x 2 ) — > (dx\ + bx 2 ) expf sinwi") for all (t, x 1: x 2 ) G M. 3 , 



F 2 : (t,x 1 ,x 2 ) — > arctan 



Vb 



.1 ■) 



>bd 



V~dx l 



cos ujt for all (t,x 1 ,x 2 ) G fl 



on any domain Q, from the set {(i, x^, x 2 ) : x 1 ^ 0} C M 3 . 
Example 2.8. The linear homogeneous differential system 

*^p- = 3(tanhi + t 2 ) x 1 + tanht (x 2 - x 3 ), ^j 2 - = - 2tanhtx 1 + 3t 2 x 2 , 



dx 3 
~dt 



(2.27) 



= tanhi (x l + x 2 ) + (tanhi + 3t 2 ) x 3 



is the Lappo-Danilevskii system of form (2.20) with the constant matrices 





1 










- 2 


2 


- 1 


B l = 





1 





and B 2 = 


- 1 


1 


- 1 










1 




1 









and the scalar functions a ± : t — > tanh t + 3t 2 for all t G M, a 2 : t — ) 
Using the linearly independent common eigenvectors 

^ = (1,0,-1), » 2 = (i,i,l), u 3 = (-i,-i,l) 

of the matrices B 1 , B 2 and the corresponding eigenvalues 



tanhi for all t G 



Ai = l, A? 



1. 



A 



1, A 



A3 



1, A| 



we can construct the integral basis of the Lappo-Danilevskii differential system (2.27) 
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F 1 : (t, x) -»• 3 1 (x 1 - x 3 ) for all (t, x) G R 4 (by Theorem 2.5), 
e* cosh i 

F 2 : (t, x) -> — 3-^ — — (x§ + (x x + x 2 ) 2 ) for all (t, x) G R 4 (by Theorem 2.6), 
e 2 ' cosh t 

F 3 : (t,x) — >■ arctan— !— — — + In cosh i for all (t,x) G (by Theorem 2.6) 

x 3 

on any domain $7 from the set {(t, x) : x 3 / 0} C R 4 . 
Lemma 2.5. Lei i/ie following conditions hold: 

(i) z/° is a common eigenvector of Bj corresponding to the eigenvalues A- 7 , j = 1, . . . , m; 

(ii) i/, = l,...,s— 1, are generalized eigenvectors of the matrix corresponding to 

the eigenvalue with elementary divisor of multiplicity s ^ 2; 

fl?x 

(iii) i/te linear differential system — = A^x has no the first integrals of the form 

Ff g :x^ aj ¥ e {x) forallxeX, j = 1, . . . , m, j ^ (, = l,...,s-l. (2.28) 
Then, we claim that 

a c *f(x) = 1 for all x G X, a f ^(x) =0 for all x G X, = 2, . . . , s - 1, (2.29) 

and 

a -^(x) = /i^ = const /or a// x G X, j = 1, . . . , m, j / (, 8 = 1, . . . , s — 1, (2.30) 

where the linear differential operators a^(x) = AjX d x for all x G M n , f/ie scalar functions 
: X — >■ R are i/ie solution to the functional system 



^2C p l\)*p( x ) v0 ~ Px forallxGX, 0=l,...,s-l, lc{i:A/0}. (2.31) 



Proof. From the proof of Theorem 1.8 it follows that the identities (2.29) for the linear 

dx 

differential system — = A^x are true. Since the linear differential operators of first order 

dx 

a.j, j = 1, . . . , m, are commutative and the system — = A^x hasn't the first integrals (2.28), 

we see that the system of identities (2.30) is satisfied. 

Thus there exist the functionally independent functions : X — >M., 6 = 1, . . . , s— 1 (the 
solution to the system (2.31)) and these functions satisfies (2.29) and (2.30). I 

Theorem 2.7. Under the conditions of Lemma 2.5, we get first integrals of the Lappo- 
Danilevskii differential system (2.20) are the scalar functions 

t „ 



= !,...,*- 1, (2.32) 



/no 
fJ- j e C dr for all (t, x) G J x X, 9 = 1 

where t is a fixed point from the interval J, a domain X from the set jx: v°x / 0} C R n . 
Proof. From the identities (2.29) and (2.30), on the domain J x X we get the following 

\. rn m 

%F e (t, x) = - d t / 4 C a j( T ) dT + a ^U x ) = °> e = i,..., s -i. 

Therefore the functions (2.32) are functionally independent first integrals on the domain 
J x X of the Lappo-Danilevskii differential system (2.20). I 
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Example 2.9. Lappo-Danilevskii differential systems with non-diagonal coefficient mat- 
rices of the second order are the linear differential systems of the form [58] 

dec dx 

—± = (a l (t) + b l a 2 (t))x l +a 2 (t)x 2 , —j2- = b 2 a 2 (t)x 1 + a 1 (t)x 2 , (2.33) 

where the functions a 1 : J — > R and a 2 : J — > R are continuous, b 1 and b 2 are real numbers. 
Consider the number D = b\ + 4b 2 . We have three possible cases. 

Let D > 0. Then, using the common real eigenvectors v 1 = (1, — A 2 ), v 2 = (1, — A 2 ) 
and the corresponding eigenvalues A} = f , A 2 = (b\ — v^D) /2, A 2 = 1 , A 2 = (61 + V^D) /2, 
we can build (by Theorem 2.5) the integral basis of system (2.33) 

t 

F k : {t,x 1: x 2 ) -> (x x - \l_ k x 2 )exj>(^ - J (a^r) + \ 2 k a 2 (r)) drj 

for all (t, xi,x 2 ) G J x R 2 , i G J, A; = 1,2. 

Let D < 0. Then, using the common complex eigenvector u 1 = (Af , 1) and the corre- 
ct; * ~ 
sponding eigenvalues A} = 1, A 2 = A - Ai, where A = &i/2, A = V — D/2, we can construct 

(by Theorem 2.6) the basis of first integrals for system (2.33) 

t 

F x : (t,x 1} x 2 ) -»• ((*\x 1 + x 2 ) 2 + (\ Xl ) 2 ) exp(- 2 j (a^r) + Aa 2 (r)) dr^) 

for all (t,x 1 ,x 2 ) G J x R 2 , t Q G J, 
Ax /*- 

F2 : (t, x 1; x 2 ) — > arctan — + / Xa 2 (r)dT for all (t, x 1 , x 2 ) G J x X, 

A x x + x 2 ^ 

where X is a domain from the set {(xi,X2) : A^j +s 2 / 0} C R 2 . 

Let D = 0. Then, using the common real eigenvector v° = (A 2 , 1), the 1-st order real 
generalized eigenvector u 1 = (1,0), and the corresponding eigenvalue A 2 = b 1 /2 with ele- 
mentary divisor of multiplicity s = 2, we can find (by Theorems 2.5 and 2.7) the functionally 
independent first integrals of system (2.33) 

t 



x R 2 , 



F 1 : (t, Xi, x 2 ) — > (A 2 x 1 + x 2 ) exp^ — J (a x (r) + A 2 a 2 (r)) dr^j for all (t, x 1: x 2 ) & J 

( 

F 2 : (t,x 1 ,x 2 ) — > 2 Xl / a 2 (r)dr for all (t,x 1 ,x 2 ) G J x X, t G J, 

1 "^2 ^ 

where X is a domain from the set {(x 1 , x 2 ) : \\ x 1 + x 2 / 0} C R 2 . 

in the complex case, from the complex-valued first integrals (2.32) of the Lappo-Danilev- 
skii differential system (2.20), we obtain the real first integrals 

* TO 

: (t, x) -> Re*J(x) - / J^Re^a^r)^ for all (t,x)eJxI, 9 = 1,. ..,8-1, 



* TO 



/Iff 
^ Im nf aj[r) dr for all (t, x) G J x X, = 1, . . . , s - 1. 
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Autonomous first integrals. The base of building autonomous first integrals for the Lappo- 
Danilevskii differential system (2.20) is the following basic proposition. 

Lemma 2.6. A function F: X — > R is an autonomous first integral of the Lappo-Da- 
nilevskii differential system (2.20) if and only if this function is a first integral of the linear 
homogeneous system of partial differential equations 

a j (x)y = 0, j = l,...,m, (2.34) 

where the linear differential operators of first order dj(x) = A-xd x for all x G R n . 

Proof [54, p. 470]. A function F: X — > R is an autonomous first integral of the Lappo- 
Danilevskii differential system (2.20) if and only if the following identity hold 

m 

(*j(t) A j x d x F{x) = for all (t, x) G J x R n . 

3=1 

Since the functions a^: J — > R, j = 1, . . . , m, are linearly independent on the interval J, 
we see that this identity is equivalent to the system of identities 

Aj x d x F(x) = for all x G R™, j = 1, . . . , m. 

This yields that the function F: X — > R is a first integral of system (2.34). ■ 

By Theorems 2.8, 2.9, 2.10, and 2.11, we can construct autonomous first integrals of the 
Lappo-Danilevskii homogeneous differential system (2.20). 

Theorem 2.8. Suppose v k are common real eigenvectors of the matrices B- corres- 
ponding to the eigenvalues A^, j = 1, . . . ,m, fc = l,...,m+l. Then the Lappo-Danilevskii 
differential system (2.20) has the autonomous first integral 

m+l 

F:x^ n\v k x\ hk forallxGX, IcDF, (2.35) 
fc=i 

where the real numbers h k , k = 1, . . . ,m + 1, are an nontrivial solution to the system 

m+l 

Yl X i h k = °> j = l,...,m. 
k=l 

Proof. First note that the Lappo-Danilevskii system (2.20) is induced the linear differential 

dx 

systems — = A-x with the operators a. Ax) = A-xd x for all x G R n , j = 1, . . . ,m. 
at J j j 

If v k are common real eigenvectors of the matrices Bj corresponding to the eigenvalues 
X 3 k , j = 1, . . . ,m, fc=l,...,m+l, then the linear homogeneous functions (by Lemma 1.1) 
p k : x — > v k x for all x G R n , fe=l,...,m + l, are partial integrals of the linear autonomous 



dx 

differential systems — = A-x, j = 1, . . . ,m, and the following system of identities hold 



a jPk {x) = \{p k (x) forallxGR™, j = l,...,m, k = 1, . . . , m + 1. (2.36) 

We obviously have 

m+l m+l m+l 

a j F(x)= Yl\p k (x)\ s & n Pk( x ) h k II \Pi( x )\ a jPk( x ) fora11 x£X, j = l,...,m. 

k=l k=l l=l,l^k 

Now taking into account the identities (2.36), we obtain 

m+l 

(Xj F(x) = Y2 tfc ^k F(x) for all x G X, j = 1, . . . , m. 
k=i 
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If the real numbers h k , fc=l,...,m+l, are an nontrivial solution to the linear homoge- 

m+l 

neous system Yl ^k ^k = ^' 3 ' = 1> • • • ' m -> then the scalar function (2.35) is a first integral 
k=i 

for the linear homogeneous system of partial differential equations (2.34). 

By Lemma 2.6, the function (2.35) is an autonomous first integral of system (2.20). ■ 

Example 2.10. Consider the third-order Lappo-Danilevskii differential system 



dx 1 
~~dt 



e* cos e* x 1 + 2 sinh t (x 2 + x 3 ) , 



dx 2 



sinh t x 1 + e* cos e* x 3 , 



^3 
(ft 



(2.37) 



sinh i x x + e* cos e* x 2 . 








1 


1 




1 








The matrices S x = 


2 








and B 2 = 








1 




2 













1 






have the common real eigenvectors 



= (1, — 1, — 1), v 2 = (0, 1, — 1), and = (1, 1, 1) corresponding to the real eigenvalues 
A{ = - 2, A? = 1, A 2 = 0, \l = - 1, and A^ = 2, A§ = 1, respectively. 

From the linear homogeneous system — 2h 1 + 2h 3 = 0, h 1 — h 2 + h 3 = 0, we have, for 
example, h 1 = 1, h 2 = 2, h 3 = 1. By Theorem 2.8, the scalar function 

F: (x 1 ,x 2 ,x 3 ) ->• (x\ - (x 2 - x 3 ) 2 )(x 2 - x 3 ) 2 for all (x 1 ,x 2 ,x 3 ) G IR 3 

is an autonomous first integral on space KL 3 of the Lappo-Danilevskii system (2.37). 

By Theorem 2.5, using the functions a 1 : t — > sinht and a 2 : t — > e t cose t for all t G R, 
we can build the nonautonomous first integrals of system (2.37) 



2 cosh t— sin e 



F 1 : (t, x-y, x 2 , x 3 ) — > e 
and 

F 3 : (t, x 1 , x 2 , x 3 ) —7- e 



(xi x 2 x 3 ), F 2 : (t,Xi,x 2 ,x 3 ) —7- e (x 2 x 3 ), 

2cosht— sine 4 , . „ . . A 

(x l + x 2 +x 3 ) for ail [t, x 1 ,x 2 ,x 3 ) G K . 



Thus every set of the functionally independent first integrals {F, F 1: F 2 }, {F, F ± , F 3 }, 
{F, F 2 , F 3 }, and {F 1: F 2 , F 3 } is an integral basis on the space IR 4 of system (2.37). 

Corollary 2.3. Let v k be real common eigenvectors of the matrices B- corresponding to 
the eigenvalues \ 3 k , j = 1, . . . , m, k = 1, . . . , m + 1. Then an autonomous first integral of the 
Lappo-Danilevskii differential system (2.20) is the scalar function 



F 12...m(m+1) ! X ~^ 



n 

k=l 



v k x\ 



\v m+1 x\ 



-A 



for all x G X, 



where the determinants A^, k = 1, . . . ,m are obtained by replacing the k-th column of the 
determinant A = lA^I by colon ( A^ +1 , A™ +1 ) , respectively. 

Example 2.11. Consider the fourth order Lappo-Danilevskii differential system 
dx 



dt 



= Q; 1 (t)A 1 x + a 2 (t)A 2 x, x G 



with linearly independent continuous functions a x : J 



matrices A 




1 

2 
2 



6 

- 7 

- 6 
10 



1 
1 

3 
2 



A, 



I, a 2 : J - 

1 4 
1-4 1-3 
0-4 1-2 

2 6-25 



(2.38) 

, and the commutative 
2 
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The matrices B 1 = Af and B 2 = A 2 , where T denotes the matrix transpose, have the 
linearly independent common real eigenvectors 



v 1 = (0,2,0,1) 



(2,2,1,1), v* = (1,0,1,0), ^ 4 = (-l,l, -1,1) 



corresponding to the real eigenvalues 

A} = — 2, A3, = 1, A3 = 2, A4 = 3, Af 
The determinants 



— 1, A| — 1, A3 — 1, A4 — 2. 



A = 



- 2 1 

- 1 1 

A 



= -1, 



A 



11 



2 1 
1 1 



1, 



A 2 i = 



-22 
- 1 1 



= 0, 



12 



3 1 
2 1 



A 



22 



2 3 
1 2 



By Corollary 2.3, the scalar functions 

F 123 : x — > (2x 2 + x A )(x 1 + £ 3 ) for all x G 



and 



F l24 : x 



(2x 2 ~\~ x^)(y x-^ + x 2 £3 + x^j 



for all 



where X is a domain from the set {x : 2x± + 1x 2 + x 3 + x 4 7^ 0} C R 4 , are autonomous first 
integrals of the Lappo-Danilevskii differential system (2.38). 

Theorem 2.9. Suppose v k = u k + u k i (v k = Rei' k , v k = lm.v k ), k = 1, . . . , s, s ^ (m+l)/2 
(this set hasn't complex conjugate vectors), and u , 9 = s + 1, . . . , m + 1 — s, are common 
complex and real linearly independent eigenvectors of the matrices Bj, j = 1, . . . , m. Then the 
Lappo-Danilevskii differential system (2.20) has the autonomous first integral 



m+l— s 



F: x^ "[\{P k {x)) " exp( -2h k cp k (x)) ] J \v x\ " for all x G X, Xc 



(2.39) 



k=i 



=s+l 



where the scalar functions P k : x — > (v k x) 2 + [v k x) 2 for all x G W 1 , (p k : x — > arctan 



v k x 

V K X 



for all x G X, the real numbers h kl h k , h e are an nontrivial solution to the system 

s m+l— s 

2 E(Hk-~ X iK)+ E K h e=°> 3 = l>->™> (2.40) 



k=l 



e=s+i 



the numbers X k = \ k + \ k i (X k = ReA^, X k = ImA^) and X J e are eigenvalues of the matrices 
Bj corresponding to the eigenvectors u k and v®, respectively. 

Proof. By Lemma 1.1, Properties 1.4 and 1.5, the linear homogeneous functions 

p^ : x — > v^x for all x G M n , £ = 1, . . . , m + 1, 

dx 

are partial integrals of the linear differential systems — = A - x, j = l,...,m, and the 

at J 

following system of identities hold 

a j P k (x) = 2*\{P k (x) forallxGK™, j = l,...,m, k = 1, . . . , s, 
<Xj(p k (x) = \ k for all x G X, j = l,...,m, k = l,...,s, (2-41) 
dj v e x = \ 3 e v e x for all x G W l , j = 1, . . . , m, 9 = s + 1, . . . , m + 1 — s. 
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Since the function F: X — > R is given by (2.39), it follows that 



K- 1 



^ k=l k=l l=l,tyk 

s h s \ m +l — s h 

fc=l fc=l ' 0=s+l 

s J m+l—s ^ y m+l—s m+1— s 

+ u( p M fc ex P (-2^^(x)) n w e «gn(^)^ n i^k-(^) 

fc=l 6»=s+l 6»=s+l J=s+l,J^0 

for all i£X, j = 1, . . . ,m. 
Using the system of identities (2.41), we get 

m+l—s 



S .F(x)= ( 2^(A^ fe -A^J+ J] forallxeX, j = l,..., 

^ fc=l 6»=s+l ' 



m. 



If the numbers h k , h k , k = 1, . . . , s, and h g , = s + 1, . . . , m + 1 — s, are an nontrivial 
real solution to the system (2.40), then (by Lemma 2.6) the scalar function (2.39) is an 
autonomous first integral of the Lappo-Danilevskii differential system (2.20). I 

Example 2.12. Consider the fourth order Lappo-Danilevskii differential system 
dx 



dt 



= a l (t)A l x + a 2 (t)A 2 x, x 6 R , 



(2.42) 



where linearly independent functions a 1 : J — > R and a 2 : J — > R are continuous, the 

and Ar, = 



matrices A 



- 1 


2 


- 2 





- 6 


5 


-4 


2 


-3 


2 


- 2 


1 


2 


- 1 


2 






and 


q 2 : 


J 


-> 


-4 


2 


2 


4 


- 5 


1 


6 


6 





- 1 


3 


1 


-3 


2 





3 



commute. 



The matrices B l = Aj and £? 2 = A 2 ', where T denotes the matrix transpose, have 
the common eigenvectors z^ 1 = (1 + i, — i, — 1 + i, — 1), z/ 2 = (1 — z, i, — 1 — i, — 1), 
v 3 = (1, — 1, 2, 0), ^ 4 = ( — 1, 0, 2, 2) corresponding to the eigenvalues Aj = l + i, A 2 = 1 — i, 
A* = - 1, A 4 = 1, and \\ = i, A 2 . = - i, A§ = 1, A| = 2. 

The linear homogeneous systems 

2h u - 2h u - h 31 = 0, - 2h u + h 31 = and 2fc 12 - 2/i 12 + /i 42 = 0, - 2h 12 + 2/i 42 = 

have, for example, the solutions h\\ = 2, /i n = 1, h 31 = 2, and /i 12 = 1, h l2 = 2, h 42 = 2. 
By Theorem 2.9, the scalar functions 

ry I i~Y* 

1 x 2 ~r x 3 



F x : x — > (x 1 — x 2 + 2x 3 ) ((x x — x 3 — x 4 ) 2 + (x x — x 2 + x 3 ) 2 ) exp ( — arctan 



and 



F 2 : x ->■ ( — x 1 + 2x 3 + 2x 4 ) 2 ((x 1 — x 3 — x^ 2 -!-^ — x 2 +x 3 ) 2 ) expf -4 arctan— — x 2 + x 3 j 

V x : — x 3 — x 4 y 

are autonomous first integrals of the Lappo-Danilevskii differential system (2.42) on any 
domain X from the set {x : x 1 — x 3 — x 4 / 0} C R 4 . 
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Example 2.13 (continuation of Example 2.8). From the linear homogeneous system 
h 1 + 2/i 2 = 0, — h l — 2h 2 = it follows that, for example, h 1 = — 2, hi = h 2 = 1. 

Then, using the common eigenvectors v l = (1, 0, — 1), v 2 = (i, i, 1) and the corresponding 
eigenvalues \\ = \\ = 1, Af = — 1, A 2 , = i, we can build (by Theorem 2.9) the autonomous 
first integral of the Lappo-Danilevskii differential system (2.27) 

(00 I j f j HQ "\ 

F: (xi,io,x,) -> — 7 — — - expl — 2arctan— -] for all (s, , x 9 , x~) € X, 

(X X - X 3 ) 2 V » 3 / 

where X is a domain from the set {(x ± , x 2 , x 3 ) : x 1 / x 3 , x 3 / 0}. 

Theorem 2.10. Suppose v T =u T +V T i, v s+T = v T — D T i, r= 1, . . . , s, s ^ m/2, z/ 2s+1 = 
_ ^2s+i _|_ ^2s+i ^ ffln ^ ^ = 2s + 2, . . . , m + 1, are common complex and real linearly 
independent eigenvectors of the matrices j = 1, . . . , m. Then the Lappo-Danilevskii diffe- 
rential system (2.20) has the autonomous first integrals 



F V .x^ Y[{P T (x)) e W (-2{h T -h s+r )cp T (x))- 

r=l 

(2.43) 



m+1 2h a 



( P 2s+i( x )) ex P v - 2 ^ 2s +i^2 S +i( x )J 11 (" x ) forallxGX 



0=2s+2 



and 



F 2 :x^ f[{P T (x)) T 3+T exp(2(h T -h s+r )^ T (x))- 
r=l 

(2.44) 

ft * m+l 2fc 

• (^+i( x )) 23+1 ex p( 2/i 2,+i ^s+i^)) II K*) * for all x G X, X C DF 1 n DF 2 , 

6»=2s+2 

where the scalar functions 

~ T 

* 2 2 V X 

P T : x — >■ (z/ T x) + (z7 T x) , : x — >■ arctan /or a// x G X, r = 1, 2, . . . , s, 2s + 1, 

v T x 

i/ie numbers h k = h k + h k i, k = 1, . . . , m + 1, are an nontrivial solution to the system 

m+l 

X k h k =°> i = 1, 

fc=l 

and Ai are eigenvalues of the matrices Bj corresponding to the eigenvectors u k . 
Proof. We form two complex-valued functions of real variables 

* 2s h h m+l , 

F-. x -> n(^) fc 2s+i n 9 f ° r a|uei 

fe=l 6»=2s+2 



and 

^ 2s j j m+l j 

F:x^ n(i/ fe x) fc (^+!x) +1 [] (*/x) S forallxGX, 

fe=l 6»=2s+2 

where / fc , A; = 1, . . . , m + 1, are some complex numbers, a domain X C W 1 . 

40 



V.N. Gorbuzov, A.F. Pranevich 



First integrals of ordinary linear differential systems 



We have 



m+l 



ajF(x) 



<XjF(x) = ^2 K h-k F( x ) f° r all x G X, j = 1, . . . , m, 
fc=i 

2s m+l s 

: \ 

A fe * fc + A L+i z 2 s+i + 2 A e l e ) for allieI ' j = l,...,m. 



fc=i 



(9=2s+2 



Let h k = h k + h k i, k = 1, . . . , m + 1, be an nontrivial solution to the linear system 

m+l 

J]^* = l>. j = l,...,m. 



k=i 



Then l k =h s+k -h s+k i, l g+k 



h k~ h k i ' k = 1 ,---, s , 1 
= 2s + 2, . . . ,m + 1 is an nontrivial solution to the linear system 
2s m+l 

E A i^ + A i+i^+i+ A^ = 0, j = l,...,m, 

fc=l 0=2s+2 

and the functions F: X — )■ C, F: X — ^ C are autonomous first integrals of the system (2.20). 

Since F x (x) = F(x) F(x) for all x G X and F 2 (x) = F % (x) F % (x) for all x G X, we 
see that the scalar functions (2.43) and (2.44) are autonomous first integrals (by Lemma 2.6) 
of the Lappo-Danilevskii differential system (2.20). ■ 

Example 2.14. The small oscillations of gyrocompass in the precession gyroscope theory 
are described by the following fourth-order linear differential system [61] 



2 s +1 ^2s+l ^2s+l h 



h n — h n i, 



- = uj q x 2 + Q(t)x 4 , 



dx 
~~dt 



tt(t)x 3 , 



dx 2 
~dt 



(2.45) 



where a continuous on an interval Jet function f2 : J — > R is the projection of the absolute 
angular velocity of the gyroscope's sensitive element onto the direction of geocentric vertical 
line, the constant oj q = y/g/R (g is the acceleration due to gravity, R is the Earth's radius). 
The system (2.45) is a Lappo-Danilevskii differential system [8, pp. 37 - 38]. 
By Theorem 2.10, using the linearly independent common eigenvectors u 1 = — 
(i, u 3 = (l, —i, —1, — i), v A = { — i, —1, —i,l) and the corresponding eigen- 



,2-1 



values \\ = \ 2 = — ujq i, A3 = A4 = ujq i, Xf = — i, X 2 = A3 = i, A 4 



B 1 






-u 








UJ Q 




















u 








-w 






and Br, 







1 

1 



2 = — i of the matrices 

- 1 

1 





we can build the autonomous first integrals of the Lappo-Danilevskii differential system (2.45) 

F 1 : x — > (x 1 — x 3 ) 2 + (x 2 + x 4 ) 2 , F 2 : x — > ( — x 2 + x 4 ) 2 + (x x + x 3 ) 2 for all x G IR 4 . 

By Theorem 2.6, using the functions a 1 : t — > 1 for all t G J, a 2 : i — O(t) for all t £ J, 
the numbers Ai = A? = 0, A{ = — ujq, \\ = — 1, \\ = \\ = \\ = — ujq, A 2 = 1, and the 

common eigenvectors v 1 and v 2 , we can construct the first integrals of system (2.45) 
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F 3 : (t, x) — > arctan 
and 



x 2 



+ J (oj + ri(r)) dr for all (t, x) G JxX l5 X 1 C {x: x 1 — x 3 / 0}, 



z 

F 4 : (t, x) — >■ arctan — — + (u> - fi(r)) for all (i, x) G JxX 2 , X 2 C : x 2 — x 4 / 0}. 

— x 2~ >r X A J 

The functionally independent first integrals F ly . . . ,F 4 are an integral basis of the Lappo- 
Danilevskii system (2.45) on any domain J x X, where X C {x : x 1 — x 3 7^ 0, x 2 — £ 4 7^ 0}. 
Example 2.15. The Lappo-Danilevskii differential system 



dx-\ . * , , dx 2 

— = ai (t) Xl+ a 2 (t)x 3 , — 



dx 3 
dt 



-a 2 (t) Xl +a 1 (t)x 3 , 



dx 4 
~~dt 



= - a^f)^ + a 2 (t)x 4 , 
= - a 2 {t)x 2 - a 1 (t)x 4 



with linealy independent on an interval J C R functions a 1 : J 



and a 2 : J 



(2.46) 



such 



that the matrices B 1 



10 
0-10 
10 
-1 



and B, 



have the linealy 



0-1 

-1 

10 

10 

(0,i,0,l), ^ 3 = (-i,0,l,0), */ = (i,0,l,0) 



independent eigenvectors u 1 = (0, — i,0, 1), V 
corresponding to the eigenvalues A{ = A 2 = — 1, A3 = \\ = 1, Af = A§ = — i, A 2 = A 4 = i. 

From the linear system — h 1 — h 2 + h 3 = 0, — ih 1 + ih 2 — ih 3 = it follows that, for 
example, h 1 = 0, h 2 = h 3 = 1. Then, by Theorem 2.10, the Lappo-Danilevskii differential 
system (2.46) has the autonomous first integrals 



F 1 : x — > (x\ + x 3 )(x 2 + x\) for all x G 



F 2 . x ^ ^2+^4 forallxeX) 



^ 2 ^3 

where X is any domain from the set {x: x x x 4 — x 2 x 3 / 0} C R . 

Example 2.16. Consider the fifth-order Lappo-Danilevskii differential system 

dx 



dt 



= a 1 (t)A l x + a 2 (t)A 2 x + a 3 (t)A 3 x, xG 



(2.47) 



where linearly independent functions a : J — > R, j = 1, 2, 3 are continuous, the matrices 



5, 







4 




2 


2 




7 


11 




4 -5 


7 









2 


2 




7 


7 




4 -5 


5 


AJ = 







4 







2 


2 




B 2 = Al = 


4 









4 




8 - 


8 




6 


- 2 




-2 6-6 6 - 


2 






4 




4 - 


4 




4 







- 2 2 - 2 2 


2 









6 




- 6 




- 9 


17 













6 




18 




3 


- 3 






3 — ^3 









12 




24 




18 - 


- 18 


(T denotes the matrix transpose). 








8 


20 




20 




14 - 


- 14 










8 


8 




- 8 




8 


-8 







Using the linearly independent common eigenvectors v 1 = (1,0, 0, i, i), v 2 = (1,0, 0, —i, —i), 
u 3 = (1 + 2i, 1 + 2i,2, 2,0), v A = (1 - 2i, 1 - 2i, 2, 2, 0), i/ 5 = (0, 1, 1,0,0) of the matrices B l5 
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B 2 ,B^ and the corresponding eigenvalues X\ = 4+4i, A 2 = 4— 4i, A3 = 4i, A 4 = — 4i, Ag = 4, 
\\ = 4 + 2i, A| = 4 - 2i, A§ = 2 + 4i, \\ = 2 - 4i, \j = 4, A? = 8i, Af = - 8i, A| = 12z, 
A 4 = — 12i, A 2 = — 12, we can construct (by Theorem 2.10) the autonomous first integrals 
of the Lappo-Danilevskii differential system (2.47) 

F . x (xf + (x 4 + x 5 ) 2 ) 2 /_ 2 arctan 2(x 1+ x 2 ) A 

1- (x 2 +x 3 ) 2 ((x 1 +x 2 +2x 3 +2x 4 ) 2 + 4(x 1 +x 2 ) 2 ) P \ xi+x 2 + 2x 3 + 2x 4 / 

and 

((x x + x 2 + 2x 3 + 2x 4 ) 2 + 4(x x + x 2 ) 2 ) 5 

F 2 - X -> , , ,10 • 

/ X4 + Xr 2(Xl + x 2 ) 

• exp 12 arctan — 10 arctan 



x-y X-^ ~t~ x 2 ~t~ 2X3 "I - 2x 4 

on any domain X from the set {x : x 2 + x 3 7^ 0, x 1 + x 2 + 2x 3 + 2x 4 7^ 0, x x 7^ 0} C M 5 . 
Theorem 2.11. Lei f/te following conditions hold: 

(i) i/ w are linearly independent common real eigenvectors of the matrices B- correspond- 
ing to the eigenvalues A^, I = 1, . . . , r, j = 1, . . . , m; 

(ii) i/« = 1,..., s t —l, are linearly independent real generalized eigenvectors of the matrix 
B^ corresponding to the eigenvalue A^ with elementary divisor of multiplicity s t , I = 1, . . . ,r, 

r 

such that Yl s i ^ m + 1; 

dx 

(iii) i/ie linear differential system — = A^x has no the first integrals of the form 

F j8l : x ~^ a J ^ei^) forallxGX, j = l,...,m, j / (, 9 = 1, . . . , - 1, l = l,...,r. 
Then autonomous first integrals of the Lappo-Danilevskii system (2.20) are the functions 

k ££ 

F: x —?■ JJ(z/° 5 x) V exp^ #^(x) for all x € X, X C DF, (2.48) 
S=i <?=i 

where the scalar functions Vl/^ : X — >■ R are £/je solution to the functional system 
q 

^x = Y,( q p Zl)^^(x)^-f^xforallxeX, q = l,...,e^ £ = l,...,k, (2.49) 
P =i 

fe 

and ^ = m — k + 1, ^ — 1, £ = 1, . . . , A;, k ^ r. Afeo, Ziere 

S'^( x ) = ' u 5' /«J = const ; i = !,••• g=l,...,££, f = l,...,fc, (2.50) 
£/ie numbers h q ^, q = 0, . . . , e^, £ = 1, . . . , fc, are an nontrivial solution to the system 

k £ i 

E( A IV + E^M=°- J' = l,---,m, (2.51) 
€=1 5=1 

Proof. The functional system (2.49) with any fixed index £ G {1, ... , fc} has the deter- 
minant (^°^x) e? for all x G M n such that (z/°^x) £5 7^ for all x G X, where a domain 
X C {x: z/°£x / 0, £ = 1, . . . , k}. Therefore for any fixed index £ G {1, . . . there exists 
the solution Vl/^ : X — > R, g = 1, . . . , e^, on the domain X of the functional system (2.49). 
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For any fixed index £ G {1, . . . , k} from the results of Lemma 2.5 it follows that the 
identities (2.50) are satisfied. 

k 

Consequently there exist Yl £ f functionally independent on the domain X scalar func- 
tions \lK : X — > R, g = 1, . . . , e^, £ = 1,. . . ,k, such that the identities (2.50) hold. 

Since the functions p^: x — > u°^x for all x G W 1 , £ = l,...,k, are partial integrals of 

dx 

the linear differential systems — = A ■ x, j = 1, . . . , m, we have 

at J 



a j F ( x ) = ( ( X l V + Y Kd ) F ( x ) for all a; G X, j = 1, . . . , m. 



If the numbers h q ^, q = 0, . . . , e^, £ = 1, . . . , A;, are an nontrivial solution to the linear 
homogeneous system (2.51), then the function (2.48) is a first integral of system (2.34). 

By Lemma 2.6, the scalar function (2.48) is autonomous first integral on the domain X 
of the Lappo-Danilevskii differential system (2.20). I 

Example 2.17. Consider the fourth-order Lappo-Danilevskii differential system 



dx 

— =a 1 (t)A l x + a 2 {t)A 2 x, x G R 4 , 



where linearly independent functions a 1 : J 
interval JcK, the commuting matrices 



and a 2 ■ J 



(2.52) 

are continuous on an 



A- 



1 

2 - 1 - 1 

10 0-1 

10 2 2 



and A 



2 





- 1 





- 1 


2 





1 


- 1 





3 


1 





1 


- 3 


1 



The matrix B 1 = A±, where T denotes the matrix transpose, has the eigenvalue \\ = 1 
with elementary divisor (A 1 — l) 4 corresponding to the eigenvector z> 01 = (—1,1, — 1, 0) and 
to the generalized eigenvectors u 11 = (1, 0, - 1, - 1), u 21 = (1, - 1, 3, 0), u 31 = ( - 3, 0, 9, 9). 

From the functional system (2.49), we get the scalar functions 



X i Xq 



x 4 



for all x £ X, 



^ : x -»• 



X-^ ~\~ 

( — x 1 + x 2 — x^)(x 1 — x 2 + 3x 3 ) — (x 1 — £3 — x 4 ) 2 



( X-y -\- 3?2 



X q 



for all 



*3i :x_> 7 A3 

V ~~ x l + x 2 ~ X 3) 



(( — 3x 1 + 9x 3 + 9x 4 )( — x 1 + x 2 — x 3 ) 



3( x 4 + x 2 x 3 )(x 4 x 3 x 4 )(x 1 x 2 + 3x 3 ) + 2(x 4 x 3 



for all x G X, 



where X is any domain from the space {x : x 1 — x 2 + x 3 ^ 0} of the space R 4 . 

The Lappo-Danilevskii differential system (2.52) is induced on the space M 4 the linear 
differential operators of first order 



«i (x) = x 2 d Xi + (2x 2 - x 3 - x 4 ) d X2 + (x 1 - x 4 ) d X:j + ( - x 1 + 2x 3 + 2x 4 ) d x 



and 



2 (x) = (2xj - x 3 ) a a . i + ( - x 4 + 2x 2 + x 4 ) a x2 + ( - x x + 3x 3 + x 4 ) d Xa + (x 2 - 3x 3 + x 4 ) d x ^ . 
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Using the identities 

a 1 u 01 x = u 01 x for all x G M 4 , o 2 v m x = 2v m x for all x G M 4 , 

tti*!!^) = 1, a^^x) = 5 a i^li( x ) = for an ^ G X, 

a 2 V\ 1 (x)=-l, a 2 ^^(x) = 0, a 2 *^(a;) = 6 for all x G X, 

we obtain the linear homogeneous systems 

h u + h 21 = 0, 2/i n — /i 21 = and h 12 + /i 22 = 0, 2h 12 — h 22 + 6/i 32 = 0. 

Prom these systems it follows that h u = h 21 = 0, h 31 = 1, and /i 12 = 2, /i 22 = — 2, h 32 = — 1. 
Then, by Theorem 2.11, the scalar functions 

F 1 : x ->• * 21 (x) and F 2 : x ->• ( - x 1 + x 2 - x 3 ) 2 exp( - 2^\ 1 {x) - *3i(x)) for all x G X 

are autonomous first integrals of the Lappo-Danilevskii differential system (2.52). 

h 



If the matrices B, have some complex common eigenvectors u 01 corresponding to the 



eigenvalues A*} with elementary divisors s t , then the proof of Theorem 2.11 is also true. 
Let the set V of m + 1 functions be given by 

k 

V = {v°*xfoT allx G R n , *^(x) for allx G X: q= l,...,e^ £=l,...,k, J2 £ t = 

5=1 

In the complex case, we shall have two logical possibilities: 

1. Any function from the set V has the complex conjugate function in the set V. 

2. At least one function from the set V has not the complex conjugate function in the V. 
In these cases the system (2.20) has the following autonomous first integrals. 

Case 1. The Lappo-Danilevskii differential system (2.20) has the autonomous first integral 



fcl * ~ Q£ 

F: x —?■ [J ( (^) 2 + ( * **) 2 ) exp ( - 2 ^ arctan + 



«=1 



£ « ki 

* * <■ ~ ~ ^ ■ > ^ " ■ " I hr. 



+ 2 - II l^l ^exp^ h qp *< p (x) for all x G X C DF, 

q=l p=l q=l 

where h q £, h q £, and h qp , q = Q,...,s k , k = £ or k = p, £ = 1,...,^, p=l,...,k 2 , are 
a real nontrivial solution to the linear homogeneous system 

fcl e £ k 2 £ P 

£=1 g=l p=l q=l 

Here z/ ^ = u ^ + v°£ i are complex common eigenvectors of the matrices B- correspond- 
ing to the eigenvalues A| = A| + A| i, j = 1, . . . , m, £ = 1, . . . , fcj^, respectively; z^ 0p are real 
common eigenvectors of B- corresponding to the eigenvalues Xp, j = 1, . . . , m, p = 1, . . . , k 2 , 

respectively; the functions 5^ = yp^ + v]/^ i and are the solution to the functional sys- 
tem (2.49); the real numbers 

h j q \ = Re ( aj *^(x)) , = Im ( aj V< ( (x)) , $ p = aj *| p (x) for a H x G ^ 

g = l,...,£ fe , fc = £ or fc = p, £ = 1,...,^, p=l,...,/c 2 , j = l,...,m; 
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the natural numbers and e p such that 2 £§+ ^ e p = m—2k 1 — k 2 + l with 2k 1 +k 2 ^ r, 



e 5 ^ s f 



i, e 



,ki, E p ^ S p 1, p 



,/c 2 , where fcj is an number of pairs of 



complex-conjugate common eigenvectors of the matrices and fc 2 is an number of real 
common eigenvectors of the matrices B y 

Example 2.18. Consider the sixth-order Lappo-Danilevskii differential system 

dx 
~dt 

R, j = 1,2,3, are linearly independent continuous functions on an interval 
JcR, the commuting matrices 



a 1 (t)A 1 x + a 2 (t)A 2 x + a 3 (t)A 3 x, x£ 



(2.53) 



where aij : J 



3 
-4 
4 
1 

2 



and £? 3 = A 3 



-1 -3 3 5 -2 

3 5-6-5 5 

-3 -5 4 8 -4 

0-1 4 3-3 

-2 -2-1 3 1 

-3 -4 5 6 -2 



3 2 3-3 -3 
2-3-3 2 3 

4 3 5-6 -6 
3 3 4-4 -4 
0-1 0-1 -1 



B 2 =A\- 




-4 
2 
1 

-1 
1 



1 

3 


1 

-1 




6 
-2 
-1 
2 
-1 



2 
-6 

2 

3 
-4 

2 



-2 
4 

-3 

-3 
2 

-2 



(T denotes the matrix transpose). 



Using the eigenvalue A} = 1 + 2i with elementary divisor (A 1 — 1 — 2i) 3 of the matrix B 1 
and the corresponding complex eigenvector v 01 = (1, 0, 1 + i, 1, i, 1), the generalized eigen- 
vector of the 1-st order u 11 = (1,1 + i, 0, 0, i, i), the generalized eigenvector of the 2-nd order 
v 21 = (2 + 2i, 0, 2 + 2i, 0, 2i, 2i), we can construct the scalar functions 

£ 1 . x (gl + X 2)(^l + X 3 + X 4 + gg) + (X 3 + X 5 )(X 2 + X 5 + gg) ^ x ^ x 

11 ' (^1 + ^3 + x 4 + x 6) 2 + ( x 3 + x 5) 2 

vj/ ^ ■ x (x 1 + x 3 + x 4 + x 6 )(x 2 + x 5 + x 6 ) - (x 1 + x 2 )(x 3 + x 5 ) ^ ^ ^ j 



(x x + x 3 + x 4 + x 6 ) 2 + (x 3 + x 5 ) 



*21 : « 



(X 1+ X 3 +X 4 + X 6 )2 + (X 3 + X 5 )2 + -3 + ^4 + X 6 ) 2 + (X 3 + X 5 f) ■ 

■ (( X i+x 3 )( y x 1 +x 3 + x 4 + x 6 ) + (x 3 +x 5 )( Xl +x 3 + x 5 + x 6 )) - ({x 1 +x 2 )(x x +x 3 + x 4 + x 6 ) + 
+ (x 3 + x 5 )(x 2 + x 5 + x 6 )) 2 + ((x x + x 3 + x 4 + x 6 )(x 2 + x 5 + x 6 ) - (x 1 + x 2 )(x 3 + x 5 )) 2 ) , 

*21 : x ->• 7 T T T 2 ^2 T 7 T ^2 f(( X l + X 3 + x 4 + x sY + ( x 3 + ^5) 2 ) ' 

(x 1 + x 3 + x 4 + XqY + (x 3 + x 5 y v v ' 

■ ((Xi + Xg + g^ggXg^gg + gg + gg) - (^l + ^3)( x 3 + x 5)) + (i x l + x 2) ( x l + x 3 + x 4 + x o) + 

+ ( x 3+ x 5)( x 2+ x 5+ x (i)) (i x 3+ x 5)( x l+ x 2) ~ ( x l+ x 3+ X A + x e)( x 2+ x 5+ x g))) for a11 x G 
where X is any domain from the set {x: x 1 + x 3 + x 4 + x 6 / 0} of the space M 6 . 
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Then, the Lappo-Danilevskii differential system (2.53) has the autonomous first integrals 

(X ~~\~ X ^ ~~ \ 

- 4 arctan ^ ^ h 6^ (x) + 2^ (x) ) , 
x 1 + x 3 + x 4 + x 6 J 



- 2 arctan ^ ^ h (x) - ^ (x) I . 

X-i "T" Xo H - X,i -|- Xr- / 



X-^ "T Xg ~t~ ~t~ Xg 

and 

F 3 : x 2* 1 1 1 (x) - 2\i/ 2 1 1 (x) - * 21 (x) for all x G X, 

where the polynomial P: x — >■ (x x + x 3 + x 4 + x 6 ) 2 + (x 3 + x 5 ) 2 for all x 6 M 6 . 

Case 2A. A common complex eigenvector of the matrices Bj, j = 1, . . . ,m, hasn't the 
complex conjugate vector. Autonomous first integrals of the system (2.20) are the functions 



F, : x 



JJ(P € ( X) ) Wo.^*) eX p( - 2( V " WoW*) + 



l 

• (^ 1+ i(^))"°' (2kl+1) exp( - 2^ (2 , i+1) ^ 2fei+1 (x)) fj(^x) 2 ^ exp(2X; h qp ^ qp {x)) 

p=l q=l 

and 

K _ _ 

F 2 : x ^ n( P ?^)) V+/l °' (fcl+ ° ex P ( 2 ( V " Wi+O) ^ + 
£ « 

p=l g=l 

for all x € X, X C DF 1 n DF 2 , 

where the polynomials P^: x — )■ (^°^x) 2 + (z7°£x) 2 for all x £ IR n , the scalar functions 

v ^x * ~ 

: x — >■ arctan for all x G X, £ = 1, . . . , A; 1 , £ = 2^ + 1, the numbers /i^ = h q ^ + h q ^ i, 

v°€x 

h q P = h qp + h qp i, q = 0,...,e k , k = i or k = p, £ = 1, . . . , 2^ + 1, p=l,...,k 2 , are an 
nontrivial solution to the linear homogeneous system 

2fe 1 ^£ fc 2 £ p 

E ( A ? V + £ ^ h tf) + A 2fe 1+ i /l 0,(2fc 1 +l) + J2 ( A p V + E 4p Kp) = °> J = 1, • • • , m. 

5=1 9=1 p=l 9=1 

Here = ^ v^+O =1M, and u°^ 2k ^ = £°>(2*i+i) + ~0,(2fc 1+ i)- are 

complex common eigenvectors of the matrices B- corresponding to the eigenvalues 

\\ = \\ + \\ i, = X 3 V £ = 1, ■ ■ • , k lt and X J 2ki+1 = A J 2fel+1 + X J 2kl+1 i, j = 1, . . . , m; 

i/ 0p are real common eigenvectors of the matrices corresponding to the real eigenvalues 
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Ap, j = 1, . . . , m, p = 1, . . . , k 2 , respectively; the functions \P ^ 
the solution to the functional system (2.49); the numbers 



' all(l £U " 



tit = a i ^^ x ^ Hi = Re ^i vfe = lm ^i & = a i for a11 x G x ' 

q=l,...,£ k , k = £ or fc = p, £ = l,...^, p=l,...,k 2 , j = l,...,m. 

the natural numbers and e p such that 2 ^ &t+ 2^ £ p = m—2k 1 — k 2 with 2/c 1 + l+/c 2 ^ r, 

^ — 1, £ = 1, . . . , k 1 , e p ^ s p — 1, p = 1, . . . , £; 2 , where fc x is an number of complex 
common eigenvectors (this set hasn't complex conjugate vectors) of the matrices Bj, and k 2 
is an number of real common eigenvectors of the matrices B-. 

Example 2.19. Consider the sixth-order Lappo-Danilevskii differential system 

4 



dx 
lit 



(Xj(t)Aj x, x £ 

3=1 



t G J C 



(2.54) 



with linearly independent continuous functions ctj : J 


-> 


IK, J 


= 1,- 


..,4, 


and the matrice 




1 


- 2 


2 





1 




1 









2 








1 


1 







2 


- 2 





- 2 




2 






- 1 


-3 








- 1 


- 1 


A 1 = 





3 


- 2 





- 2 




2 


> ^2 




- 1 


- 3 


- 1 





- 2 


- 2 





- 4 





2 


- 2 




2 




2 


2 





1 





4 




2 


-3 


4 


2 


2 




4 






3 


3 


2 


2 


1 


4 




- 1 


3 


- 2 


- 2 


1 




1 






- 2 


- 1 


- 2 


- 2 


- 1 


-4 




3 











- 1 


- 1 








1 


- 2 


4 





2 


2 




- 1 


2 








1 


1 








- 2 


- 1 


-4 





-4 


-4 


A 3 = 


- 2 
1 


- 1 

2 


2 

- 2 



1 


1 

- 1 


1 

- 1 




; i i ' 1 A 4 = 




-3 
3 


2 

- 4 


- 7 
10 



2 


-5 
7 


-5 
7 




2 


1 


- 1 








- 1 








3 


- 2 


9 





7 


5 




- 1 


- 1 


1 





1 


2 








1 


4 


-5 





-4 


- 2 



The system (2.54) has the eigenvalue \\ = 1 + i with elementary divisor (A 1 — 1 — i) 2 
and the eigenvalue A 2 = 2i with elementary divisor A 1 — 2i. 

Since the eigenvalue \\ = l+i corresponding to the eigenvector v 01 = (1,1 + i, 0, 0, i, i) 

1, ... ,4 (T denotes the matrix transpose) and to the 1-st 



of the matrices B- 



AJ, j 



order generalized eigenvector u 11 = (1 + i, 0, 1 + i, 0, i, i) of the matrix B 1 , we have 
(x 1 + x 2 )(x 1 + x 3 ) + (x 2 + x 5 + x 6 )(x 1 + x 3 + x 5 + x 6 ) 



i\ : x 



P x (x) 

(gi + a 2 )(g 1 + x 3 + x 5 + x 6 ) - (gi + x 3 )(x 2 + x 5 + Zg) 

^ /V,'-) 



for all x 6 X, 



for all x € X, 



where X is a domain from the set {x : x 1 + x 2 / 0, x x + x 3 + x 4 + x 6 / 0} of the space M 6 , 
the polynomial P 1 : x — > (x 1 + x 2 ) 2 + (x 2 + x 5 + x 6 ) 2 for all x 6 M 6 . 

Using the common eigenvector z^ 02 = (1,0, 1 + i, 1, «, 1) of the matrices B 1 , . . . , i? 4 cor- 
responding to the simple eigenvalue A 2 = 2i, we can build autonomous first integrals of the 
Lappo-Danilevskii differential system (2.54) 

F 1 : x^ P 1 (x)(P 2 (x)) 2 exp(-10( / 9 1 (x)+84i 1 1 (x) + 6$ 1 1 1 (x)) for all x G X, 
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F 2 : x -> (P 1 (x)) 3 exp( - 10^ (x) - 4(p 2 (x) + 12 ^(x) + 14 ^(x)) for all x G X, 
where the polynomial P 2 : x — > (x 1 + x 3 + x 4 + x 6 ) 2 + (x 3 + x 5 ) 2 for all x G R 6 , the scalar 
functions ip 1 : x — > arctan %2 — Xs — X6 , <£> 2 : x — > arctan Xs — Xs for all x £ X. 

X-^ ~\~ X 2 X-^ ~t~ X^ ~t~ X^ -\- Xq 

Case 2B. Suppose a function , 7 € {1, ... , A;!}, ft G {1, . . . ,e 7 } hasn't the complex 
conjugate function. Then the system (2.20) has the autonomous first integrals 



1 1 



9=1 

+ 
and 



+ 2^(1 - v**7)((fctf + V(* 1+ o) + ( h i,(k 1+ o - \t) + 

9=1 

2 (fc^ fy 7 *! 7 0r))) n(^) 2/l ° P exp(2^V^ p (x)j for all x G X 



p=l ^ g=l 



F 2 : x n(P f (x)) V +ft o,(, 1+ eX p(2( V " h ,(k 1+ o) M x ) + 

+ 2^(1 - VW((^ + + Utf " + 

+ 2(^ 7 ^ 7 (x) + V^ 7 (x))) n(^x) 2 V ex p(2X;V^(x)) for all x G X, 



where X is a domain from the set DF 1 PlDF 2 , the polynomials P^ : x — > (^°^x) 2 + ( z/°^x) 2 

V°^x 

for all x G W 1 , the scalar functions (f^ : x — > arctan — for all x G X, £ = 1, . . . , k ± , the 

numbers h qi = h ql . + h q ^ i, h qp = h qp + h qp i, q = 0,...,e k , k = £ or k = p, f = 1, . . . , 2k 1: 
p = l,...,k 2 , are an nontrivial solution to the linear homogeneous system 

X ( A | V + S ^ h qt) ~ ^(fc 1+7 ) + X] ( A P V + J] V) = 0, J = 1, . . . , m, 

5=1 g=l p=l g=l 

and 5 is the Kronecker symbol. 

Here z/ ^ = is°£ jy°'( fe i+5) = i/OS, are CO mplex common eigenvectors of the matrices 

Bj corresponding to the eigenvalues A| = a| + A| i, A{ i+ £ = A|, j = 1, . . . , m, £ = 1, . . . , 
respectively; i/ 0p are real common eigenvectors of the matrices corresponding to the eig- 



envalues X P , j = 1, . . . , m, p = 1, . . . , k 2 , respectively; the scalar functions = VI/^ + Vl/^ 
and are the solution to the system (2.49); the numbers 

f^H = S ^ = Re K?' ^ = Im ^' ^ = a i for a11 x G x ' 

g=l,...,e fc , A; = £ or fc = p, £ = 1, . . . , 2k 1: p=l,...,k 2 , j = l,...,m; 
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the natural numbers and e p such that 2 J2 s^+J2 £ p = rn—2k 1 — k 2 +2 with 2k l +k 2 ^ r, 

^ — 1, £ = 1, . . . , k-y, e p ^ s p — 1, p = 1, . . . , fe 2 , where fc x is an number of complex 
common eigenvectors (this set hasn't complex conjugate vectors) of the matrices and k 2 
is an number of real common eigenvectors of the matrices J3-. 

Example 2.20. Consider the sixth-order Lappo-Danilevskii differential system 



dx . . 



3-44102 
-1 3-3 -2 -3 
-3 5 -5 -1 -2 -4 
3 -6 4 4 -1 5 
5 -5 8 3 3 6 
-2 5-4-3 1 -2 



x + a 2 (t) 



0- 4 2 1-1 1 
1 3 1 -1 
6-2-1 2-1 
2-6 2 3-4 2 

1- 6 3 2-2 2 
-2 4 -3 -3 2 -2 



x, xeM 6 , 



where linearly independent functions a 1 : J — > R and a 2 : J — >■ R are continuous. 

Using the complex eigenvalue \\ = 1 + 2i with elementary divisor (A 1 — 1 — 2i) 3 of the 
matrix B 1 = A\ (T denotes the matrix transpose) corresponding to the common eigenvector 
u° — (1,0, 1), to the 1-st order generalized eigenvector v 11 = (1, 1 + i, 0, 0, i, i), and 
to the 2-nd order generalized eigenvector z^ 21 = (2 + 2i, 0, 2 + 2i, 0, 2i, 2i), we can construct 
the autonomous first integrals of this Lappo-Danilevskii differential system 

F x :x^ P(x)exp( - <p(x) - *n(x)), F 2 : x ->■ P(x) exp( - 2<p(x) + 2 ^(x)) , 

P 3 : x ->■ P 2 (x) exp ( - 2(p(x) - ^ 21 (x)), and F 4 :x->- ^ 21 (x) for all x G X, 

where the polynomial P: x — >■ (x : + x 3 + x 4 + x 6 ) 2 + (x 3 + x 5 ) 2 for all x 6 R 6 , the scalar 

x 3 + x 5 



functions <p : x — >■ arctan 



for all x G X, 



*n : x 



X-^ "t" X 3 X 4 -(- Xg 

(xi + x 2 )(x 1 + x 3 + x 4 + x 6 ) + (x 3 + x 5 )(x 2 + x 5 + x 6 ) 
(x x + x 3 + x 4 + x 6 ) 2 + (x 3 + x 5 ) 2 



for all x 6 X, 



^ • X ->■ ( X l + X 3 + x 4 + x 6)( x 2 + x 5 + x e) ~ ( x l + x 2)( x 3 + x 5) £ Qr & jj ^ g j 



* 21 : X 



(x x + x 3 + x 4 + x 6 ) 2 + (x 3 + x 5 ) 2 

1 



, — 7 ttt ( 2((xi + x 3 + x 4 + x 6 ) 2 + (x 3 + x 5 ) 2 ) • 

(x 1 +x 3 + x 4 + x 6 ) 2 + (x 3 + x 5 ) 2 V vv 1 3 6; V3 s>> 

■ (( x i+x 3 )( y x 1 +x 3 + x 4 + x 6 ) + (x 3 + x 5 )(x 1 + x 3 + x 5 + x 6 )) - ((x 1 +x 2 )(x 1 +x 3 + x 4 + x 6 ) + 
+ ( x 3 + x 5)( x 2 + x 5 + x e)) 2 + (i x i + x 3 + x 4 + x 6 )(x 2 + x 5 + x 6 ) - (x x + x 2 )(x 3 + x 5 )) 2 ) , 

~ 2 / 

* 2] : X ->■ rr3 7 rrr I ((Xi + X 3 + X 4 + X fi ) 2 + (x 3 + X 5 ) 2 ) • 

21 (x 1 +x 3 + x 4 + x 6 ) 2 + (x 3 +x 5 ) 2 V vv 1 3 4 b/ V3 b> ' 

■ (( x l+ x 3 + x 4 + x 6)( x l+ x 3 + x 5 + x 6) - ( x l + x 3)( x 3 + x 5)) + (( x l + ^2) ( x l + x 3 + x 4 + x %) + 

+ { x Z+ x b)( x 2+ x b+ x &)){{ x Z+ x b)( x l+ x 2)- ( x l+ x 3+ X A + x e)( x 2+ x 5+ x e))) for a11 x G X ) 
where X is any domain from the set {x: x x + x 3 + x 4 + x 6 / 0} of the space R 6 . 
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2.3.2. Linear nonhomogeneous differential system 

Consider an nonhomogeneous Lappo-Danilevskii differential system 



dx 
~dt 



S ^a j {t)A j x + f(t), teJ, f G C(J), (2.55) 

i=i 

where linearly independent on an interval JcM functions ay. J — > R are continuous, real 
constant n x n matrices A ■ such that A -A, = A,A-, j = 1, . . . , m, k = 1, . . . , m. 

The corresponding homogeneous differential system of system (2.55) is the system (2.20). 

Using Theorems 2.12, 2.13, and 2.14, we can build first integrals of the nonhomogeneous 
Lappo-Danilevskii differential system (2.55). 

Theorem 2.12. Suppose that the conditions of Lemma 2.3 hold. Then the Lappo-Dani- 
levskii differential system (2.55) has the first integral 



(2.56) 



(2.57) 



F : (t, x) -»• vx <p(t) - J vf{r) <p(r) dr for all (t, ijeJxl", 
*o 

where t is a fixed point from the interval J, the exponential function 

t m 

ip: t — ^ expl — / A- 7 <x, (t) dr j for all t G J. 

Proof. By Lemma 2.3, it follows that the Lie derivative of the function (2.56) by virtue 
of the Lappo-Danilevskii differential system (2.55) is equal to 

t 

< BF(t,x)= vxd t ip(t)-d t J vf{T)ip{T)dr + (%vx + f{t)d x vx)ip(t)= for all (i, x) G JxIR n , 

to 

where the linear differential operator %$(t,x) = 2l(i, x) + f(t)d x for all (t, x) G J x IR n is 
induced by the Lappo-Danilevskii differential system (2.55). ■ 

Example 2.21. The second-order Lappo-Danilevskii differential system 
^ = - 3(j + t 2 )x 1 -^x 2 + t (sin t + 24 e* 2 ) e"* 3 , 

■x 1 + 2,( K --t 2 ^ x 2 + t^cosht + 6arctan e - * 3 
such that the coefficient matrix A(t) = 3t 2 A 1 + (j + 6i 2 )^4 2 for all i G J C {t: i / 0}, 



(2.58) 



cte 2 4 
~~dt 



5 4 


and A 2 = 


- 3 - 2 






-8 -7 




4 3 



where the constant matrices A 1 

The matrices 2^ = Aj and S 2 = A 2 , where T denotes the matrix transpose, have 
the eigenvalues A} = 1, A 2 = — 3, and A 2 = — 1, A 2 = 1 corresponding to the linearly 
independent common real eigenvectors u 1 = (2, 1), v 2 = (1, 1). 

Using the scalar functions 

a x : t ->■ 3t 2 for all t G J, a 2 : t ->• i + 6t 2 for all i G J, 
f^.t-t t (sin 4 + 2ie* 2 ) e - ' 3 , / 2 : i — )■ t^cosht + 6arctan e"* 3 for all t G J, 
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ip^t) = exp^ - J (Ajai(t) + \\a 2 {t)) dtj = exp J + 3t 2 ) dt = texpt 3 for all t G J, 
<p 2 (t) = exp^ - J (A^ai(t) + A|a 2 (*)) dt^ = exp ^ ^3t 2 - eft = i expt 3 for all t G J, 
/!(*) = y (2f 1 (t)+f 2 (t))(f 1 (t)dt = J (2t 2 smt+4t 3 e t2 +t 2 cosht+6t 2 arctan ^jdt = -3i 2 + 

2 t 

+ 4t sin t - 2(t 2 - 2) cos i + 2(t 2 - l)e* - 2i cosh i + (i 2 + 2) sinh i + 2t 3 arctan - + 27 ln(t 2 + 9) , 
I 2 (t) = J (f^t) + f 2 (t))<p 2 (t)dt = J (sint + 2te f2 + cosh t + 6 arctan ^ (it = 

2 t 

= - cos i + e* + sinh t + 6t arctan - - 9 ln(i 2 + 9) for all t G J, 

we can build (by Theorem 2.12) first integrals of the Lappo-Danilevskii system (2.58) 
F 1 : {t,x 1 ,x 2 ) -»• te ii {2x 1 + x 2 ) + 3t 2 -4tsint + 2(t 2 - 2) cost - 2(t 2 - l)e* 2 + 

+ 2t cosh t - (t 2 + 2) sinh t - 2t 3 arctan - - 27 ln(i 2 + 9) for all (i , x 1 , x 2 ) G J x R 2 



and 



F 2 : (t , x 1 , x 2 ) — > j e* 3 (x x + x 2 ) + cos t — e f2 — sinh i — 6t arctan ^ + 9 ln(t 2 + 9) 



for all (t,x l ,x 2 ) G J x R 2 . 

The functionally independent first integrals F 1 and F 2 are an integral basis of the Lappo- 
Danilevskii differential system (2.58) on any domain Jxl 2 . 

Theorem 2.13. Suppose that the conditions of Lemma 2.4 are satisfied. Then the Lappo- 
Danilevskii differential system (2.55) has the first integrals 

t 

F p : (t, x) -> 7p (t, s) - J lp (t, f(t)) dt for all (t, x) G J x R n , p = 1, 2, (2.59) 

*o 

where t is a fixed point from the interval J, the scalar functions 

7 X : (i , x) — )■ I fx cos / ^2 ^ a j (t) dr + sin / ^ A J (r) dr ) exp I — / ^ A 7 Oj (r) dr ) , 

7 2 : (t , x) — > [vx cos / A J (r) dr — vx sin / A J (r) dr J exp ( — / A 17 otj (r) dr ) . 

t i =1 ^ i= l t i =1 

Proof. The idea of the proof of Theorem 2.13 is analogous to the proof of Corollary 1.6 
or Theorem 2.2. Formally using Theorem 2.12, we obtain the complex-valued function (2.56) 
is a first integral of system (2.55). Then the real and imaginary parts of this complex- valued 
first integral are the real first integrals (2.59) of the Lappo-Danilevskii system (2.55). ■ 

Example 2.22. Consider the second-order Erugin system of differential equations (right- 
hand side of system satisfy the Cauchy — Riemann equations) [62; 48, pp. 152 - 153] 

doc dx 

= *i + <* 2 (*) *2 + h(t), ~i= ~ Mt) *i + «i (t) x 2 + f 2 (t), (2.60) 

where functions ay. J — > R and /• : J — > R, j = 1,2 are continuous on an interval J C R. 
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Since u l = and v 2 = (1, — i) are common eigenvectors of the matrices B l 



1 








h- 1 



and Br. 



- 1 

1 



corresponding to the eigenvalues A^ = 1, A 2 = —i, and A 2 = 1, A 2 = i, 



respectively, we see that the Erugin system (2.60) has the integral basis (by Theorem 2.13) 

t 

F p : (t,x 1 ,x 2 ) -»• j p (t,x 1 ,x 2 ) - J 'Y ft (T,f 1 (T),f 2 (T))dT for all (t, x 1 ,x 2 ) G J x IR 2 , p=l,2, 

*o 

where t is a fixed point from the interval J, the scalar functions on the domain JxM 2 

t t t 

j 1 : (t, x 1 , x 2 ) — > ^x 1 cos Ja 2 (t) dt — x 2 sin Ja 2 (t)dtj exp^— Ja 1 (t)dtj , 



t t t 

7 2 : (t, x 1 , x 2 ) — > (x 1 sin J a 2 (t) dt + x 2 cos J a 2 (t) dtj exp ^— y a x (t) dtj . 



For example, if a ± : t — , - , a 2 : t -»• 2t, f^.t^t, f 2 :t-> 2t 2 for all 4 G Jc{t: t/0}, 

then the Erugin system (2.60) has the functionally independent first integrals 

t 

F 1 : (t,x 1 ,x 2 ) — > j (x x cost 2 — x 2 sin t 2 ) — cost 2 — J cos t 2 dr for all (t, x l ,x 2 ) £ JxM 2 

*o 

and 

t 

F 2 : (t,x 1 ,x 2 ) — , i (x x sin t 2 + x 2 cost 2 ) — sin t 2 — J sin t 2 dr for all (t, x 1 ,x 2 ) £ J x IR 2 . 

*o 

Note that these first integrals of the Erugin system (2.60) are nonelementary functions. 
Lemma 2.7. Under the conditions of Lemma 2.5, we have 
9 

a 3 v ° x = J20 $v e ~ Px forallxel, j = l,...,m, 9 = 1, . . . , s - 1, (2.61) 

p=0 



where the numbers fiff = A J , fifi = a •^'g(x), j = 1, . . . 



,m, 



Proof. The proof of Lemma 2.7 is by induction on s. 
Let s = 2. Using the functional system (2.31), we get 

v x x = ^(x) z^x for all i£l 

Then, taking into account Lemmas 2.3 and 2.5, we obtain 

z^x = \i^v x x + /Jfis°x for all x£l, j = 1, . . . , m. 

Therefore the system of identities (2.61) for s = 2 is true. 

Suppose that the assertion of Lemma 2.7 is valid for s = e. Then, from the functional 
system (2.31) with s = e + 1 and 9 = e, we have 

e e e—p 

a j v£% = 2 C-i) $ v £ ~ Px + Yj (p-i) Yj ( £ 7) $ v e ~ P ~ Px for all x G X, j = l, . . . ,m. 



p=i 



p=i 



(3=0 
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By the induction hypothesis, so that 

s- - £x = EC=D + E(V) ^ EC^i 1 ) - {£ ~ 5) ^ = 

P=l <5=0 t?=1 

e £—1 e 

= E(?l) /4 C " e_Pii: +E( e « 1 ) ^^"^ = EO fora11 xGX ' i = i.---.^- 

p=l <5=0 p=0 

Consequently the statement (2.61) for s = e+1 is true. Thus by the principle of induction, 
the system of identities (2.61) is true for every natural number s ^ 2. ■ 

Theorem 2.14. Suppose the system (2.55) satisfies the conditions of Lemma 2.5. Then 
the Lappo-Danilevskii system (2.55) on the domain J x X has the first integrals 

e 

F e : (t,x)^ v d x<p(t) - Y, Ke P -^) F P -^ x ) ~ C e( t )^ 9 = 0,. ..,8-1, (2.62) 

p=i 

where X is a domain from the set {x: v°x ^ 0} C R n , the exponential function ip: J — > R 
is given by the formula (2.57), the scalar functions 

1 e- P 



-i : *- >> /((A)w-^( r )+E(!)^( r )<-?( r ))^ p=i,...,9, o = i,. ..,s-i, 

t Q 

C e :t^ ! {y e f{r) V {r) + ^{ 6 p )4{T)C e _ p {r)\dT, 9 = 0,. ..,8-1, for all t € J, 

l o 

m 

4 : 1 E^ Cft j(*) /° r all teJ = aflfe), j = l,...,m), 9 = l,...,s-l, t € J. 

3=1 

Proof. The proof is by induction on s. The case s = 1 was considered in Theorem 2.12. 
Suppose s = 2. Using Theorem 2.12 and Lemma 2.7, we obtain 

*8F 1 (t,x) = i4(t)(v°x<p(t)-C (t)-F (t,x))=0 for all (i,i)eJxX. 

Therefore the scalar function F 1 : J x X — > R is a first integral of system (2.55). 

Suppose that the assertion of Theorem 2.14 is valid for s = e, i.e., the scalar functions 
Fq: T xX — 7- R, 9 = 1, ...,£— 1, are first integrals of the Lappo-Danilevskii system (2.55). 
Then, from Lemma 2.7, we get on the domain J x X 

<BF £ (t, x) = £ (p a4(*) ( (" £ " P ^ " EX-i (*) Vi (*> *) " " - *)) = °- 

p=l ^ n=l ' 

Consequently if s = e + 1, then the scalar function F £ : J x X — > R is a first integral on 
the domain J x X of the Lappo-Danilevskii differential system (2.55). 

Now from the method of building scalar functions (2.62) it follows that the Lappo-Dani- 
levskii differential system (2.55) has the functionally independent first integrals (2.62). ■ 

Example 2.23. The third-order Lappo-Danilevskii differential system 

^ = tx, - (2t + l)x 2 + (t + l)x 3 + \ , ^ = (2t + l) Xl - x 2 + (2t + l)x 3 - 2t 2 , 

(2.63) 

^= -(t + l)x 1 + (2t + l)x 2 -{t + 2)x 3 + t 
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such that the coefficient matrix A(t) = A 1 +tA 2 for all tgR, where the constant matrix 



- 1 

1 - i 
i i 



and A n 



1-2 1 
2 2 
1 2-1 



The matrix B l = A\ , where T denotes the matrix transpose, has triple eigenvalue 
X\ = — 1. The rank of the matrix B 1 — X\E is equal 2. Therefore the eigenvalue A} = — 1 
has x 1 = 3 — 2 = 1 elementary divisor (A 1 + l) 3 . 

The matrices B 1 and B 2 = A 2 have the common real eigenvector v° = (1,0, 1) corres- 
ponding to the eigenvalues X\ = — 1 and A 2 = 0. Also, the matrix B 1 has the 1-st order 
generalized eigenvector u 1 = (0, 1, 0) and the 2-nd order generalized eigenvector v 2 = (0, 2, 2) 
corresponding to the eigenvalue \\ = — 1. 

Now using the common eigenvector v>° and the generalized eigenvectors i/ 1 , u 2 , we can 
build the functions z/°x: x — > x 1 + x 3 , u l x: x — > x 2 , v 2 x: x — > 2(x 2 + x 3 ) for all x G R 3 , 



ty{ : x 



, x 



2(x 1 + x 3 )(x 2 + x 3 )-x 2 l 
(x 1 + x 3 ^ 2 



for all x G X C {x : x l + x 3 / 0}. 



x 1 + x 3 ^ -r 

The real numbers fi\ = a 4 ^(x) = 1, /ii 2 = a± ^(x) = 0, and [i\ = a 2 ^\(x) = 2, 
/x| = a 2 ^>l(x) = 2, where the linear differential operators of first order 

ai ( x ) = (-x 2 + x 3 )d Xi + (x x - x 2 + Xg)^ - (x 4 - x 2 + 2x3)^ for all x G R 3 , 

a 2 (x) = {x l - 2x 2 + x 3 ) d Xi + 2(x x + x 3 ) d x<i - (x 1 - 2x 2 + x 3 ) d Xa for all x G M 3 . 
Using the scalar functions 

a^t-^l, a 2 :t->t, / x : t -)• ^ , f 2 : t - 2t 2 , f 3 : t -> t, <p: t e l for all t G R, 



C„: t-> (*-^)e*, n\:t^2t + l, Kq-. t -¥t(t + 1), C^i-^-^e*, t ^ 2t, 

K$: t -»• t 2 (t 2 + 2t + 2), if 2 : t -»• 2t(t + 1), C 2 : t -»• - (2t 2 - 3t + 4) e* for all t G R, 
we can build (by Theorem 2.14) the basis of first integrals on the space R for system (2.63) 

F :(t,z)-> (x 1 +x 3 -t+^e t , F ± : (t,x) ^ (x 2 + ^)e* - t(t + l)F (t, x), 

F 2 : (t, x) -»• (2x 2 + 2x 3 + 2t 2 - 3t + 4)e* - t 2 {t 2 + 2t + 2)F (t, x) - 2t(t + 1)^(4, x). 

The proof of Theorem 2.14 is true both for real eigenvectors (common and generalized) of 
the matrix B^, and for complex eigenvectors (common and generalized) of the matrix B^. 

In the complex case, from the complex-valued first integrals (2.62) of the Lappo-Danilev- 
skii differential system (2.55), we obtain the real first integrals 

Fg 1 : (t,x) -> ReF d (t,x), F e 2 : (t, x) -> TmF e (t, x) for all (t, x) G J x X, 6 = 0, . . . , 5 - 1. 
Example 2.24. The fourth-order Lappo-Danilevskii differential system 



dxi 1 

—jj- = - (x 1 - x 2 ) - x 4 + t COS t, 



dx 2 
~~dt 



1 

x 1 + — x 2 — x 3 — t sin t, 



dx 3 
~~dt 



= (l + i)x 2 + Ix 3 + x 4 + t 2 , ^=(3-I)x 1 + (2-i)x 3 + ix 4 + 



(2.64) 



d,t 
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has the coefficient matrix of the form A(t) = A 1 + — A 2 for all t G J C {t : t ^ 0}, where 



the constant matrix A 1 





1 
1 

3 



- 1 

1 
1 

2 



and A 



1 


- 1 











1 











1 


1 





- 1 





- 1 


1 



The matrix B 1 = Af , where T denotes the matrix transpose, has two double complex 
eigenvalues \\ = i and A 2 = — i. The rank of the matrix B 1 — X{ E is equal 3. Therefore 
the eigenvalue X\ = i has >c x = 4 — 3 = 1 elementary divisor (A 1 — i) 2 . 

Thus the matrix B l has two elementary divisors (A 1 — i) 2 and (A 1 + i) 2 . 

The matrices B l and B 2 = A 2 have the common complex eigenvector v° = (1, —i, 1,0) 
corresponding to the eigenvalues X\ = i and X\ = 1. Also, the matrix has the generalized 
eigenvector of the 1-st order u 1 = ( — i, 1, 0, 1). 

From the eigenvector and the generalized eigenvector v 1 , we obtain the function 



<S>\ : x 



ix, + x 2 + x 4 for &ll xeXc{x . (xi + X3) 2 + x 2 Q} 



The numbers = i, = 1, = ^^^x) = 1, = o 2 ^'j;(x) = — 1, where 
«i 0*0 = - ^4 ft^ - { Xl + x 3 ) d X2 + (x 2 + x 4 ) ft^ + (3x 1 + 2x 3 ) ft^ for all x G M 4 , 
a 2 (x) = (x x - x 2 ) ft^ + x 2 8^ + (x 2 + x 3 ) d Xs - {x 1 + x 3 - x 4 ) 0^ for all x G M 4 . 

Using the scalar functions 

a ± : t — >■ 1, a 2 : * — ^ J > /1 : * ~ ^ t cost, / 2 : t — ^ — tsint, / 3 : t — ^ t 2 , / 4 : t — > 2t, 
ip: t — > — (cost — isint), C : t — > t + cost + tsint — i(sint — t cost), /4 : t — > 1 — - , 

t 2 1 f cos t 

i^o : t ->■ t - In |t| , : t ->■ - t + — + cos t + 4 sin t + - cos 2t - t cos t - / dt + 

2 2 J t 

( 1 /" sint \ 

+ i ^4 cos t — sin t — - sin 2t + t sm t + / — — at J for all t G J, 

we can build the integral basis of the Lappo-Danilevskii differential system (2.64) 

FQ-.(t,x)— > - (cost (x 1 + x 3 ) — sin tx 2 ) — t — cost — t sint for all (t,x) S J x I 4 , 



F 2 : (t , x) — > — — (cos tx 2 + sin t (x x + x 3 )) + sin t — t cos t for all (t , x) G J 
Fj 1 : (t,x) ->■ - (cost (x 2 + x 4 ) - smtX]) — (t — In ItDFg 1 ^, x) + 



+ t — - — cos t — 4 sm t — — cos 2t + t cos t + 



cost 



(it for all (t, x) G J x 



F x 2 : (t,x) ->■ - (costxi + sint(x 2 + x 4 )) + (t - In |t|) F 2 (t, x) + 



+ 4 cos t — sin t — — sin 2t + t sin t 



sint 



(it for all (t, x) G J x 
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3. Integrals of reducible ordinary differential systems 



3.1. Linear homogeneous differential system 



Consider a system of n ordinary real linear differential equations 

^ = A(t)x, (3.1) 

where x = colon (xi, x n ) from the arithmetical phase space W 1 , the real nxn coefficient 
matrix A: t — > A{t) for all t G J is continuous on an interval Jet. 

One of the most efficient methods for investigation of linear nonautonomous differential 
systems is the method of reducibility. In this method linear nonautonomous systems are redu- 
ced to linear systems with constant coefficients by various transformation groups [63; 8; 10]. 
The idea of this method is due to A.M. Lyapunov [57]. He studied linear periodic systems 
and showed that there exists a transformation which does not change the character of the 
growth of solutions and reduces a system with periodic coefficients to a system with constant 
coefficients. Lyapunov called systems having this condition reducible systems. The general 
theory of reducible systems was developed by N.P.Erugin in his article [58]. 

Later on, the notion of reducible system was given for systems of difference equations [64] 
and for multidimensional differential systems (see [65 - 69; 9, pp. 72 - 82, 242 - 246]). 

Let G be a multiplicative group of real continuously differentiable on the interval J 
nonsingular matrices of order n. The differential system (3.1) is called reducible with respect 
to the nonautonomous transformation group G if there exist a constant matrix B and a 
matrix g G G such that the linear transformation y = g(i) x reduces the nonautonomous 
differential system (3.1) to the system with constant coefficients 



By Theorems 3.1 - 3.4, using eigenvectors and corresponding eigenvalues of the matrix 
C = B T , where T denotes the matrix transpose, and a transformation matrix g £ G, we 
can construct first integrals of the reducible system (3.1). The following basic statements 
(Lemmas 3.1 and 3.2) are base for the method of building integral basis of system (3.1). 

Lemma 3.1. Suppose the system (3.1) is reducible to the system (3.2) by a transformation 
matrix gGG, and v is a real eigenvector of the matrix C corresponding to the eigenvalue A. 
Then the linear function p : (t, x) — > v g(i) x for all (t, x) G J x W 1 is a partial integral of 
the reducible differential system (3.1) such that 



where the linear differential operator 2l(t,x) = d t + A(t)xd x for all (t,x) G J x M n . 
Indeed, using the matrix identity (3.3), we get 

2lp(i,x) = d t p(t,x) + A(t)xd x p(t,x) = vg'(t)x + A{t)xvg{t) = 

= u(B g(t) - g(t) A(t)) x + A(t) x v g(t) = vB g(t) x + (A(t)x vg(t) - ug(t) A(t)x) = 

= Cvg{t)x = \ i>g{t)x = Xp(t,x) for all (t,x) G J x W 1 . ■ 

Lemma 3.2. Suppose the system (3.1) is reducible to the system with constant coefficients 
(3.2) by a transformation matrix g G G, and v = v + vi (Rev = u, Imu = v) is a complex 
eigenvector of the matrix C corresponding to the eigenvalue A = \+\i (Re A = A, ImA = A). 




(3.2) 



(3.3) 



2lp(t, x) = Xp(t, x) for all (t, x) G J x W 



n 
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Then the Lie derivatives of the scalar functions 

P: (t,x) -»• {tg{t)x) 2 + (ug(t)x) 2 for all (t,x) £ JxM" 

and 

(p: (t,x) -» arctan ; &v ; for all (t,x) e A, A C { (t, x) : t G J, 1/ g(t) x / 0} C J x i n , 
fg(i)x 

6y virtue of the reducible system (3.1) are eouai to 

2lP(t,x) = 2AP(t,x) /orsII(J,i)eJxR n and 2l^(t,x) = A /or all (t, x) G A. 

Proof. Formally using Lemma 3.1, we obtain the linear function p: (t,x) — > vg(t)x for 
all (t,x) £ JxP is a complex-valued partial integral of the reducible linear differential 
system (3.1) and the following identity holds 

2L(vg(t)x + iVg(t)x) = (A + Xi){vg{t)x + ivg{t)x) for all (t,x)GJxR n . 
This complex identity is equivalent to the real system of identities 

%Lvg(t)x = \vg(t)x- Xvg{t)x for all (t,x) £ J xW\ 

'&vg(t)x=\vg{t)x + \vg{t)x for all (i,x)GJxR n . 
Using this system of identities, we have 

%P(t,x) = 2l((^g(t)x) 2 + (ug(t)x) 2 ) = 2^g(t)x2l^g(t)x + 2vg(t)x2LVg(t)x = 

= 2tg(t)x(\vg{t)x-\Vg{t)x) + 2Vg{t)x(\vg{t)x + \tg{t)x) = 

= 2*\((vg(t)x) 2 + (ug(t)x) 2 ) = 2\P(t,x) for all (i,i)eJxl n ; 

= 21 arctan Mil = . ^MlilM^ = 

ug(t)x 1 , ("g(*)*r (^g(t)x) 2 



(i^g(t)x) 2 

^g(t)x (A^g(t)x + A^g(t)x) - z?g(i)x (\i> g(t)x -Xvg{t)x) 



A for all (t, x) G A. 



(i/g(t)x) 2 + (i/g(t)x) 2 

Now we can state the Theorems 3.1 - 3.4 for building first integrals of system (3.1). 
Theorem 3.1. Let the conditions of Lemma 3.1 be satisfied. Then the scalar function 

F:(t,x)-> i/g(t)xexp(-Ai) for all (£,x)gJxR™ 
is a _/irsi integral on the domain J x R ra of the reducible differential system (3.1). 
Proof. From Lemma 3.1, we get 

QLF(t,x) = 21 (vg(t)x exp( - Xt)) = exp( - Xt) QLvg(t)x + z/g(f)x2lexp( - Xt) = 

= Xvg(t)x exp( - At) + ug(t)xd t exp( - Xt) = for all (t,x)GJxR n . 

Therefore the function F : J x R n — > R is a first integral of the reducible system (3.1). ■ 
Example 3.1. The second-order nonautonomous differential system [8, p. 155] 

^± = (t 2 +t + 2)x 1 -(t 3 + t 2 + t-l)x 2 , ^ = (t + l)x 1 -(t 2 + t-l)x 2 (3.4) 
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is reducible to the system with constant coefficients — ^ = y 1 , —j^ = 2x 2 by the polynomial 

-t 1 + 1 2 

group 1 P(2) with the transformation matrix g(t) = for all tel. 

Using the real eigenvectors v 1 = (1, 0), v 2 = (0, 1) of the matrix C = diag(l, 2) and the 
corresponding eigenvalues X 1 = 1, A 2 = 2, we can build (by Theorem 3.1) the basis of first 
integrals on space R 3 for the reducible differential system (3.4) 

F 1 : (t, Xl ,x 2 ) ( - tx x + (1 + t 2 )x 2 )e"*, F 2 : (t, x x , x 2 ) -> (x 1 - tx 2 ) e~ 2t . 

Theorem 3.2. Let the conditions of Lemma 3.2 be satisfied. Then the reducible differential 
system (3.1) has the first integrals 

F 1 : (t,x) -> ((^g(t)x) 2 + (ug(t)xf) exp( - 2 At) for all (t,x)eJxR n 

and 

F 2 : (t,x) ->■ arctan " g ^ X _ At for all (t, x) G A, A C {(t, x) : t G J, £ g(t)x / 0}. 
fg(t)x 

Proof. It follows from Lemma 3.2 that 
Sli^x) = exp(-2At)2l((^g(t)x) 2 + (?g(t)x) 2 ) + ((£ g(t) x) 2 + (?g(t) x) 2 ) 2lexp(-2 A t) = 

= 2Ai ? i(t,x) + ((z^g(t)x) 2 + (?g(t)x) 2 ) 9 4 exp( - 2 At) = for all (t,x)eJxR Tt , 

2lF 2 (t,x) =St arctan — gt(At) = A-d t (At)=0 for all (t,i)eA. 

i/g(t)x 

Therefore the scalar functions F 1 : J X R n — >■ R and i 7 ^ : A — >■ R are first integrals of the 
reducible linear homogeneous differential system (3.1). I 

Example 3.2. The linear differential system [70, pp. 125 - 126] 



~r^~ = — x 2 + a/2 cos2tx 3 + a/2 sin2tx 4 , = x 1 + y/2 sin 2tx 3 — a/2 cos 2tx 4 , 



(3.5) 



dx Q . _ /- dx 



-7 1 = - a/2 sin2tx 1 + \/2 cos2tx 2 - x 4 , — - = a/2 cos 2tx x + a/2 sin2tx 2 + 



dt 1 ' v z ~ 4 ' dt 1 ' z ' Xs 

is reducible by the group 2 P 4 (27r). Indeed, the 27r-periodic nondegenerate transformation 
y l = costx 1 +sintx 2 , y 2 = — sintx 1 +costx 2 , y 3 = costx 3 +sintx 4 , y A = — sintx 3 +costx 4 
reduces the system (3.5) to the linear system with constant coefficients 

^1-J2 V ^1--J2 V ^3-a/2w ^i-A/2y 

dt~ V2y ^ dt ~ V2y *> dt~ V2V2 > dt~ V2yi - 

By Theorem 3.2, using the eigenvalues A 1 = 1 + i, A 2 = — 1 + i and the corresponding 
complex eigenvectors v 1 = (1 + — i, a/2, y/2i), v 2 = (1 +i, — y/2i, — a/2) of the 

matrix C, we can build the basis of first integrals for the reducible system (3.5) 



1 P(n) is the multiplicative group of n x n polynomial matrices with nonzero constant determinants. 

2 P n (uj) is the multiplicative group of c^-periodic invertible continuously differentiable square matrices of 
order n. At the same time P n {uS) is a sub-group of the Lyapunov group L(n). 

L(n) is the multiplicative group of invertible continuously differentiable on T — (0; + 00) square matrices 
of order n such that these matrices and their inverse matrices are bounded on the interval T. 
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F 1 : (t, x) — > ^((cos t — smt)x l + (cost + sint)x 2 + a/2 costx 3 + y/2 sintx 4 ) 2 + 
+ ((cost + s'mt)x 1 + (sint — cost)x 2 — a/2 sintx 3 + a/2 cos t x 4 ) 2 J e~ 2t for all (t, x) G 



F 2 : (t, x) — > arctan 



(cos t + sin t) x 1 + (sin t — cos t) x 2 — \/2 sin t x 3 + a/2 cos t x 4 
(cos t — sin t) x 1 + (cos t + sin t) x 2 + a/2 cos t x 3 + a/2 sin t 



X , 



for all (t, x) G A, 

F 3 : (t, x) — > ^((cos t + sint) x 1 + (sint — cost) x 2 + \/2 sintx 3 — \/2 cos t x 4 ) 2 + 
+ ((cost — sint) x 4 + (cost + sint) x 2 — a/2 cost x 3 — \/2 sint x 4 ) 2 ^e 2 ' for all (t,x) G 



F 4 : (t,x) — )■ arctan 



(cos t — sin t) x 1 + (cos t + sin t ) x 2 — \/2 cos t x 3 — y/2 sin t x 4 



(cos t + sin t) x 1 + (sin t — cos t ) x 2 + a/2 sin t x 3 — a/2 cos t x 4 
for all (t, x) G A, 

where a domain A C {(t, x) : (cos t — sint) x 1 + (cos t + sint)x 2 + a/2 cos tx 3 + a/2 sintx 4 7^ 0, 
(cost + sint) x 4 + (sint — cost) x 2 + a/2 sint x 3 — a/2 cost x 4 / 0} of the space M 5 . 

Example 3.3. The third-order Euler differential system 



dx 1 



dx<, 



dxi 



1 



dt ^ dt 

is reducible to the linear system with constant coefficients 
dyi 



dt "t3 Xl t2 X2 



• X, 



(3.6) 



dr 



dy 2 , dy 3 

V2, ~^: = y2 + y^ ~dr =Vl ~ y2 



by the exponential group 1 Exp(3) with the transformation 







Vi = 


x \1 










1 


The matrix C = 


1 


1 


- 1 







1 






- f2. 



'-.3- 



t = e T . 



has the eigenvalues X 1 = 1, X 2 = i, and A 3 = — i corres- 



ponding to the eigenvectors v x = (1, 1, 1), v 2 = (1, — and z/ 3 = (1, — 1, — i). 

An integral basis of the reducible system (3.6) is the scalar functions 



F 1 : (t, x) — > - x 1 + x 2 + tx 3 for all (t, x) G (0; + 00) x I 
F 2 : (t, x) ->• (x 4 - tx 2 ) 2 + t 4 x| for all (t, x) G (0; + 00) x 



(by Theorem 3.1), 



(by Theorem 3.2), 



and 



F 3 : (t, x) — > arctan 



t 2 Xn 



- lnt for all (t,x) G A (by Theorem 3.2), 



where A is any domain from the set {(t, x) : t > 0, x l — tx 2 / 0}. 



Exp(n) is the multiplicative group of invertible continuously differentiable on the interval (0; + 00) 
square matrices g of order n such that lim — ||g ±1 (t)|| = 0. 

t— > + oo t 

At the same time the Lyapunov group L(n) is a sub-group of the exponential group Exp(n). 
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Theorem 3.3. Suppose the system (3.1) is reducible to the system (3.2) by a transforma- 
tion matrix g G G, and A is an eigenvalue with elementary divisor of multiplicity m ^ 2 
of the matrix C corresponding to a real eigenvector v° and to a real generalized eigenvector 
u 1 of the 1-st order. Then the reducible system (3.1) has the first integral 

v l g(i) x 



F: (t,x) 



t for all (t, x) G A, A C {(t, x) : t G J, i/°g(t) x / 0}. (3.7) 



i/°g(i)x 

Proof. From the matrix identity (3.3) it follows that 

2lz/g(t)x = fyi/ 1 g(t)x + A(t)a;a i .i/ 1 g(t)a: = z^g'(i)x + A^xz^g^) = 
= v l (B g(t) - g(t) A(t)) x + A(t) x i/ 1 g(t) = g(t) x = CV 1 g(t) x = (Ai/ 1 + i/°)g(i) x 

= A zv 1 g(t)x + i/°g(t)x for all (t, x) G J x R ra . 
Using this identity and Lemma 3.1, we get on the domain A 
v 1 g(t)x 



2lF(i,x) = 21 



v°g{t)x 



(A z/g(t)x + i/°g(t)x) v°g(t)x - A v°g(t)x v l g(t)x _ 
(^g(i)x) 2 



0.1 



In the complex case, if A is a complex eigenvalue, then from the complex-valued first 
integral (3.7) of the reducible differential system (3.1), we obtain the real first integrals 



and 



F x : (i,x) 



F 2 : (t, x) 



u°g(t)x v l g(t)x + z^°g(t)x u 1 g(t)x 



{^g{t)xf + (v°g(t)x) 2 

u°g(t)x V l g{t)x — V°g(t)x t l g(t)x 
(>g(t)x) 2 + pg(i)x) 2 



- t for all (t, x) G A 



for all (t, x) G A 



where A is any domain from the set {(i,x): t G J, (^°g(i)x) 2 + (?°g(i)x) 2 / 0} of the 

space R™ +1 , the real vectors u k = Keb ,k , u k = lmu k , k = 0, 1. 

Example 3.4. Linear Hamiltonian systems of second order reducible by orthogonal trans- 
formation group 1 are the linear differential systems of the form [8, pp. 142 - 143] 



dx 
~dt 



IA(t)x, xG 



A: t 



a(t) m - i3{t) 



for all t G J, 



(3.8) 



where the simplectic matrix I 



1 

1 



, the functions a: t — > acos2(p(t) — bsin2ip(t), 



/3: t — > asin2tp(t) + 6 cos 2ip(t), Y>: t — > iff(t) + c for all t G J are continuously differentiable 
on the interval JcR, and a, 6, c are some real numbers. 

The orthogonal transformation y x = cosip(t)x 1 — sin <£>(£) x 2 , y 2 = sinip(t)x 1 + cos (p(t) x 2 
reduces the linear differential system (3.8) to the linear autonomous Hamiltonian system 
dy 



—— = IBy with the constant matrix B 
dt 



c + b a 
a c — b 

Consider the real number D = a? + b 2 — c 2 . We have three possible cases for building 
first integrals of the reducible linear Hamiltonian system (3.8). 



1 O(n) is the multiplicative group of continuously differentiable on an interval T C K squarte orthogonal 
matrices of order n. At the same time if T = (0; + oo), then the orthogonal group of transformations 0(n) 
is a sub-group of the Lyapunov group L(n). 
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Let D = 0. Then, using the real eigenvector v° = (c+b, a), the real 1-st order generalized 
eigenvector v l = (c + b, a — 1), and the corresponding double eigenvalue A 1 = 0, we can 
build (by Theorems 3.1 and 3.3) the integral basis on a domain A of system (3.8) 

F-y : (t, x) — > ((c + b) cos ip(t) + a sin <p(t)) x x + (a cos ip(t) — (c + b) sin (fi(t)) x 2 , 

{(c+b) cos ip(t) + (a— 1) sin ip(t))x 1 + ((a — 1) cos tp(t) — (c+b) sin <p(t))x 2 
((c + 6) cos ip(t) + a sin (p(t)) x x + (a cos ip(t) — (c + 6) sin ip(t)) x 2 

where A is any domain from the set 

M = [(t, x) : t G J, ((c + b) cos ip(t) + a sin <£>(£)) + (a cos </?(£) — (c + 6) sin<£>(i))x 2 7^ 0}. 

Let D > 0. Then, using the real eigenvectors v k = (c + b,a — X k ), k = 1,2, and the 
corresponding real eigenvalues A x = y/D, A 2 = — y/D, we can construct (by Theorem 3.1) 
the basis of first integrals for the Hamiltonian system (3.8) 

F k : (t, x) -+ ({(c+b) cos + (a- A fc ) sin <£>(£)) x : + ((a- A fe ) cos ^(t)-(c+6) sin <p(i))x 2 ) e~ Afe * 

for all (£,x) G Jxl 2 , fc = 1, 2. 

Let D < 0. Then, using the eigenvectors z^ fe = (c + 6, a — A fe ), A; = 1, 2, and the corres- 
ponding complex eigenvalues A x = y 7 — D i, A 2 = — \/ — D i, we can find (by Theorem 3.2) 
the functionally independent first integrals of the Hamiltonian system (3.8) 

Fy : (i, x) — )■ (((c + 6) cos ip(t) + a sin <£>(£)) x x + (a cos <p(t) — (c + b) sin </?(£)) x 2 ) 2 — 

- D (sirup (t)x x + cos (p(t)x 2 ) 2 for all (£,x)gJxR 2 , 



P /. w . (sin y (t) X! + cos <g (t) x 2 ) ^— - 

F 2 : (i, x) -)• arctan 7- — - — — - 1 — — h V - D t 

((c + b) cos ip(t)+asm ip(t))x 1 + (a cos <p(t) — (c + b) sin </?(£)) x 2 

for all (i, x) G A C M, M C R 3 . 

Theorem 3.4. Suppose the system (3.1) is reducible to the system (3.2) fry a transforma- 
tion matrix g G G, anc? A is t/ie eigenvalue with elementary divisor of multiplicity m ^ 2 
0/ t/ie matrix C corresponding to a real eigenvector u° and to real generalized eigenvectors 
u k , k = 1, . . . , m — 1. T/ien £/ie system (3.1) /jas i/ie functionally independent first integrals 

F^: (t,x) -+ * { (t,s) for all (t,x) e A, £ = 2, ...,m-l, (3.9) 

where the functions ^ : A — >■ R, £ = 2, . . . , m — 1, are £/te solution to the functional system 
k 

u k g(t)x = J2 (rZ\) * T (*. *) ^~ T gW* /or all (t,x)eA, k = l,...,m- 1, (3.10) 

r=l 

and A is any domain from the set {(t,x): t G J, z^°g(i)x / 0} C J x R n . 

Proof. From the matrix identity (3.3) and the notion of generalized eigenvectors for matrix 
(Definition 1.1) it follows that 

Vlv k g(t) x = d t i> k g(t) x + A(t) x d x i> k g(t) x = v k g'(t) x + A(t) x i> k g(t) = 

= v k (B g(t) - g(t) A(t)) x + A(t) x v k g(t) = v k B g(t) x = Cv k g(t) x = (\v k + ku k - 1 )g(t) x = 

= A v k g(t) x + kv k_1 g(t) x for all (t, x) G J x R n , k = 1, . . . , m - 1. 
Hence using Lemma 3.1, we obtain the system of identities 
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(3-11) 



X.v°g(t)x = \i>°g(t)x for all (t, x) G J x W n , 

$lv k g(t)x = \v k g(t)x + kv k ~ l g(t)x for all (t, x) G J x R n , k = 1, . . . , m — 1. 

The functional system (3.10) has the determinant (v°g(t) x)" 1 ^ 1 / for all (t, x) G A, 
where a domain A C {(t,x): t G J, v°g(t) x 7^ 0}. Therefore there exists the solution 
* r , r = 1, . . . , m — 1, on the domain A of the functional system (3.10). Let us show that 



m k (t,x) 



1 for all (t,x) G A, k = 1, 

for all (t, x) G A, fe = 2, ...,m — 1, 



(3.12) 



The proof of identities (3.12) is by induction on m. 

For m = 2 and m = 3, the assertions (3.12) follows from the identities (3.11). 
Assume that the identities (3.12) for m = e is true. Then, using the system of identities 
(3.11) and the identities (3.10) for m = e + 1, m = e, we get 



e-l 



T=l 



T=l 



+ i^- 1 g(*)a; + i/°g(t)x 2l* £ (t,x) for all (i,x) G A. 

Now taking into account the system (3.10) with k = £ — 1 and k = e, the identity (3.11) 
with k = £, and z^°g(t) s/fl for all (i, x) G A, we have 

a* e (t,s)=0 for all (t, x) G A. 

This implies that the identities (3.12) for m = e + 1 are true. So by the principle of 
mathematical induction, the statement (3.12) is true for every natural number m ^ 2. 

Now from the method of building scalar functions (3.9) it follows that the reducible dif- 
ferential system (3.1) has the functionally independent first integrals (3.9). I 

The proof of Theorem 3.4 is true both for the case of the real eigenvalue A and for the 
case of the complex eigenvalue A (ImA / 0). In the complex case, from the complex- valued 
first integrals (3.9) of system (3.1), we obtain the real first integrals of system (3.1) 

F£: (t,x) -> Re* € (t,s), F| : (t,x) -»• ImVfax) for all (t, x) G A, f = 2,...,m-l, 

where A is any domain from the set {(£, x): t G J, (^°g(t)x) 2 + (z?°g(t)x) 2 7^ 0}. 
Example 3.5. The third-order Euler differential system 



dxi 



x 2i 



d,Xn 



dXr, 



1 



1 



dt ~ Z1 dt 3 ' dt 

is reducible to the linear differential system with constant coefficients 



(3.13) 



dVi 
dr 



d V2 

V2, -^r = V2 + y^ 



dy-y 

-jf = yi-V2 + 2 % 



by the exponential group Exp(3) with the transformation 



The matrix C 



Vi ~ x ii V2 ~ t- x 2' % — ^ x 3i t — e T . 

has the triple eigenvalue X 1 = 1 with the elementary divisor 



1 

1 1-1 
1 2 

(A — l) 3 corresponding to the eigenvector u°= (1,-1,1), to the generalized eigenvector of the 
1-st order u 1 = (0, 0, 1), and to the generalized eigenvector of the 2-nd order v 2 = (0, 2, 0). 
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A basis of first integrals for the reducible system (3.13) is the scalar functions 
F 1 : (t, x) ->■ j x 1 - x 2 + tx 3 for all (t, x) G (0; + oo) x IR 3 (by Theorem 3.1), 
t 2 x 

F 2 : (t, x) -»• ^ — Jn/ for all (t, x) G (0; +oo)xM 3 (by Theorem 3.3), 

and (by Theorem 3.4) 

F 3 :(t,x)^ 2 ^2(^i-^2 + ^3)- * 4 s§ forall (t a;)€A 
d v 7 (x x - tx 2 + t 2 x 3 ) 2 v 7 

where A is a domain from the set {(t, x) : t > 0, x 1 — tx 2 + t 2 x 3 / 0}. 

3.2. Linear nonhomogeneous differential system 

A real linear nonhomogeneous differential system of the n-th order 

- = A(t)x + f(t), feC(J), (3.14) 

such that the system (3.1) is the corresponding homogeneous system of system (3.14). 

The system (3.14) is called reducible with respect to the nonautonomous transformation 
group G if there exist a constant matrix B and a matrix g G G such that the linear 
transformation y = g(i)x reduces the system (3.14) to the system 

^ = By + g(t)f(t), y G M. n . (3.15) 

Let us remark that if the nonhomogeneous system (3.14) is reduced to the system (3.15), 
then homogeneous system (3.1) is reduced to the system with constant coefficients (3.2). 
The reducible system (3.14) is induced the linear differential operator of first order 

»(t, x) = d t + (A(t) x + f(t)) d x = 2l(t, x) + f{t) d x for all (t, x) G J x R n . 

3.2.1. Case of simple elementary divisors 

If the matrix C is diagonalizable, then using Theorem 3.5 and Corollary 3.1, we can build 
first integrals of the reducible nonhomogeneous differential system (3.14). 

Theorem 3.5. Let the system (3.14) be reducible to the system (3.15) by a transformation 
matrix g G G, let u be a real eigenvector of the matrix C corresponding to the eigenvalue A. 
Then a first integral of the reducible system (3.14) is the scalar function 

t 

F: (t,x) -> vg(t)x exp( - At) - J z/g(C)/(C)exp( - \()d( for all (t,x) G J x R n , (3.16) 

*o 

where t Q is a fixed point from the interval J. 

Proof. From Lemma 3.1 it follows that on the domain J x R n : 

t 

«8F(t,x) = 5S^g(t)x exp( - Xt) - 53 J i/g(C)/(C) exp( - AC) = - A^g(t)xexp( - At) + 

t 

+ %vg{t)x exp( - At) + f(t) d x vg{t)x exp( - At) - d t J ^g(C)/(C) exp( - AC) d( = 0. 

Therefore the function (3.16) is a first integral of the reducible system (3.14). ■ 
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Corollary 3.1. Let the system (3.14) be reducible to the system (3.15) by a transformation 
matrix g G G, and let v = v + v i (Re v = v, liau = u) be an eigenvector of the matrix C 

corresponding to the complex eigenvalue A = \ + \i (Re A = A, ImA = A). Then first integrals 
of the reducible system (3.14) are the scalar functions 

t 

F g : (t,x) a e (t,x) - J a g ((, f(()) d( for all (t, x) £ J x W 1 , 9 = 1,2, (3.17) 

*o 

where t is a fixed point from the interval J, the functions 

a x : (t, x) — > (vg(t)x cos At + vg(t)x sin At) exp( — \t) for all (t,x) G J xl" 

a 2 : (t,x) — > (ug(t)x cos At - vg(t)x sin At) exp( - At) for all (t,x)eJxR n . 

Proof. Formally using Theorem 3.5, we obtain the complex- valued function (3.16) is a 
first integral of system (3.14). Then the real and imaginary parts of this complex- valued first 
integral are the real first integrals (3.17) of the reducible system (3.14). ■ 

Example 3.6. Consider the second-order linear nonhomogeneous differential system 1 

= (a sin ujt cos uit) x l — (b + a cos 2 uit) x 2 + /i (t) , 
1 (3.18) 
= (b + a sin 2 uit ) x 1 — (a sin uit cos oot ) x 2 + / 2 (t) , 

where a, b, and are real numbers such that uj / b and w / a + b, the scalar functions 
f 1 : J — > M. and f 2 : J — > M are continuous on an interval JcK. 

The system (3.18) is reducible to the nonhomogeneous system with constant coefficients 

dy dy 

= (uj-a-b)y 2 + f^t) cos ujt + f 2 (t) sin ujt, = (b-u)y 1 - f 1 (t)sinujt + f 2 (t) cos ujt 

by the orthogonal group of transformations 0(J, 2) with the transformation 

y 1 = cos ujtx 1 + smut x 2 , y 2 = — sinujtx 1 + cosujtx 2 

Consider the real number D = (b — uj)(uj — a — b). We have two possible cases for building 
first integrals of the reducible linear differential system (3.18). 

Let D > 0. Then, using the real eigenvectors v k = (b — uj, X k ), k = 1,2, of the matrix 



C = 



b-uj 
uj — a — b 



and the corresponding real eigenvalues A 1 = VD, X 2 = — VT), we can build (by Theorem 3.5) 
the basis of first integrals for the reducible system (3.18) 

F k : (t, x 1 ,x 2 ) — > (((6 — uj) cos uit — X k sinujt)x l + (X k cos ujt + (b — uj) sinujt)x 2 )e~ Xkt — 

t 

- J ( ((b - uj) cos uj( - X k sin cu() f x (C) + (X k cos ui( + (b - uj) sin ui() f 2 (C)) e~ x ^ d( 
*o 

for all (t JeJxR 2 , k = l,2. 



1 Note that the linear homogeneous differential system corresponding to the linear nonhomogeneous differ- 
ential system (3.18) is a simplified differential system (has not nutations), which is describing free oscillations 
of gyroscopic pendulum in twin gyrocompass [17, pp. 528 - 529]. 
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Let D < 0. Then, using the eigenvectors u k = (oj — 6, A 3 _ fc ), k = 1,2, and the cor- 
responding complex eigenvalues A 1 = \J — Di, A 2 = — \J — Di, we can (by Corollary 3.1) 
construct the functionally independent first integrals of the reducible system (3.18) 

t 

F g : (t, Xl ,x 2 ) -> a e (t, Xl ,x 2 ) - J a g (C, /^C), f 2 (0) d C for a11 (t,x 1 ,x 2 ) GJxR 2 , 9 = 1,2, 

to 

where t is a fixed point from the interval J, the scalar functions 



a 1 : (t, x l ,x 2 ) — > ( (oj — b) cos y/ — D t cos ut + V — D sin V — -D £ sinwt) x x + 
+ ( (uj — b) cos \J — D t sin u;t — \/ — D sin \/ — D t cos wt) x 2 for all (t , x 1 , x 2 ) G J x 



a 2 : (*) i x 2 ) — > ( V — D cos V — D t sin wt — (cj — 6) sin V — D t cos ut) x 1 — 
— (V — D cos V — D t cos wt + (to — b) sin V — D t sin wt) x 2 for all (t, x l , x 2 ) G J x M 2 . 

3.2.2. Case of multiple elementary divisors 

If the matrix C has multiple elementary divisors, then using Theorem 3.6 and Corolla- 
ry 3.2, we can build first integrals of the reducible nonhomogeneous differential system (3.14). 

Theorem 3.6. Let the system (3.14) be reducible to the system (3.15) by a transformation 
matrix g G G, and let A be the eigenvalue with elementary divisor of multiplicity m ^ 2 
of the matrix C corresponding to a real eigenvector u° and to real generalized eigenvectors 
v k , k = 1, . . . , m — 1. Then the system (3.14) has the functionally independent first integrals 

k-l 

F k+1 : (t,x) -> z/g(t)xexp(- At) - t k ~ T F T+1 (t, x) - C k (t) 

T=0 

(3.19) 

for all (t,x) <E J x R n , k = 1, . . . , m - 1, 

where the first integral (by Theorem 3.5) 

F : : (t,x) i/°g(t)x exp( - At) - C (t) /or a// (t,x) G J x R n , 
the scalar functions 
t 

C k :t^ y , (^g(C)/(C)exp(-AC) + A:C fc _ 1 (C)) ( iC for all t e J, k = 0,...,m-l, 

*o 

and t is a fixed point from the interval Jcl. 

Proof. The proof of Theorem 3.6 is by induction on to. 
By the system of identities (3.11), it follows that 

3S(z/ e g(t)x exp(-At)) = -\u e g{t)x exp( - At) + (2U e g(t) x + f(t) d x v e g{t) x) exp(-At) = 

= - Az; £ g(t)x exp(- At) + (\v £ g{t) x + e ^ £_1 g(t) x) exp( - At) + v £ g(t)f(t) exp( - At) = 

= (ei/ e-1 g(t)a; + i/ e g(t)/(t))exp(- At) for all (t, x) G J x M n , e=l...,m-l. 

Therefore, we have 

<8(i/ £ g(t)x exp( - At)) = (eu e ~ 1 g(t)x + v e g{t)f{t)) exp( - At) 

(3.20) 

for all (t, x) G J x W l , e = 1 . . . , m - 1. 
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Let m = 2. Using the identities (3.20) with e = 1, we get 

53 F 2 (t , x) = <B (^g(t)x exp( - At) - t F x (t, x) - C 1 (t)) = 

= (u°g(t)x + ^g(t)/(t)) exp( - At) - Fjfrx) - (^g(t)/(t) exp( - At) + C (t)) = 

= (z/°g(t)x exp( - At) - C (t)) - F x (t, x) = for all (t, x) G J x R n . 

Hence the function F 2 : J X R n — )• R is a first integral of the reducible system (3.14). 
Assume that the functions (3.19) for m = ii are first integrals of system (3.14). Then, 
from the identities (3.20) for the function F +1 it follows that 

5SF M+1 (t,x) = <8(>g(t)x exp(-At) - £ (?) t^F T+1 (t, x) - C^tj) = 



T=0 

/i-2 



= (/./-^Wx + ZgW/W) exp(- At) - M J^Cr 1 )*" T ^r+lft*) - 



T = 

- ^(t,x) - (^g(t)/(t) exp(-At) + liC^t)) = 
= ^(^- X g(t)x - J] C; 1 )t"" T " 1 J F T+1 (t, x) - C^t)) - >iF^t, x) = for all (t, x) G JxR". 

T = 

This implies that the scalar function F +1 : JxR" -> R (for m = //+ 1) is a first integral 
of the reducible linear nonhomogeneous differential system (3.14). 

So by the principle of mathematical induction, the scalar functions (3.19) are first integrals 
of the reducible system (3.14) for every natural number rn ^ 2. ■ 

Example 3.7. Consider the linear nonhomogeneous differential system 1 
^ = (<r + tp(t))x 2 + 5(asin(p(t) - ft cos <p{t))x 3 + f^t), 
dx 

—1= (cr + i;(t)) Xl + 6{a cos ^(t) + /? sin <p(t)) x 3 + / 2 (t), (3.21) 

dx 

—3 = j((3 cos ip(t) - asmip(t))x 1 - -f(acosip(t) + ft sin ip(t)) x 2 + / 3 (t), 

where ip: J — > R, /■ : J — > R, j = 1,2,3, are continuous functions on an interval Jcl, 
t 

the function (p: t f 4)(()d( for all t € J, t € J, and a, ft, j, 5, a ^ are real numbers. 

*o 

We claim that the differential system (3.21) is a reducible system by the orthogonal group 
of transformations 0(J, 3). Indeed, the orthogonal transformation 

y x = cos(^(t)x 1 — sin ip(t)x 2 , y 2 = sin(p(t)x 1 + cos ip(t)x 2 , y 3 = x 3 

reduce the system (3.21) to the system with constant coefficients 

=ay 2~ P S V3 + COSi P( t )fi( t ) ~ sin<p(t)/ 2 (t), 
= -o-y 1 + Q<5y 3 + sinv9(t)/ 1 (t) + cosv9(t)/ 2 (t), = ft^y x - a-yy 2 + f 3 (t). 



''Note that the linear homogeneous differential system corresponding to the nonhomogeneous system (3.21) 
is the linearization [8, pp. 139 - 142] of differential equations, which are describing the motion of a symmetric 
balanced nonautonomous gyrostat with one point of attachment [71, pp. 219 - 226; 72]. 
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The matrix C 



— aoo /?7 
a — cry 

- PS aS 



has the characteristic equation 



det(C - XE) = A (A 2 + (a 2 + (3 2 )j5 + a 2 ) = 0. 

Consider the real number D = (a 2 + (3 2 )^/5 + a 2 . We have three possible cases for building 
first integrals of the reducible linear differential system (3.21). 

Let D = 0. The matrix C has the eigenvalue X 1 = with the elementary divisor A 3 
corresponding to the real eigenvector v 10 = (cry, ^7, a) and to the generalized eigenvectors 

(7 7 \ / 27 27 2 \ 
(a + /3) — , (j3 — a) — , 1 1 , i/ 12 = I — ^ /?, =■ a, — ) . Using Theorems 3.5 and 3.6, we 
a cr / VcH a z a J 

can construct the basis of first integrals for the reducible differential system (3.21) 

F 1 : (t, x) — > j (a cos ip(t) + /3 sin (f(t)) x 1 + j(f3 cos (p(t) — a sin <p(t)) x 2 + ax 3 — C (t), (3.22) 



F 2 : (t,x) 



a 



((a + (3) cos c^(i) + ((3 — a) sin <^(t)) x x + 



+ ^ ((/?- a)cosc^(t) - (a + /?) sin y>(t)) x 2 + x 3 - tF^t^x) - C^t) for all (t,a;) eJxR 3 , 

27 27 
F 3 : (t, x) — > — 2 (/? cos </?(i) — a sin c^(i)) x 1 j ( a cos + P sm x 2 + 

2 

+ -x 3 -t 2 F 1 (t,x)-2tF 2 (t,x)-C 2 (t) for all (i,i)eJxR 3 , 

where the scalar functions 

t 

C :t^ J (7(acos^(C)+/3sin^(C))/ 1 (0 + 7(/3cos^(C)-asin^(C))/ 2 (C)+^/3(C))dC, 



1 

Ci= *->• 1 (^((a + ^cos^O + ^-aJsin^C))/!^) + 
+ 2( (/ 3_ a ) COSV3(C) _( a + / 3 )sin(/9(C )) /2(C) + / 3 ( C) + C (C))dC for all t G J, 

/(^ (/^MO " asiii^OjACC) - ^ (acos V (C) + /3sin ¥ »(C))/ 2 (C) + 

+ - / 3 (C) + 2^(0) d( for all t€J, t € J. 

Let D < 0. The matrix C has the linearly independent eigenvectors u w = (aj, /3j,cr), 
u 2 = (a 2 + 7 5/3 2 , -a^b-a^D ,5(-aa+(3\f^D)), u 3 = (ct 2 + 7 5/3 2 , -a^jS+ay/^D , 
— 8(cta + /3a/ — D )) corresponding to the eigenvalues A x = 0, A 2 = \/ — D, A 3 = — a/— -D- 
An integral basis of system (3.21) is the function (3.22) and the functions (by Theorem 3.5) 

F k : (t, x) -> ( ((ex 2 + 7J/3 2 ) cos <p(t) - (a/3 7< 5 + a\ k ) sin ^(t)) x x - 
- ((a/37<5 + crA fe ) cos <p(t) + (a 2 + 7<5/3 2 ) sin ip(t)) x 2 + 5(/3A fc - aa) x 3 ) e~ V - 
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t 

- J (((^ 2 + 7^ 2 ) cos <p(Q " («07<* + <?\) sm ?(()) A(C) - ((«07<* + aA fe ) cos y>(C) + 

+ (a 2 + 7 5/3 2 ) sin p(C)) / 2 (C) + <K/3A fc - twr)) / 3 (C)) e" A ^ V(t, x) G J x M 3 , fc = 2, 3. 

Let .D > 0. The matrix C has the linearly independent eigenvectors v 10 = (cry, (3j, <j), 
z^ 2 = ( - o" 2 - 7J/3 2 , a/3 7 5 + ctv 7 ^, *(a<7 - pVDi)), v 3 = (-a 2 - 7 5/3 2 , a/? 7 o" - av^Di, 

5{aa + flVDi)) corresponding to the eigenvalues A 1 =0, \ 2 =\fDi, and A 3 = — \fT)i. An 
integral basis of system (3.21) is the function (3.22) and the functions (by Corollary 3.1) 

t 

F k : (t,x 1 ,x 2 ,x 3 ) -»• a k (t,x 1 ,x 2 ,x 3 )- J a fc (C,/i(C)> /2(C). /3(C)) <*C> fc = 2,3, 

*o 

where t is a fixed point from the interval J, the scalar functions 

a 2 : (t,x) -> (-(a 2 + 7<5/3 2 )cosv A Dtcosv3(t) + (Q/375cos \fDt + a\/D sinv^Dt) sin</?(t)) x x + 
+ ( (a(3j5 cos v^D t + o\f~b ■sm\[Dt) cos^(t) + (<t 2 + 7<5/3 2 ) cos v^Dt sin </?(*)) x 2 + 
+ 5(aa cos v^D t - p^D sin v^D t) x 3 for all (t, x) G J x M 3 , 
a 3 : (t, x) ->■ ((cr 2 + 7<5/3 2 ) sin v^D t cos </>(t) + (av^D cos \f~D 1 - af3^8 sin \f~D t) sin <^(t)) x ± + 
+ ((crv^D cos v^Dt - a/?7o" sin \f~D t) cos tp(t) - (cr 2 + <yd/3 2 ) sin v^D t sin <p(i)) x 2 - 
- <5(/3\/Z> cos v^D t + aa sin v^D t) x 3 for all (t, x) G J x R 3 . 

Corollary 3.2. Lei t/ie differential system (3.14) be reducible to the system (3.15) by a 

transformation matrix g G G, and Zei A = A + Ai (Re A = A, ImA = A ^ 0) be the complex 
eigenvalue of the matrix C with elementary divisor of multiplicity m ^ 2 corresponding 
to an complex eigenvector is = u° +V°i (Rez/° = v°, Im^° = V°) and to generalized 
eigenvectors v k = v k + u k i (Rev k = v k , \m.v k = v k ), k = 1, . . . , m— 1. Then first integrals 
of the reducible system (3.14) are the functions 

k-i 

F e k+1 : (t, x) -> a, fc (t, x) - £ (*) t fe ~ T F, >T+1 (t , x) - (t) 

T = 

(3.23) 

/or a// (t, x) G J x R n , k = 1, . . . , m - 1, = 1, 2, 
where the first integrals (by Corollary 3.1) 

F : (t,x) -»• a gQ (t,x) - C gQ (t) for all (t, x) G Jx R n , 9 = 1,2, 
the scalar functions 

a lk : (t,x) ->■ (i> k g{t)x cos At + u k g(t)x sin At) exp( - At) /or a// (t,x) G J x E n , 

a 2fc : (t,x) ->■ (z? fc g(t)x cos At - v k g{t)x sin At) exp( - At) V(t,x) G J x R n , 
t 

C ek .t^ j(a ek {Q,f{0) + kC ek _ 1 {Q)dC for all t G J, 9 = 1, 2, A; = 0, . . . , m - 1, t G J. 

*o 
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Proof. Since the functions (3.19) are complex- valued first integrals (by Theorem 3.6) of 
the reducible system (3.14), we see that the real and imaginary parts of the functions (3.19) 
are the real first integrals (3.23) of the reducible system (3.14). ■ 

Example 3.8. The linear nonhomogeneous Euler differential system 



dx 1 
~~dt 



x 2 + a, 



dx 2 
~~dt 



b , 

X 3 + ^2 lnt 



dxn c 

—rr = x 4 + 77 sm hi *> 
at t d 



dx A 1 

It = ~¥ Xl 



3 9 6 d 



is reducible to the linear differential system with constant coefficients 



dVi 
dr 



y 2 + ae T , 



dy 2 

dr 



= V2 + V3 + bT ' 



dy 3 

dr 



= 2y 3 + y A + csinr, 



(3.24) 



dr 



d 



COST 



by the exponential group Exp(4) with the exponential transformation 
Dl = x li V2 = ^ x 2' ^3 = ^ x 3> Va = ^ x a-> t = e T 



The matrix C 












- 1 


1 


1 





- 3 





1 


2 


- 9 








1 


- 3 



has the complex eigenvalues A 1 



— i and A 2 = i 



with the multiple elementary divisors (A + i) 2 and (A — i) 2 , respectively. 

The eigenvalue X 1 = — i of the matrix C is corresponding to the complex eigenvector 
u°= (1, l + 2i, l + 3i, i) and to the generalized eigenvector of the 1-st order u 1 = (—i, 2+i, i, 0). 

By Corollaries 3.1 and 3.2, using the numbers Ai = 0, A x = — 1, the linearly independent 
vectors ^ = (1,1,1,0), v° = (0, 2, 3, 1), ^ = (0,2,0,0), i/ 1 = ( — 1,1,1,0), the functions 

a 1Q : (t, x) — > cos In t {x l + tx 2 + t 2 x 3 ) — t sin In t (2x 2 + 3tx 3 + t 2 x A ) for all (t, x) G J l x IR 4 , 

a 2Q : (t,x) — > t cos lnt (2x 2 + 3tx 3 + t 2 x 4 ) + sin lnt (x 1 + tx 2 + t 2 x 3 ) for all (t, x) G J t x IR 4 , 

a n : (t,x) — > 2tcoslntx 2 + sin lnt (x l — tx 2 — t 2 x 3 ) for all (t,x) £ J,x IR 4 , 

q 21 : (t, x) — >■ cos In t ( — x 1 + tx 2 + t 2 x 3 ) + 2t sin lnt x 2 for all (t, 1) £ J ( x IR 4 , 



and 



C, n (t) = / ( (a+- lnt + - sin lnt) cos lnt— lnt+^ sinlnt + ^ — — — )tsi 
iU J \ \ t t / V t z t 2 t 2 cos In t / 



sin In t ) dt 



- (cos In t + sin In t) t + 6(cos In t — 2 sin In t + (2 cos In t + sin In t) In t) 
— — (6 In t — 2 sin 2 In t — 3 sin 2 In t) + d In | cos In t| for all t G , 



C 20 (t) 



/26 3c . , d 

-77 In t + -77 sm m t + -= : — 

Vt 2 t 2 t 2 cos lnt 



^ t cos In t + (a + - In t + - sin hi ' 



^ sin In t^ 



(it 



(sin lnt — cos In t)t + 6(2 cos lnt + sin lnt + (2 sin lnt — cos lnt) lnt) + 



+ - (2 lnt + 6sin 2 lnt - sin21nt) + dint for all t € J t , 
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C u (t) = J In t cos In t + (a — j In t — j sin In t^j sin In t + C w (t)^j dt = 

= ^ (2 sin In t — cos In t) t + 6 (6 cos In t + sin In t + 4 In t sin In t) — 
- | (21nt + 61n 2 t + 3cos21nt - sin21nt) + d(\nt In | cosing + /3(lni)) for all t G J t , 



C„ (t) 



J ^ — a + ^ In i + ^ sin In tj cos In t + ^ In i sin In t + C 20 (t)^j dt = 



= — ^ (2 cos In t + sin In t ) t — b (cos In t — 6 sin In t + 4 In t cos In t) + 

c d 
+ -(61nt + 21n 2 t + cos21nt-3sin21nt + 4sin 2 lnt) + - In 2 i for all f£ J„ 
8 2 

where f3(t) = J tt&ntdt, we can build the first integrals of the reducible system (3.24) 
F n : (t,x) — > cos In t (x l + tx 2 + t 2 x 3 ) — t sin In t (2x 2 + 3tx 3 + t 2 x A ) — 

— ^ (cos In t + sinlni)i — 6(coslnt — 2 sin In t + (2 cos \nt + sin In t) hit) + 

+ |(61ni-2sin 2 lni-3sin21ni)-(iln|coslni| for all (i,x)e J ( xl 4 , 

i ? 21 : (i,x) — > t coslni (2x 2 + 3ix 3 + t 2 x A ) + sin In t (x x + tx 2 + t 2 x 3 ) — 

— ^ (sin In t — coslnt)i — 6(2 cos In t + sin In t + (2 sin In t — cos In t) hit) — 

- |(21nt + 6sin 2 lnt-sin21nt) -dint for all (t, x) G J, x M 4 , 

F 12 : (i , x) — )■ 2t cos In t x 2 + sin In t (x l — tx 2 — t 2 x 3 ) — lntF n (t,x) — 

— — (2 sin In t — cos In t) t — 6(6 cos In t + sin In t + 4 In t sin In t) + 

+ - (2 In t + 6 In 2 t + 3 cos 2 In t - sin 2 In i ) - d (in i In | cos In t | + f3 (In t )) for all (t , x) G J l ; 
8 

F 22 : (i, x) — )■ cos In £ ( — x : + tx 2 + t 2 x 3 ) + 2t sin In t x 2 — In t F 21 (t,x) + 

+ — (2 cos lnt + sinlnt)t + 6(cos hit — 6 sin In t + 4 In t cos In t) — 

c d 

(6 In t + 2 In 2 1 + cos 2 In t - 3 sin 2 In t + 4 sin 2 In i) In 2 1 for all (t, i) £ J, x 

8 2 

on any domain J, x I 4 C K 5 , I = 0,1,2, ... , where the intervals 



(0;e 2 ), J s =(e 2+W{S l) ;e 2+WS ), s = l,2,.... 



The functionally independent first integrals F u , F 21 , F 12 , and F 22 are an integral basis 
of the Euler differential system (3.24) on any domain J t x M 4 . 
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